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ABSTRACT 

Heat  transfer  performance  of  horizontal  roped  tubes  was 
studied  using  an  evaporator/condenser  test  rig.  The  condenser 
section  contained  a  single  horizontal  tube  and  was  operated  at 
both  vacuum  and  atmospheric  conditions.  A  vacuum  pump  was 
integral  to  the  system  to  ensure  the  presence  of 
noncondensible  gases  did  not  affect  the  results.  The  test 
tubes  were  Wolverine  MHT  Korodense  titanium  tubes  with  an 
outside  diameter  of  16.07  mm  and  a  corrugation  pitch  of  8  mm. 
The  effect  of  wire-wrapping  was  studied  using  titanium  wires 
ranging  from  0.1  to  1.6  mm  in  diameter  corresponding  to  a 
range  of  pitch  to  diameter  ratio  from  80  to  5.  The  outside 
heat  transfer  coefficient,  hQ,  was  determined  by  first 
experimentally  obtaining  the  overall  heat  transfer 
coefficient,  UQ,  and  then  utilizing  a  modified  Wilson  Plot 
procedure . 

Results  obtained  for  steam  condensation  indicate  that  the 
MHT  Korodense  titanium  tube  provides  negligible  enhancement 
compared  to  a  titanium  tube.  The  addition  of  wire -wrapping 
provided  minimal  increase  in  the  outside  heat  transfer 
coefficient,  about  8%  at  a  fractional  coverage  of 
approximately  5%  corresponding  to  a  pitch  to  wire  diameter 
ratio  near  20. 


in 


// 

TABLE  OF  CONTENTS 

I.  INTRODUCTION   1 

A.  BACKGROUND 1 

B.  CONDENSATION  APPLICATIONS  FOR  ENHANCED  TUBES   .  3 

C.  NAVAL  POSTGRADUATE  SCHOOL  CONDENSATION  RESEARCH  5 

D.  OBJECTIVES 6 

II.  LITERATURE  SURVEY   7 

A.  INTRODUCTION   7 

B.  VAPOR  SIDE  CONSIDERATIONS  8 

1.  Horizontal  Smooth  Tubes   8 

2..  Wire-Wrapped  Smooth  Tubes 10 

3 .  Corrugated  Tubes  12 

4.  Wire-Wrapped  Corrugated  Tubes   16 

C.  COOLANT  SIDE  CONSIDERATIONS  17 

1.  Smooth  Tubes 17 

2 .  Corrugated  Tubes  19 

3.  Tubes  With  Inserts 23 

III.  EXPERIMENTAL  APPARATUS   25 

A.  SYSTEM  OVERVIEW 2  5 

B.  SYSTEM  INSTRUMENTATION   3  0 

C.  TUBES  TESTED 31 


IV 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY  CA    93943-5101 

IV.  EXPERIMENTAL  PROCEDURES  AND  DATA  ANALYSIS   ....  34 

A.  TUBE  PREPARATION 3  4 

B.  SYSTEM  START-UP  AND  SHUTDOWN  PROCEDURES  ....  35 

C.  DATA  COLLECTION  PROCEDURE 3  9 

D.  DATA  REDUCTION  PROCEDURE 3  9 

E.  MODIFIED  WILSON  PLOT  TECHNIQUE   43 

F.  ENHANCEMENT  RATIO  4  6 

V.  RESULTS  AND  DISCUSSION 48 

A.  HEAT  TRANSFER  CORRELATIONS 4  8 

B.  VERIFICATION  AND  REPEATABILITY  OF  RESULTS  ...  49 

C.  ANALYSIS  OF  CORRUGATED  TUBE  RESULTS 53 

D.  ANALYSIS    OF    WIRE -WRAPPED    CORRUGATED    TUBE 
RESULTS 5  9 

VI.  CONCLUSIONS  AND  RECOMMENDATIONS   79 

A.  CONCLUSIONS 79 

B.  RECOMMENDATIONS 79 

APPENDIX  A.   SYSTEM  CORRECTIONS  81 

A.   FRICTIONAL  TEMPERATURE  CORRECTIONS   81 

APPENDIX  B.   UNCERTAINTY  ANALYSIS  84 

APPENDIX  C.   DRPML  PROGRAM  LISTING   9  0 


v 


LIST  OF  REFERENCES 122 


INITIAL  DISTRIBUTION  LIST 129 


VI 


Table  I. 
Table  II. 

Table  III. 
Table  IV. 

Table  V. 
Table  VI. 

Table  VII. 

Table  VIII 
Table  IX. 
Table  A.  1 


LIST  OF  TABLES 

LISTING  OF  TUBES  TESTED   3  3 

SMOOTH  TITANIUM  TUBE  AT  ATMOSPHERIC 

PRESSURE 54 

SMOOTH  TITANIUM  TUBE  AT  VACUUM  PRESSURE   .    5  5 
0.5  mm  WIRE-WRAPPED  SMOOTH  TITANIUM  TUBE 

WITH  HEATEX  INSERT 5  6 

PLAIN  LPD  KORODENSE  TUBE 5  7 

0.5  mm  WIRE -WRAPPED  LPD  KORODENSE  TUBE  WITH 

HEATEX  INSERT   5  8 

0.5  mm  WIRE -WRAPPED  SMOOTH  COPPER  TUBE  WITH 

HEATEX  INSERT   5  9 

PLAIN  MHT  KORODENSE  TUBE 6  0 

WIRE -WRAPPED  AND  PLAIN  MHT  KORODENSE  TUBES    6  8 
FRICTIONAL  TEMPERATURE  RISE  EQUATIONS  FOR  MHT 
KORODENSE  TUBE 8  3 


VII 


LIST  OF  FIGURES 

Figure  1.    Profile  of  a  Corrugated  Tube 3 

Figure  2.    Schematic  of  the  Single  Tube  Test  Apparatus    26 

Figure  3.    Schematic  of  the  Test  Section 27 

Figure  4.    Schematic  of  the  Purging  System  and  Cooling 

Water 29 

Figure  5.    hQ  vs .  ATcf  for  a  Plain  MHT  Korodense  Tube 

Under  Vacuum 5  0 

Figure  6.    hQ  vs  .  ATcf  for  a  Smooth  Titanium  Tube  Under 

Vacuum 51 

Figure  7.    hQ  vs .  ATcf  for  a  Smooth  Titanium  Tube  at 

Atmospheric  Pressure   52 

Figure  8.    hQ  vs .  ATcf  for  a  Smooth  Titanium  and  a  Plain 

MHT  Korodense  Tubes  Under  Vacuum   61 

Figure  9.    UQ  vs .  Vw  for  a  Plain  MHT  Korodense  Tube  at 

Atmospheric  Pressure   62 

Figure  10.   UQ  vs.  Vw  for  a  Plain  MHT  Korodense  Tube  Under 

Vacuum 63 

Figure  11.   UQ  vs.  Vw  for  Wire-Wrapped  MHT  Korodense  Tubes 

Under  Vacuum  with  HEATEX  Insert 64 

Figure  12.   hQ  vs.  ATcf  for  Wire-Wrapped  MHT  Korodense 

Tubes  Under  Vacuum 65 

Figure  13.   hQ  vs.  ATcf  for  Wire-Wrapped  MHT  Korodense 

Tubes  Under  Vacuum 66 


Vlll 


Figure  14.   Enhancement  Ratio  vs.  Pw/Dw  for  Wire-Wrapped 

Korodense  Tubes  69 

Figure  15.   Enhancement  Ratio  vs.  Fractional  Coverage  for 

Wire -Wrapped  Korodense  Tubes   70 

Figure  16.  Enhancement  Ratio  vs.  Pw/Dw  for  Wire-Wrapped 
Korodense  and  Wire-Wrapped  Smooth  Tubes,  NPS 
Data 72 

Figure  17.  Enhancement  Ratio  vs.  Fractional  Coverage  for 
Wire-Wrapped  Korodense  and  Wire-Wrapped  Smooth 
Tubes,  NPS  Data 73 

Figure  18.   Enhancement  Ratio  vs.  Pw/Dw  for  Wire-Wrapped 

Korodense  and  Wire -Wrapped  Smooth  Tubes  .  .    75 

Figure  19.  Enhancement  Ratio  vs.  Fractional  Coverage  for 
Wire-Wrapped  Korodense  and  Wire-Wrapped  Smooth 
Tubes 7  6 

Figure  20.  Ci   vs.  Fractional  Coverage  for  Wire-Wrapped 
Korodense  and  Wire-Wrapped  Smooth  Tubes  at 
Atmospheric  Pressure   78 

Figure  A.l   Frictional  Temperature  Rise  for  the  MHT 

Korodense  Titanium  Tube  with  and  without  a 
HEATEX  Insert  82 


IX 


2 


NOMENCLATURE 

a  as  defined  in  equation  (4.30) 

A^  effective  inside  surface  area,  m2 

AQ  effective  outside  surface  area,  m 

b  as  defined  in  equation  (4.26) 

c  as  defined  in  equation  (2.22) 

c  specific  heat  at  a  constant  pressure,  J/kg-K 

c8  as  defined  in  equation  (2.14) 

C1  leading  coefficient  for  the  inside  heat  transfer 

correlation  used 

Cn  as  defined  in  equation  (2.8) 

d  as  defined  in  equation  (2.23) 

dQ  external  groove  diameter,  m 

D^  tube  inside  diameter,  m 

DQ  tube  outside  diameter,  m 

DTI  wire  diameter,  m 

e  groove  depth,  m 

E1  as  defined  in  equation  (2.32) 

E2  as  defined  in  equation  (2.33) 

E3  as  defined  in  equation  (2.35) 

f  as  defined  in  equation  (2.30) 

F  as  defined  in  equation  (2.4) 

g  gravitational  constant,  9.81  m/s 

hf  specific  enthalpy  of  vaporization,  J/kg 


x 


hi  inside  heat  transfer  coefficient,  W/m2-K 

hQ  outside  heat  transfer  coefficient,  W/m^-K 

hoe  outside  heat  transfer  coefficient  for  enhanced  tubes, 

W/m2-K 

hos  outside  heat  transfer  coefficient  for  smooth  tubes, 

W/m2-K 

kc  thermal  conductivity  of  the  coolant  film,  W/m-K 

kf  thermal  conductivity  of  the  condensate  film,  W/m-K 

k^  thermal  conductivity  of  the  tube  material,  W/m-K 

K1  as  defined  in  equation  (2.26) 

K2  as  defined  in  equation  (2.27) 

L  active  length  of  tube  exposed  to  steam,  m 

LQ  external  groove  depth,  m 

L1  length  of  inlet  portion  of  tube,  m 

L2  length  of  outlet  portion  of  tube,  m 

LMTD  log  mean  temperature  difference,  K 

m  Reynolds  number  exponent 

m  mass  flow  rate  of  the  coolant,  kg/s 

m1  as  defined  in  equation  (4.25) 

n  Prandtl  number  exponent 

nt  number  of  tubes  in  a  vertical  row 

Nu  Nusselt  number 

P  corrugation  pitch,  m 

Pw  wire-wrap  pitch,  m 

Psat  saturation  pressure,  Pa 

Pr  Prandtl  number 


XI 


Prw  Prandtl  number  at  the  wall  temperature 

PW  as  defined  in  equation  (2.10) 

Q  heat  transfer  rate,  W 

q  heat  flux,  W/m2 

r  tube  specific  constant  as  defined  for  equation  (2.30) 

R  radius  of  curvature  of  corrugation  ridge,  m 

Rf  fouling  resistance,  K/W 

Ri  inside  thermal  resistance,  K/W 

RQ  outside  thermal  resistance,  K/W 

Rtotal  tot:al  thermal  resistance,  K/W 

Rw  wall  thermal  resistance,  K/W 

Re  Reynolds  number 

Ref  Reynolds  number  for  the  condensate  film 

Re2A  two  phase  Reynolds  number 

Tcf  temperature  difference  across  the  condensate  film,  K 

ATcf  temperature  difference  across  the  film,  K 

Tm  mean  coolant  bulk  temperature,  K 

Tsat  vapor  saturation  temperature,  K 

Two  outside  wall  temperature,  K 

Tx  cooling  water  inlet  temperature,  K 

T2  cooling  water  outlet  temperature,  K 

U0  overall  heat  transfer  coefficient,  W/m2-K 

Um  vapor  velocity,  m/s 

Vw  coolant  velocity,  m/s 

x  Reynolds  number  exponent  in  equation  (4.15) 

X  as  defined  in  equation  (4.24) 


XII 


Y     as  defined  in  equation  (4.23) 

Z      as  defined  in  equation  (4.19  or  4.20) 

Z*     as  defined  in  equation  (2.9) 

Greek  Symbols 

oi  leading  coefficient  for  the  outside  heat  transfer 

correlation  used 
6i  helix  angle  measured  from  the  tube  axis 

0     as  defined  in  equation  (2.37) 
e  as  defined  in  equation  (2.25) 

eT     enhancement  ratio  based  on  AT  for  wire-wrapped  smooth 

tubes 
eTc    enhancement   ratio   based   on  AT   for   wire -wrapped 

corrugated  tubes 
7     as  defined  in  equation  (2.29) 
r     mass  flow  rate  of  condensate  per  unit  length  of 

horizontal  condenser  tube 
juc     dynamic  viscosity  of  the  coolant,  kg/m-s 
lif  dynamic  viscosity  of  the  condensate  film,  kg/m-s 

/iw     dynamic  viscosity  of  condensate  at  the  wall,  kg/m-s 
rj  surface  efficiency 

*c    as  defined  in  equation  (2.14) 

co     tube  specific  constant  as  defined  for  equation  (2.10) 
Q     as  defined  in  equation  (4.15,  4.16  or  4.17) 
pf    density  of  the  condensate  film,  kg/m3 
pv     density  of  the  vapor,  kg/m3 

xiii 


crf     surface  tension  of  the  condensate  at  film  temperature, 

lbf/ft 
0  severity  factor 

d*  tube  specific  constant  as  defined  for  equation  (2.9) 


xiv 


ACKNOWLEDGEMENTS 

The  author  would  like  to  express  his  sincere  appreciation 
to  Professor  Paul  J.  Marto  and  Professor  Stephen  B.  Memory  for 
their  advice,  guidance,  and  continual  support  towards  the 
completion  of  this  thesis. 

Many  thanks  is  also  given  to  Mr.  Tom  McCord  and  the 
workers  in  the  Mechanical  Engineering  Department  Machine  Shop 
for  their  timely  fabrication  of  all  the  required  materials 
necessary  to  complete  this  endeavor. 

The  author  wishes  to  extend  his  greatest  appreciation  and 
thanks  to  his  wife,  Laura.  Without  her  extreme  patience, 
love,  and  encouragement  this  project  would  not  have  been 
possible . 


xv 


INTRODUCTION 


A .   BACKGROUND 

Reductions  in  size,  weight  and  cost  have  become 
increasingly  important  design  considerations  in  the 
development  of  every  component  destined  for  use  aboard  United 
States  Navy  vessels.  Significant  reductions  in  weight  have 
already  been  achieved  for  main  propulsion  machinery  with  the 
advent  of  the  gas  turbine  propulsion  plant.  Weight  reduction 
is  also  possible  through  redesign  of  auxiliary  systems  and 
other  large  components  such  as  main  steam  condensers.  The  use 
of  heat  transfer  enhancement  techniques  for  condenser  tubing 
has  been  one  avenue  receiving  considerable  attention  in  the 
pursuit  of  this  goal.  Smaller  more  efficient  condensers  would 
result  in  reductions  in  capital  cost,  space  requirements, 
weight  and  operating  cost. 

The  four  primary  factors  that  affect  the  performance  of 
condensers  are  the  thermal  resistances  on  the  coolant  side, 
across  the  tube  wall,  the  vapor  side,  and  also  that  due  to 
fouling.  Each  can  be  reduced  to  minimize  its  effect,  which  in 
turn  enhances  the  overall  heat  transfer  rate  of  the  system. 
In  experimental  studies,  the  fouling  resistance  can  be 
neglected  if  the  tube  is  thoroughly  clean  prior  to  testing. 


Enhancement  techniques  can  be  categorized  into  two  major 
groups,  active  and  passive.  Guido  [Ref .  1]  compiled  a  fairly- 
comprehensive  listing  of  both  groups  of  enhancement  techniques 
for  the  tube  side  heat  transfer  rate.  Passive  techniques  such 
as  internal  helical  ribbing  and  displaced  promoters  (the  use 
of  wire  mesh  or  HEATEX  inserts)  will  be  discussed  in  detail  in 
a  later  section.  It  should  be  noted,  however,  that  both  of 
these  techniques  have  a  corresponding  increase  in  pumping 
power  requirements.  Wall  thermal  resistance  can  be  reduced 
passively  by  decreasing  the  wall  thickness  or  changing  the 
tube  material.  Similar  to  enhancement  techniques  for  the  tube 
side,  the  shell  side  resistance  can  also  be  reduced  by  both 
active  and  passive  means.  The  passive  enhancement  techniques 
of  wire-wrapping  and  tube  corrugation  will  also  be  covered  in 
a  later  section. 

A  variety  of  corrugated  or  roped  tubes  are  available  for 
use  in  condenser  applications.  Figure  1  shows  the  profile  of 
one  type  of  corrugated  tube  which  has  a  groove  radius  smaller 
than  the  ridge  radius .  Wolverine  Tube  Company  manufactures  a 
low  pressure  drop  (LPD)  and  a  maximum  heat  transfer  (MHT)  tube 
of  this  type  (trade  name:  Korodense) ,  available  for  commercial 
use.  The  MHT  tubes  have  a  larger  groove  depth  compared  to  the 
LPD  tubes  [Ref.  2] .  Other  types  of  corrugated  tubes  are  also 
available  which  have  a  larger  groove  radius  than  the  ridge 
radius.  Marto  ec  al .  [Ref.  3]  present  test  results  on  various 
other  corrugation  geometries. 


Figure  1.   Profile  of  a  corrugated  tube 


B.   CONDENSATION  APPLICATIONS  POR  ENHANCED  TUBES 

In  1971  Withers  and  Young  [Ref .  4]  performed  a  theoretical 
analysis  on  a  power  plant  main  condenser.  Their  findings 
indicated  a  substantial  saving  in  almost  every  condenser 
aspect  if  the  smooth  tubes  were  replaced  by  Korodense  roped 
tubes.  Similar  results  were  obtained  by  Webb  et  al .  [Ref.  5]  . 
They  indicated  a  possible  30  to  40  percent  increase  in  the 
overall  heat  transfer  coefficient  when  analyzing  a  tube-for- 
tube  replacement  of  smooth  to  Korodense  tubes.  This  benefit 
was  calculated  for  two  separate  power  plants  of  different  heat 
duty.  An  ongoing  comprehensive  study  of  14  Tennessee  Valley 


Authority  (TVA)  power  plant  condensers,  retubed  with 
corrugated  tubes,  has  yielded  significant  economic  and 
performance  benefits  [Refs.  6  and  7]  .  Of  the  14  condensers 
retubed  by  TVA,  their  Gallatin  Unit  1  has  provided  long-term 
results  with  over  12  years  of  successful  operation  since 
retubing.  In  1983  Boyd  et  al .  [Ref .  8]  studied  the 
performance  of  TVA's  Gallatin  Unit  1  plant  after  only  three 
years  of  operation.  Their  findings  yielded  up  to  a  43% 
increase  in  the  overall  heat  transfer  coefficient  when  clean 
tubes  were  analyzed.  However,  data  also  indicated  that  the 
corrugated  tubes  fouled  quicker  than  plain  tubes.  Rabas  et 
al .  [Ref.  9]  studied  12  of  the  retubed  TVA  condensers  and 
found  fouling  rates  for  the  enhanced  tubes  ranging  from  about 
the  same  as  to  almost  twice  that  experienced  in  plain  tubes. 
Although  the  fouling  rates  were  higher,  Rabas  et  al .  [Ref.  9] 
reported  the  fouling  resistance  values  for  the  LPD  Korodense 
tubes  never  exceeded  the  minimum  allowable  TEMA  standard. 

Union  Carbide  [Ref.  10]  performed  analysis  on  roped 
condenser  tubes  with  a  Thin  Film  Condensing  Promoter 
(T.F.C.P.),  which  is  essentially  a  wire-wrapped  roped  tube. 
Their  findings  indicated  an  increase  in  the  heat  transfer 
coefficient  by  a  factor  of  between  2  to  3  times  that  found 
with  a  smooth  tube  under  condensate  loaded  conditions. 


C.   NAVAL  POSTGRADUATE  SCHOOL  CONDENSATION  RESEARCH 

Titanium  has  a  number  of  properties  which  make  it  a 
desirable  alternative  to  the  conventional  copper-nickel  tubing 
currently  in  use  today  on  board  Naval  vessels.  Its  high 
strength  to  weight  ratio  enables  the  tubes  to  be  designed  with 
a  much  thinner  wall  thickness,  while  still  maintaining  the 
durability  of  the  overall  system.  Titanium  also  has  a  high 
resistance  to  erosion  and  sea  water  corrosion.  This  leads  to 
reduced  shutdown  time  for  tube  cleaning  and  replacement.  One 
of  the  primary  disadvantages  of  titanium  tubing  is  its  poor 
thermal  conductivity.  Titanium  has  approximately  50%  the 
thermal  conductivity  of  copper-nickel  [Ref .  11] .  This 
disadvantage  can  be  partially  offset  by  a  decreased  wall 
thickness  and  also  through  the  use  of  enhanced  surfaces.  The 
Naval  Postgraduate  School  (NPS)  ,  in  conjunction  with  the  Naval 
Surface  Warfare  Center,  Carderock  Division,  Annapolis 
Detachment,  has  been  testing  the  enhancement  characteristics 
of  various  types  of  titanium  tubing. 

In  order  to  examine  the  heat  transfer  characteristics  of 
various  types  of  enhanced  tubes,  a  suitable  test  platform  must 
be  utilized.  While  the  ideal  test  platform  would  be  a  full 
scale  model  of  the  actual  condenser  used  in  the  field,  the 
cost  associated  with  such  an  apparatus  would  be  prohibitive. 
The  condensation  test  platform,  constructed  by  Krohn  [Ref.  12] 
at  NPS,  has  been  in  use  since  1982.  Although  several 
modifications  have  been  made  by  various  researchers  since  that 


time,  the  overall  apparatus  and  the  operation  remain  basically 
unchanged.  Van  Petten  [Ref .  13]  lists  a  comprehensive  summary 
of  the  work  completed  by  eight  NPS  researchers  from  1982 
through  1987.  Most  of  the  work  concentrated  on  the 
performance  of  various  types  of  externally  finned  tubes. 
Mitrou  [Ref.  14],  in  addition  to  studying  finned  tubes,  also 
performed  experiments  on  wire -wrapped  smooth  copper  tubes. 
More  recent  work  conducted  by  Coumes  [Ref.  15],  Guttendorf 
[Ref.  16]  and  Swensen  [Ref.  17]  dealt  with  further  studies  of 
finned  copper  tubes.  O'Keefe  [Ref.  18]  studied  the 
performance  of  various  titanium  tubes.  His  tests  included 
smooth  tubes,  wire -wrapped  smooth  tubes,  LPD  Korodense  tubes, 
and  wire-wrapped  LPD  Korodense  tubes.  Other  researchers 
including  Kanakis,  Brower  and  Mazzone  [Refs.  19-21]  have 
studied  tube  bundle  configurations  for  smooth  and  corrugated 
tubes,  both  with  and  without  wire -wrapping. 

D .   OBJECTIVES 

The  main  objectives  of  this  study  are  as  follows: 

1.  To  verify  the  repeatability  of  the  calculated 
coefficients  for  both  coolant  side  and  vapor  side  heat 
transfer  correlations  obtained  by  previous  researchers  on 
wire-wrapped  smooth  and  corrugated  titanium  tubes. 

2.  To  manufacture  and  collect  condensation  data  on  smooth 
and  MHT  Korodense  titanium  tubes  with  varying  fractional 
coverage  due  to  wire -wrapping. 

3 .  To  determine  if  any  relationship  exists  for  optimal 
fractional  coverage  or  pitch  to  diameter  ratio  between 
wire -wrapped  smooth  tubes  and  wire -wrapped  MHT  and  LPD 
Korodense  tubes. 


II.   LITERATURE  SURVEY 


A.   INTRODUCTION 

Condensation  on  horizontal  tubes  has  been  an  intensely 
studied  heat  transfer  mechanism  throughout  the  years.  New 
developments  in  condenser  tubing  require  continued 
experimental  and  theoretical  research  in  this  field. 

There  are  two  fundamental  modes  of  condensation  on 
horizontally  mounted  condenser  tubes,  filmwise  and  dropwise. 
Dropwise  condensation  can  produce  heat  transfer  rates  as  much 
as  an  order  of  magnitude  larger  than  those  associated  with 
filmwise  condensation.  Various  surface  coatings  have  been 
developed  in  order  to  promote  dropwise  condensation  and  take 
advantage  of  the  higher  heat  transfer  rates.  However,  with 
time  these  coatings  lose  their  effectiveness  and  filmwise 
condensation  eventually  returns.  For  this  reason,  condenser 
design  calculations  are  normally  based  on  the  more 
conservative  filmwise  mode  of  condensation  [Ref .  22] . 

Internally,  the  primary  factors  affecting  thermal 
resistance  are  the  intensity  of  turbulence  and  the  thermal 
conductivity  of  the  coolant  medium.  The  degree  of  turbulence 
can  be  changed  by  altering  the  internal  tube  wall  geometry  or 
through  the  use  of  displacement  promoters. 


B.   VAPOR  SIDE  CONSIDERATIONS 

Condensation  on  horizontal  tubes  consists  of  a  complex  two 
phase  process.  A  suitable  theoretical  model  must  be  utilized 
in  order  to  accurately  predict  these  phenomena. 

1.   Horizontal  Smooth  Tubes 

Nusselt  [Ref.  23]  developed  the  foundation  for  the 
study  of  filmwise  condensation  on  horizontal  tubes  in  1916. 
His  work  was  formulated  for  a  quiescent  vapor  condensing  on  a 
single  horizontal  tube.  Due  to  the  increase  in  the  thickness 
of  condensate  as  gravity  forces  it  around  the  sides  of  the 
tube,  the  local  heat  transfer  coefficient  decreases  around  the 
tube  circumference.  Nusselt  developed  a  theory  for  the 
average  heat  transfer  coefficient  around  the  tube,  accounting 
for  the  lower  resistance  at  the  top  of  the  tube  where  the  film 
was  the  thinnest  and  the  higher  resistance  at  the  bottom  of 
the  tube  where  the  film  was  the  thickest.  As  a  result,  he 
developed  the  following  equation: 


hQ  =   0.72: 


k£  g  pf  (pf-py)  hfg 
»f  Do   (Tsat-Two) 


U/4 


(2.1) 


where   the   fluid  properties   are   evaluated  at   the   film 
temperature,  given  by: 

Tf   =  —  T    +  —  T  (2.2) 

f  r>      sat  ~>      wo 

Nusselt 's  theory  has  been  extensively  studied  since 
1916,  and  with  the  imposed  assumptions,  has  been  found  to  be 


generally  valid  [Refs.  24,  25] .  It  has  also  been  found  to  be 
accurate  for  cases  which  do  not  conform  to  Nusselt's  original 
assumptions,  such  as  variable  wall  temperature  [Ref .  26]  .  One 
of  the  major  problems  encountered  in  applying  Nusselt's  theory 
in  the  design  of  condensers  arises  from  his  assumption  of  a 
quiescent  vapor.  While  in  theory,  and  in  some  limited 
practical  applications,  the  assumption  of  a  stationary  vapor 
can  be  justified,  most  steam  condensers  operate  under 
conditions  where  the  vapor  is  travelling  at  some  constant 
velocity . 

With  the  vapor  in  motion,  sufficient  shear  forces  can 
be  developed  which  significantly  thin  the  condensate  film  and 
increase  the  outside  heat  transfer  coefficient.  Early 
theoretical  work  by  Shekriladze  and  Gomelauri  [Ref.  27]  ,  which 
accounted  for  vapor  shear,  resulted  in  the  following  equation 
to  approximate  the  mean  Nusselt  number. 

=  0.64(1  +(1  +  1.69F)1/2)1/2  (2.3) 


Dp   1/2 
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where : 


F  -  9D°»'  h"  (2.4) 


V-    kf    (Tsac-Tu 
and  the  two  phase  Reynolds  number  is  given  by 


Re      -    0'  "•  D°  (2.5) 


Lee  and  Rose  [Ref.  28]  compared  vapor  shear  models  with 
experimental  results  and  found  that  the  Shekriladze-Gomelauri 
results  were  more  conservative  than  more  rigorous  developments 
performed  by  other  researchers  due  to  their  simplified 
approximation  for  the  interfacial  shear  stress.  Fujii  et  al . 
[Ref.  29]  ,  in  a  more  recent  study,  performed  an  extensive 
development  on  the  outside  of  a  horizontal  tube.  From 
experimental  data,  they  also  developed  an  empirical 
formulation  for  the  condensation  of  steam  on  a  horizontal  tube 
which  includes  vapor  velocity  effects: 

=  0.96  F1/5  (2.6) 


Re2   ] 


where  F  and  Re2(p  are  defined  in  equations  (2. 4)  and  (2.5)  .  For 
situations  where  the  surface  shear  forces  dominate,  equation 
(2.6)  tends  to  more  accurately  predict  the  vapor  side  heat 
transfer  coefficient  for  steam. 

As  discussed  previously,  the  primary  modes  of  external 
tube  enhancement  covered  in  this  study  are  wire -wrapping  and 
the  addition  of  helical  grooves.  Both  methods  have  associated 
advantages  and  disadvantages. 

2.   Wire-Wrapped  Smooth  Tubes 

The  benefits  of  wire-wrapping  smooth  tubes  have  been 
attributed  to  two  primary  effects:  the  thinning  of  the 
condensate  film  between  the  wires  due  to  surface  tension 
forces  and  the  better  drainage  of  the  condensate  along  the 
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wire,  (Fujii  et  al .  [Ref.  30]  and  Marto  et  al .  [Ref .  31]). 
Thomas  [Ref.  32]  conducted  tests  on  vertical  tubes  with  wires 
attached  vertically  along  the  tubes  and  demonstrated  this 
latter  enhancement  effect.  Horizontal  smooth  tubes  with 
coiled  wire  were  later  studied  by  Thomas  et  al .  [Ref.  33]  and 
their  findings  indicated  a  two  fold  increase  in  the  heat 
transfer  coefficient  over  that  predicted  by  Nusselt  theory  for 
a  smooth  tube.  Marto  et  al .  [Ref.  31]  determined  an  optimum 
wire  pitch  to  diameter  ratio  between  5  and  7  for  steam 
condensation.  They  achieved  an  enhancement  between  1.5  and 
1.8  for  the  heat  transfer  coefficient  over  plain  tubes. 

For  low  surface  tension  fluids  (R-ll  and  ethanol), 
Fujii  et  al .  [Ref.  30]  obtained  results  2  to  3  times  those 
predicted  by  Nusselt  theory,  with  an  optimum  wire  pitch  to 
diameter  ratio  of  2.  O'Keefe  [Ref.  18]  compared  the 
enhancement  ratio  versus  fractional  coverage  (the  percentage 
of  the  tube's  surface  covered  by  wire)  on  a  titanium  smooth 
tube.  Together  with  data  obtained  by  Mitrou  [Ref.  14]  and 
Sethumadhavan  and  Rao  [Ref.  34],  he  determined  that  the 
optimum  fractional  coverage  was  between  10  and  30%.  This 
fairly  broad  range  of  values  was  attributed  to  the  possible 
effects  of  the  different  tube  material  thermal  conductivities 
used  in  the  studies. 
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3 .   Corrugated  Tubes 

In  1971,  Withers  and  Young  [Ref.  35]  performed  a 
design  study  of  two  condenser  bundles  with  heat  duty  of  25 
million  Btu/hr,  one  with  bare  (smooth)  tubes  and  the  other 
with  corrugated  (roped)  tubes.  Coolant  flow  pressure  drop  was 
one  of  the  matching  parameters  for  the  comparison.  This 
requirement  made  it  necessary  to  use  a  larger  corrugated  tube 
(1.00  inch  outside  diameter  (OD) ) ,  than  the  smooth  tube  (0.75 
inch  OD) .  They  found  an  enhancement  of  88%  for  the  outside 
bundle  heat  transfer  coefficient.  Catchpole  and  Drew  [Ref. 
36]  tested  single  corrugated  tubes  and  determined  that  the 
steam  side  heat  transfer  coefficient  increased  with  increasing 
groove  density  and  passed  through  a  maximum  with  increasing 
groove  depth.  Their  conclusion  was  that  the  smaller  the  helix 
angle  (measured  from  the  tube  axis) ,  the  worse  the  condensate 
drainage,  resulting  in  a  reduction  of  the  outside  heat 
transfer  coefficient.  Similar  results  were  reported  by  Marto 
et  al .  [Ref.  3]  and  Mehta  and  Rao  [Ref.  37] .  Eleven  different 
corrugated  tubes  of  three  primary  types  were  tested  by  Marto 
et  al .  [Ref.  3] .  Within  each  corrugation  type,  the  largest 
value  for  the  outside  heat  transfer  coefficient  was  obtained 
at  small  pitch  to  diameter  ratios  (P/D) .  As  P/D  increased, 
the  enhancement  ratio  decreased.  In  tests  conducted  by  Mehta 
and  Rao  [Ref.  37],  they  found  the  condensing  coefficient 
enhancement  ratio  decreased  from  1.26  to  1.1  when  the 
corrugation  groove  pitch  was  increased  from  6.35  mm  to  12.5  mm 
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with  a  constant  groove  depth  of  0.63  mm.  They  also  tested  the 
effect  of  holding  the  groove  pitch  to  groove  depth  ratio 
constant.  A  maximum  enhancement  ratio  of  1.38  was  achieved 
with  the  smallest  groove  depth  and  pitch.  Zimparov  et  al . 
[Ref.  38]  confirmed  the  same  trend,  using  27  mm  outside 
diameter  corrugated  tubes,  with  the  enhancement  ratio 
increasing  from  0.99  to  1.22  for  a  corrugation  groove  pitch 
decreasing  from  16.9  mm  to  6.5  mm. 

Cunningham  and  Bahjernejad  [Ref.  39]  examined  the 
effects  of  condensate  inundation  on  corrugated  tube  bundles. 
They  tested  two  different  multistart  tubes,  one  with  two 
helical  starts  and  the  other  with  six  starts.  Their  results 
indicated  severe  limitations  in  performance  due  to  inundation 
effects  for  tube  bundles  more  than  a  few  rows  deep. 
Additional  inundation  tests  were  conducted  on  multistart  tubes 
by  Cunningham  and  Boudinar  [Ref.  40] .  They  reported  similar 
results  as  those  of  Cunningham  and  Bahjernejad  [Ref.  39]  ,  and 
additionally  found  that  the  more  severe  the  indentation 
(groove  depth)  on  the  roped  tube,  the  greater  was  the 
reduction  in  vapor  side  heat  transfer  for  the  tube  bundle. 

Rabas  [Ref.  41]  studied  an  extensive  data  bank  of 
results  published  by  numerous  researchers.  He  found  that  the 
enhancement  in  the  outside  heat  transfer  coefficient  was 
negligible  for  tubes  with  grooved  geometries  that  are  suitable 
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for  ball  cleaning1  of  the  inside  surface.  He  also  compared 
the  collected  data  with  the  results  predicted  by  three 
different  correlations.   The  Wolverine  correlation  [R-f .  2]  : 


ho     .  0.725  Cn  kS»?*hf 

n\nt  ixfD0    (TsaC-Two) 


-11/4 


(2.7) 


where 


Cn  =  Z*    [nc 


PW 


(2.8) 


Z*  =  38  (ot) 


0  .65 
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(2.9) 


and: 


PW  =  [(5.31E-7)     +  [0.195  (o>)  °-»2'»»f  ]1/6 


(2.10) 


•  Values   for  oj  and   0*   are   listed   for   specific   tube 
geometries  by  Wolverine  in  Reference  2 . 

•  The  variable  Cn  in  Equation  (2.7)  must  be  applied  using 
English  units. 


and  the  correlation  developed  by  Mehta  and  Rao  [Ref .  37] : 


(V  P/  glv-f 


=  1.54  e205e 


4r 


-0.33e"7be 


(2.11) 


for  0  <  6    <    15xl0"4,  and 


1ball  cleaning  is  a  method  of  removing  deposits  from  the 
inside  of  condenser  tubes  by  periodically  recirculating  sponge 
rubber  balls  through  the  tubes  with  the  coolant. 
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(k/  pf2  g/\Lf2) 
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ilV0,3 


0.33e' 
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for  15xl0~4  <  0  <  200xl0~4,  were  both  found  to  overpredict  the 
outside  heat  transfer  coefficient.  Values  similar  to  those 
obtained  from  the  experimental  data  were  predicted  by  the 
Yorkshire  correlation  [Ref.  42]  given  by: 


h. 


cc 


kf2   p 


n,r 


1/3 


(2.13) 


where : 


cs   =   1  -  #c 


fL°l 

p 

■0  .75 


(2.14) 


and  $c  varies  with  tube  profile  and  can  be  found  in  Reference 
42.  However,  Rabas  [Ref.  41]  cautioned  its  use  due  to  the 
results  being  close  to  unity.  In  conclusion,  he  recommended 
using  smooth  tube  correlations  and  neglecting  the  minimal 
enhancement  associated  with  corrugated  tubes  for  design 
calculations . 

Dreitser  et  al .  [Ref.  43]  determined  a  correlation  for 
the  enhancement  ratio  of  corrugated  tubes  based  on  three 
primary  geometric  ratios:  1)  The  ratio  of  the  radius  of 
curvature  of  the  corrugation  ridge  to  the  outside  diameter, 
R/DQ;  2)  the  ratio  of  groove  pitch  to  the  outside  diameter, 
P/DQ;  and  3)  the  ratio  of  the  external  groove  diameter  to  the 
outside  diameter,  d0/D0. 


15 


h. 
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=  2 .469 
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1-0  .  379 —  |exp 


3.651 


(2.15) 


They  found  that  if  any  of  the  three  ratios  were  increased 
while  keeping  the  others  constant,  the  enhancement  ratio  would 
decrease . 

Cunningham  and  Fidler  [Ref.  44]  made  a  comparative 
study  of  the  effects  of  vapor  shear  on  multistart  corrugated 
tubes  and  smooth  tubes.  They  found  that  the  corrugated  tubes 
maintained  their  enhancement  advantage  over  the  smooth  tube 
when  subjected  to  increased  vapor  velocity. 
4.   Wire-Wrapped  Corrugated  Tubes 

Wire -wrapping  of  corrugated  tubes  has  been  shown  by 
various  researchers  to  significantly  reduce  the  effect  of 
condensate  inundation  (Memory  et  al .  [Ref.  45],  Union  Carbide 
[Ref.  10] ,  and  Marto  and  Wanniarachchi  [Ref.  46] ) . 

Union  Carbide  [Ref.  10]  tested  steam  condensation  for 
a  wire-wrapped  roped  tube  in  a  bundle  arrangement.  They 
reported  nearly  identical  performance  for  both  the  unloaded 
and  the  condensate  loaded  conditions  tested.  Condensing 
coefficients  approximately  2.5  times  that  of  plain  tube 
results  were  achieved  for  the  condensate  loaded  condition. 
Memory  et  al .  [Ref.  45]  obtained  a  60  to  68%  enhancement  over 
smooth  tube  results  when  condensing  R-113  on  a  wire-wrapped 
roped  tube  bundle  arrangement  containing  four  tubes  in  a 
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vertical  column.   Their  tests  indicated  an  optimal  wire  pitch 
to  diameter  ratio  of  7.5. 

O'Keefe  [Ref .  18]  tested  the  enhancement  effects  of 
three  different  wire  diameters  (1.6  mm,  1.0  mm,  and  0.5  mm) 
when  used  to  wrap  an  LPD  Korodense  tube.  Enhancements  ranging 
between  20  and  30%  above  smooth  tube  results  were  obtained 
with  the  maximum  obtained  for  the  0.5  mm  wire  with  a  wire 
pitch  to  wire  diameter  ratio  of  14. 

C.   COOLANT  SIDE  CONSIDERATIONS 
1.   Smooth  Tubes 

Numerous  correlations  have  been  developed  for  use  in 
the  analysis  of  the  coolant  side  heat  transfer  correlation 
associated  with  the  turbulent  flow  regime  inside  a  pipe  (Re  > 
10,000)  .  Many  of  the  earlier  correlations  have  been  presented 
in  the  form: 

Nu  =  C1RemPrn  (2.16) 

This  form  has  been  used  for  several  well-known 
correlations  including  those  developed  by  Dittus  and  Boelter 
[Ref.  47] : 

Nu   =  0.023  Re°-8Pr°-i  (2.17) 

and  Colburn  [Ref.  48] 
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Nu  =   0.023  Re°-aPrl/2  (2.18) 

A  correction  factor  for  equation  (2.18)  was  developed 
by  Sieder  and  Tate  [Ref .  49] : 

Nu  =   0.027  Re°-aPr1/2(\ic/\iw)°-li  (2.19) 

This  modification  compensates  for  the  variation  in  the  coolant 
viscosity  when  large  temperature  differences  exist  between  the 
bulk  coolant  and  the  inner  tube  wall  temperature.  All  fluid 
properties  for  equations  (2.17),  (2.18),  and  (2.19)  are 
evaluated  at  the  mean  coolant  bulk  temperature,  Tm. 

(T    +  T  ) 
T    =   — 1 111  (2.20) 

2 

In  more  recent  years,  Sleicher  and  Rouse  [Ref.  50]  and 
Petukhov  [Ref.   51]   developed  equations   (2.21)   and   (2.24) 
respectively,  which  are  applicable  over  a  wider  range  of 
Prandtl  numbers. 


where 


Nu  =  5   +  0.015  Refc  Prwd  (2.21) 


c  =  0.88   -       °'24  (2.22) 

4  +  Prw 


d  =   —   +0.5  exp(-0.6   Prw)  (2.23) 


and, 


NU  =   (£/8)    Re  Pr     ■ (2.24) 

Kr  +  K2(c/S)1/2    (Pr2/3-l) 
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where : 

e  =  (1.82  log  (fie)  -  1.64)~';-  (2.25) 

/Cj.  =  1  +  3  .  4  e  (2.26) 

K2   =  11.7  +  1.8  Pr'1/2  (2.27) 

Lorenz  et  al  .  [Ref .  52]  compared  their  experimental  results 
for  cold  water  (i.e.  large  Prandtl  number)  with  seven  of  the 
more  common  inside  heat  transfer  correlations.  They  found 
that  the  Petukhov- Popov  and  Sleicher- Rouse  correlations  agreed 
to  within  5%  of  the  experimental  data.  The  other  correlations 
typically  underpredicted  the  data  by  between  5  and  15%. 
2 .   Corrugated  Tubes 

The  internal  ridges  of  the  tube  wall  produced  when 
manufacturing  corrugated  tubes  cause  two  primary  effects.  The 
irregularity  of  the  wall  surface  results  in  an  increased 
frictional  effect  and  causes  the  flow  to  swirl,  thereby 
requiring  additional  pumping  power  to  achieve  flow  rates  equal 
to  similar  diameter  smooth  tube  flow.  The  increase  in 
turbulence  also  has  the  effect  of  increasing  the  inside  heat 
transfer  rate.  In  the  study  by  Withers  and  Young  [Ref.  35], 
the  total  inside  enhancement  ratio  of  the  tube  bundle  tested 
was  56%  above  that  obtained  for  a  smooth  tube  bundle  designed 
with  an  equal  pressure  drop.  Catchpole  and  Drew  [Ref.  36] 
found  the  inside  heat  transfer  coefficient  to  increase  with 
increasing  groove  depth  and  groove  frequency.   However,  they 
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stressed  that  the  increase  in  performance  was  obtained  at  the 
expense  of  significantly  greater  friction  factors  and  therefor 
increased  pumping  power  and  operating  cost.  In  comparing 
their  results  with  the  Prandtl  analogy  [Ref.  53]  they  found 
that,  due  to  the  turbulence  caused  by  the  presence  of  grooves, 
the  inside  heat  transfer  coefficient  followed  the  trend 
predicted  by  the  Prandtl  analogy. 

Mehta  and  Rao  [Ref.  37]  presented  their  findings  in 
terms  of  the  severity  factor,  the  ratio  of  groove  depth 
squared  to  the  product  of  groove  pitch  and  inside  diameter. 
This  dimensionless  group  was  chosen  because  it  has  been  found 
to  be  a  useful  criterion  for  evaluation  of  heat  transfer  and 
frictional  characteristics  of  plain  and  spirally  enhanced 
tubes.  An  improvement  of  260%  was  obtained  above  smooth  tube 
results  with  the  higher  inside  heat  transfer  coefficient 
corresponding  to  higher  severity  factors.  Of  the  eleven  tubes 
tested  by  Marto  et  al .  [Ref.  3]  ,  inside  enhancement  was 
reported  for  all  the  tube  geometries  examined,  with  the 
magnitude  ranging  from  1.6  to  5.1  times  a  smooth  tube. 
Zimparov  et  al .  [Ref.  38]  tested  nine  tubes  with  differing 
corrugation  geometries  and  obtained  results  similar  to 
previous  researchers.  They  achieved  inside  enhancement  ratios 
between  1.77  and  2.73,  with  an  increase  in  the  friction  factor 
of  2  to  4  fold. 

For  corrugated  tubes,  Withers  [Ref.  54]  modified 
existing  correlations  and  developed  the  following: 
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Nu   = 


/TT8  Re  Pr 

7  .22(P/Diyl/2  y/Pr  {(e/DJRe/TTsf-^+y 


(2.28) 


where 


Y  = 


2.5  In 


2  — 


(2.29) 


and 


\ 


f 


2.46  ln[r  +  (1  /Re) 


(2.30) 


The  variables  m  and  r  were  determined  for  each  tube  by  a 
curve- fitting  procedure.  He  obtained  agreement  between 
experimental  and  predicted  data  to  within  10%.  Wolverine 
currently  utilizes  a  form  of  this  correlation  in  their 
Engineering  Data  Book  II  [Ref.  2]. 

In  a  comparison  of  the  LPD  and  MHT  Korodense  tubes 
manufactured  by  Wolverine  [Ref.  2],  they  report  a  25%  higher 
inside  heat  transfer  coefficient  for  the  MHT  Korodense  tube 
than  for  the  LPD  Korodense  tube.  The  MHT  Korodense  tube  also 
has  twice  the  pressure  drop  of  the  LPD  Korodense  tube. 

After  compiling  an  extensive  collection  of  published 
data,  Rabas  et  al .  [Ref.  55]  recommended  a  minor  modification 
to  the  leading  coefficient  in  the  Yorkshire  correlation  [Ref. 
42]  to  better  fit  the  experimental  data.  The  Yorkshire 
correlation  is  similar  in  form  to  equation  (2.17)  with  an 
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additional  complex  leading  coefficient  (E±)    to  compensate  for 
the  corrugation  tube  geometry. 


where 


Nu   =  0.027  El  Re°-S  Pr0-*  (2.31) 


^  =  1.0  +  1.182  £2  Ej  (P/Z^)0"7  (2.32) 


E2  =   cos (a)  (2.33) 


a  =  tarrMP/nDj)  (2.34) 

ln(E3)    =  -  t   aj[ln(e/D1)]J  (2.35) 

J=0 


and 


a0   =   5.0766,      a1    =  4.4224,      a2    =   1.1240,      <33    =  6.700x10 


-2 


The  modification  by  Rabas  et  al .   [Ref .   55]   changed  the 
exponent  in  equation  (2.32)  from  0.7  to  0.406. 

Panchal  and  France  [Ref.  56]  conducted  a  study  of 
large  pitch  spirally  indented  tubes  for  possible  application 
to  pressure  drop  sensitive  systems.  Their  data  resulted  in 
internal  enhancements  as  high  as  100%  over  smooth  tube  data 
and  they  reported  agreement  within  a  20%  error  band  for  the 
results  predicted  by  the  correlation  proposed  by  Mehta  and  Rao 
[Ref.  57] ,  given  by: 
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Nu  =   JT/2  Re  Pr (2>36) 

P  Pr°-55(^e/772)0-11  +  y 


where : 


P  =  9.4(P/Di)-°-15(e/Z?i)0-11  (2.37) 

and  7  is  defined  as  in  equation  (2.29)  .  It  should  be  noted 
that  (3  in  equation  (2.37)  is  constant  for  any  given  tube 
configuration . 

3.   Tubes  With  Inserts 

Numerous  types  of  inserts  are  available  for  use  in 
enhancing  the  inside  heat  transfer  coefficient.  Some  of  the 
types  used  by  previous  researchers  at  NPS  include:  twisted 
tape,  wire -wrap,  and  HEATEX  (a  wire  mesh  insert  supplied  by 
CAL  GAVIN) .  A  number  of  correlations  have  been  developed  for 
determining  the  inside  heat  transfer  coefficient  when  using 
these  various  inserts  [Refs.  58-60]  .  These  correlations  for 
HEATEX,  however,  are  for  specific  loop  densities  of  the  wire 
matrix  and  are  not  universally  valid.  HEATEX  inserts  not  only 
significantly  increase  the  overall  heat  transfer  rates,  but 
also  inhibit  fouling.  Gough  and  Rogers  [Ref.  61]  presented 
results  collected  from  a  condenser  fitted  with  HEATEX  inserts 
which  was  in  operation  for  over  four  years  without  being  shut 
down  for  cleaning.  Their  findings  indicated  the  actual 
fouling  resistance  was  over  96%  less  than  the  allowable  limits 
after  a  service  life  of  26  times  that  of  the  original  unit 
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with  no  insert.  However,  it  should  be  noted  that  the  use  of 
an  insert  has  a  pronounced  effect  on  increasing  the  pressure 
drop  across  the  tube.  Inserts  are  used  in  experimental 
applications  to  realize  a  larger  temperature  drop  in  the 
coolant  flow,  thus  decreasing  the  uncertainty  of  the  heat 
transfer  calculations. 
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III.   EXPERIMENTAL  APPARATUS 

A.   SYSTEM  OVERVIEW 

The  apparatus  used  for  this  research  was  originally 
constructed  by  Krohn  [Ref .  12] .  The  initial  design  has  been 
modified  by  several  previous  researchers  (in  particular 
Swensen  [Ref.  17])  and  has  been  used  successfully  by  Coumes 
[Ref.  15],  Guttendorf  [Ref.  16],  Swensen  [Ref.  17]  and  O'Keefe 
[Ref.  18]  to  test  condensation  on  single  horizontal  tubes  of 
various  configurations.  Figure  2  contains  a  general  schematic 
of  the  overall  system.  The  boiler,  a  0.30  meter  diameter 
Pyrex  glass  cylinder,  contains  ten  4  kW,  440  Volt  Watlow 
immersion  heaters  which  are  used  to  generate  steam  from 
distilled  water.  From  the  boiler,  steam  passes  up  through  a 
cylindrical  section  of  Pyrex  glass  with  an  inside  diameter  of 
0.15  meters  and  a  length  of  2.13  meters.  Two  90  degree  Pyrex 
glass  elbows  redirect  the  steam  back  down  a  second  similar 
cylindrical  section  of  1.52  meters  in  length.  The  steam  then 
enters  the  stainless  steel  test  section  containing  the 
horizontal  condenser  tube  mounted  as  shown  in  Figure  3.  The 
test  section  also  contains  a  circular  glass  viewing  port  so 
the  condensation  process  can  be  observed  while  conducting 
tests.  Excess  steam,  not  condensed  on  the  sample  tube,  passes 
through  the  test  section  and  into  the  final  Pyrex  glass 
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Figure  2.   Schematic  of  the  Single  Tube  Test  Apparatus 
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Figure  3.   Schematic  of  the  Test  Section 
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cylinder  containing  an  auxiliary  condenser.  The  auxiliary 
condenser  section  collects  all  the  condensate,  where  it 
returns  to  the  boiler  via  a  gravity  drain  in  the  baseplate. 

Cooling  water  for  the  single  horizontal  tube  in  the  test 
section  was  supplied  by  two  centrifugal  pumps  connected  in 
series  and  drawing  suction  from  an  adjacent  cooling  water 
sump.  A  throttling  valve  and  calibrated  flowmeter  allowed  for 
control  of  the  cooling  water  flow  rate.  The  auxiliary 
condenser  cooling  water  was  supplied  directly  from  the 
building  water  main.  After  setting  the  desired  coolant  flow 
rate  through  the  single  horizontal  tube,  the  saturation 
temperature  in  the  test  section  was  controlled  by  adjusting 
the  throttle  valve  on  the  auxiliary  condenser  coolant  flow 
line . 

A  suction  port  at  the  base  of  the  auxiliary  condenser 
section  allowed  for  the  removal  of  non- condensible  gases. 
Figure  4  shows  a  schematic  of  the  associated  system.  Non- 
condensibles  drawn  from  the  auxiliary  condenser  passed  through 
an  internal  condensing  coil  where,  any  associated  steam  that 
might  be  present  was  condensed  and  collected  in  a  plexiglass 
vacuum  chamber  for  later  removal.  The  non-condensibles  then 
passed  through  the  vacuum  pump  and  were  expelled  to  the 
atmosphere . 


28 


Q-  <-* 


Figure  4.   Schematic  of  the  Purging  System  and  Cooling  Water 
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B.   SYSTEM  INSTRUMENTATION 

The  440  Volt  immersion  heater  power  supply  was  controlled 
by  a  silicon  controlled  rectifier  mounted  on  the  electrical 
switchboard.  Poole  [Ref.  62]  details  the  power  calculation 
for  the  input  into  the  data  acquisition  system. 

Three  separate  methods  were  used  to  monitor  the  pressure 
of  the  system: 

1.  Setra  model  204  pressure  transducer. 

2.  System  saturation  temperature  converted  into  pressure. 

3.  Heise  solid  front  pressure  gage  (visual  reference  only) . 

Temperature  measurements  were  obtained  using  several 
devices.  The  vapor  temperature  of  the  system  was  measured  by 
a  Teflon  coated  and  a  metal  sheathed  type-T  copper/constantan 
thermocouple  located  just  upstream  of  the  horizontal  tube 
being  tested.  Condensate  return  and  ambient  temperatures  were 
also  measured  with  Teflon  coated  type-T  copper/constantan 
thermocouples.  The  inlet  and  outlet  coolant  temperature  of 
the  condenser  test  tube  were  measured  with  three  separate 
probes : 

1.  A  Teflon  coated  type-T  copper/constantan  thermocouple. 

2.  A  ten  junction,  Teflon  coated  type-T  copper/constantan 
thermopile . 

3 .  An  HP  2804A  quartz  crystal  thermometer. 
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The  three  temperature  probes  at  the  outlet  of  the  horizontal 
tube  were  located  downstream  of  a  mixing  chamber  in  order  to 
obtain  a  representative  average  coolant  discharge  temperature. 
Figure  3  contains  the  locations  of  these  temperature  probes. 

An  HP-3497A  data  acquisition  system,  linked  to  a  HP-9826A 
computer,  was  used  to  collect  the  voltage  readings  from  all 
the  measurement  devices.  The  data  collection  and  reduction 
program  "DRPML",  listed  in  Appendix  C,  allowed  for  the 
processing  and  storage  of  all  the  data  by  the  HP-9826A 
computer. 

Swensen  [Ref.  17]  provides  calibration  data  for  the 
pressure  transducer,  thermocouples,  quartz  thermometer  and 
flowmeter. 

C.   TUBES  TESTED 

Eight  enhanced  tubes  were  manufactured  for  this  thesis. 
Tests  were  also  conducted  on  tubes  studied  by  previous 
researchers  to  ensure  repeatability  of  results.  Table  I 
contains  a  listing  of  all  the  tubes  tested  and  additional 
tubes  used  in  the  comparison  of  results.  All  of  the  tubes 
assembled  for  this  study  were  made  from  a  single  titanium  MHT 
Korodense  tube  manufactured  by  Wolverine.  Seven  of  the  tubes 
were  then  wrapped  with  wire  having  diameters  ranging  from  0.1 
to  1.6  mm  to  test  enhancement  effects.  To  attach  the  wire  to 
the  tubes,  two  different  methods  were  used.  The  larger 
diameter  (stiffer)  wires  (a  0.5  mm)  were  wrapped  around  a  rod 
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of  a  smaller  diameter  than  the  tube,  this  formed  the  wire  into 
a  tight  coil.  The  wire  was  then  fed  onto  the  tube  by  twisting 
the  coil  and  advancing  it  along  the  groove  of  the  tube.  The 
wire  was  then  spot  welded  in  place  at  each  end  of  the  active 
condensing  surface.  For  the  wires  <  0.5  mm,  one  end  was  spot 
welded  to  the  tube  at  the  edge  of  the  active  surface.  The 
wire  was  then  tightly  wrapped  around  the  tube,  by  hand,  in  the 
groove  and  spot  welded  to  the  other  end  of  the  active  surface. 
Both  methods  resulted  in  a  constant  wire  pitch  equal  to  the 
corrugation  pitch  of  the  tube.  The  wire  was  not  physically 
attached  to  the  active  tube  surface  and  did  not  act  as  a  fin. 
The  additional  tubes  included  in  Table  I  were  manufactured  for 
research  conducted  by  O'Keefe  [Ref .  18]  and  Mitrou  [Ref .  14] . 
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IV.   EXPERIMENTAL  PROCEDURES  AND  DATA  ANALYSIS 

A.   TUBE  PREPARATION 

To  ensure  consistent  results  in  data  collection,  a 
continuous  film  of  condensate  must  be  present  over  the  entire 
tube  surface.  A  chemical  treatment  procedure  used  by 
Guttendorf  [Ref.  16]  and  several  other  researchers  at  NPS  was 
used  to  produce  filmwise  condensation  on  the  tube  surface. 
Prior  to  installation  into  the  test  section,  the  tube  was 
chemically  treated  as  follows: 

1.  ^he  inside  and  outside  of  the  tube  were  thoroughly 
cleaned  with  mild  soap  using  a  soft  bristle  brush.  The 
soap  was  rinsed  off  with  distilled  water,  then  the  tube 
was  rinsed  with  acetone.  Distilled  water  was  again 
poured  over  the  tube,  to  ensure  the  water  formed  a 
continuous  film  on  the  tube  surface.  If  breaks  in  the 
film  were  observed  (an  indication  of  impurities  on  the 
tube  surface)  the  cleaning  process  was  repeated.  Once  a 
continuous  film  was  present,  the  active  surface  was  not 
handled,  thus  minimizing  the  possibility  of 
contamination. 

2.  The  tube  was  then  placed  in  a  steam  bath. 

3 .  A  solution  of  50%  by  weight  sodium  hydroxide  and  ethyl 
alcohol  was  mixed  for  treating  the  tube  surface.  The 
solution  was  kept  warm  to  ensure  a  watery  consistency. 

4.  The  solution  was  then  applied  over  the  tube  surface,  with 
a  small  brush,  in  10  minute  intervals  for  one  hour. 
While  applying  the  solution,  the  tube  was  rotated  so  the 
entire  surface  could  be  treated.  If  the  tube  was  treated 
previously,  the  solution  was  only  applied  every  5  minutes 
over  a  20  minute  period. 

5.  After  the  tube  surface  was  treated,  it  was  removed  from 
the  steam  bath  and  rinsed  with  distilled  water  to  remove 
any  excess  alcohol/sodium  hydroxide  solution.  To  verify 
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a  continuous  film  will  be  produced  in  the  presence  of 
steam,  the  tube  was  returned  to  the  steam  bath,  to  ensure 
the  entire  surface  wets  easily. 

Without  touching  the  active  surface,  the  tube  was  removed 
from  the  steam  bath  and  immediately  mounted  in  the  test 
section  of  the  boiler  apparatus.  If  tests  were  to  be 
conducted  with  an  insert,  the  insert  was  installed  prior 
to  securing  the  coolant  flow  piping  to  the  test  section. 


The  oxide  layer  formed  by  the  chemical  treatment  promotes 
a  continuous  film  on  the  tube  surface.  Due  to  the  negligible 
thickness  of  the  oxide  film,  its  thermal  resistance  was 
considered  insignificant  compared  to  the  overall  thermal 
resistance  of  the  tube. 

B.   SYSTEM  START-UP  AND  SHUTDOWN  PROCEDURES 

After  the  tube  has  been  installed  in  the  test  section,  the 
system  can  be  brought  on-line  using  the  following  procedures: 

1.  The  boiler  section  must  be  filled  with  distilled  water  to 
a  level  approximately  4  to  6  inches  above  the  top  of  the 
heater  elements.  The  boiler  can  be  filled  by  attaching 
a  hose  from  the  distilled  water  tank  to  the  piping  at  the 
boiler  drain/fill  valve.  Opening  the  vent  valve  on  the 
auxiliary  condenser  allows  for  gravity  filling  of  the 
boiler.  To  drain  the  boiler,  disconnect  the  distilled 
water  tank  hose  and  attach  the  hose  leading  to  the  drain 
in  bilge  area  below  the  boiler.  Then  open  the  drain/fill 
valve  and  allow  the  water  to  gravity  drain  from  the 
boiler. 

2.  When  the  appropriate  water  level  has  been  achieved,  shut 
both  the  fill  /drain  valve  and  the  vent  valve  on  the 
auxiliary  condenser. 

3.  Energize  the  data  acquisition  system,  computer,  printer, 
and  quartz  thermometer  power  supply.  Load  the  software 
program  DRPML  and  check  for  proper  operation.  Then 
verify  that  the  thermocouple  outputs  all  correspond  to 
the  same  ambient  temperature. 
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4.  Open  the  fill  valve  to  the  cooling  water  sump  and  adjust 
the  flow  rate  to  ensure  the  drain  box  does  not  overflow. 
The  fill  valve  for  the  sump  is  located  to  the  left  of  the 
boiler  switchboard. 

5.  Perform  a  leak  test  on  the  cooling  water  systems. 

a.  Energize  the  cooling  water  supply  pumps  and  adjust 
the  rotameter  flow  rate  from  20%  to  80%  while 
checking  for  leaks  in  the  test  section.  Then  secure 
the  control  valve  and  de-energize  the  cooling  water 
supply  pumps . 

b.  Open  the  control  valve  for  the  auxiliary  condenser 
and  adjust  the  flow  rate  to  at  least  30%  and  check 
for  leaks  in  the  system.  After  completing  the  leak 
check,  reset  the  flow  rate  to  the  auxiliary 
condenser  to  the  desired  position,  minimum  of  10%. 

6.  Place  the  system  under  a  vacuum.  To  do  this,  shut  the 
drain  valve  in  the  plexiglass  cold  trap  and  energize  the 
vacuum  pump.  Check  the  pressure  gage  on  the  vacuum  pump 
to  ensure  proper  operation  and  open  the  suction  valve 
located  on  the  side  of  the  auxiliary  condenser.  After 
the  vacuum  pump  has  reduced  the  system  pressure  to  below 
3  psia,  shut  the  suction  valve  and  secure  the  vacuum 
pump. 

7.  With  the  system  under  a  vacuum,  the  heaters  may  be 
energized  .  Three  switches  must  be  placed  in  the  "ON" 
position  to  energize  the  heaters. 

a.  Switch  #3  on  panel  p5 ,  in  the  main  hallway  adjacent 
to  room  H- 106 . 

b.  The  heater  load  bank  circuit  breaker  on  the  left 
side  of  the  boiler  heater  switchboard. 

c.  The  condensing  rig  boiler  power  switch,  on  the  front 
of  the  boiler  heater  switchboard. 

When  the  heaters  are  energized,  the  power  level  should  be 
adjusted  to  50  Volts,  using  the  silicon  control  rectifier 
on  the  boiler  control  switchboard.  If  the  system  is  below 
2  psia,  the  power  level  should  initially  set  to  40  Volts 
to  minimize  vibrational  shock  to  the  system  caused  from 
vapor  bubble  formation  and  collapse.  As  the  system  warms 
up,  the  power  can  be  increased,  in  increments  of  10 
Volts,  until  the  desired  setting  is  reached. 
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8.  After  the  system  pressure  rises  above  4  psia,  non- 
condensible  gases  need  to  be  purged  from  the  system. 
This  can  be  accomplished  by  following  the  procedures 
outlined  in  step  6.  Throughout  the  purging  procedure, 
the  horizontal  tube  coolant  flow  should  remain  secured, 
while  the  flow  rate  to  the  auxiliary  condenser  can  be 
adjusted  until  all  the  gases  have  been  evacuated.  When 
the  auxiliary  condenser  is  warm  to  the  touch  everywhere, 
and  condensate  has  formed  in  the  plexiglass  container, 
most  of  the  non- condensible  gases  have  been  removed.  The 
initial  purging  of  the  system  may  take  between  15  and  30 
minutes.  This  process  should  be  repeated  every  few  hours 
during  extended  operation. 

9.  To  establish  filmwise  condensation  on  the  tube  being 
tested,  the  following  steps  should  be  performed: 

a.  Allow  the  system  vapor  temperature  to  increase  to 
=  90  °C  (at  least  3800  microvolts  on  channel  40) . 

b.  Increase  the  flow  rate  in  the  auxiliary  condenser 
to  50%  or  60%  and  allow  the  vapor  temperature  to 
reduce  to  =75  °C  (approximately  3200  microvolts) . 

c.  Secure  coolant  flow  to  the  auxiliary  condenser  and 
allow  the  vapor  temperature  to  increase  to  about  85 
°C  (=3700  microvolts) .  This  forms  a  steam  blanket 
over  the  tube. 

d.  Initiate  coolant  flow  of  80%  in  the  single 
horizontal  tube. 

e.  Adjust  the  auxiliary  condenser  flow  rate  to  maintain 
the  desired  temperature  and  pressure  for  data 
collection. 

If,  after  performing  steps  9  a.  through  9  e.,  some 
dropwise  condensation  persists,  repeat  step  9  again.  If 
dropwise  condensation  continues,  the  tube  should  be 
removed,  cleaned  and  retreated  in  accordance  with  the 
procedures  listed  in  section  IV.  A.  A  cutout  section  in 
the  insulation  above  the  test  section  allows  for  a  light 
to  be  shown  on  the  upper  surface  of  the  tube  in  order  to 
check  for  a  continuous  film  on  the  tube  surface. 

10.  Run  the  software  program  DRPML  by  pressing  the  "RUN" 
key  on  the  keyboard.  The  program  will  prompt  you 
with  questions  for  the  necessary  information 
required  to  run  as  listed  by  O'Keefe  [Ref.  18]  and 
in  Appendix  C. 
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11.  Ensure  the  system  has  been  operating  at  steady- state 
conditions  for  at  least  30  minutes  prior  to 
continuing  past  the  question  "Enter  the  flowmeter 
reading" . 

12.  Vacuum  runs  are  conducted  at  a  heater  setting  of  90 
Volts  and  1980  ±  10  microvolts  on  channel  40.  This 
corresponds  to  TSat  =48°C,  and  a  vapor  velocity  of  « 
2  m/s. 

13 .  The  atmospheric  runs  are  conducted  at  a  heater 
setting  of  175  Volts  and  4280  ±  10  microvolts  on 
channel  40.  This  corresponds  to  TSat  <=100°C,  and  a 
vapor  velocity  of  =  1  m/s.  Special  care  must  be  used 
when  operating  at  atmospheric  pressure  to  ensure 
over-pressurization  and  rupture  does  not  occur. 

14.  When  taking  readings,  double  check  the  flowmeter 
reading  prior  to  accepting  any  data  point.  If  both 
vacuum  and  atmospheric  pressure  runs  are  to  be 
conducted  in  the  same  day,  the  vacuum  data  should  be 
collected  first.  This  eliminates  the  long  cool  down 
time  required  after  an  atmospheric  pressure  run. 

After  data  has  been  collected,   the  system  should  be 

secured  utilizing  the  following  procedures: 

1.  Secure  power  to  the  heating  elements.  Turn  off  the 
switches  in  the  reverse  order  listed  in  step  7  of  the 
start-up  procedures. 

2.  Secure  coolant  flow  to  the  auxiliary  condenser.  If  the 
system  is  to  remain  at  vacuum  pressure  until  the  next 
data  run,  then  the  auxiliary  condenser  can  be  used  in 
assisting  to  cool  the  system  down,  provided  the  same  tube 
will  be  used  for  the  following  test. 

3.  Secure  the  coolant  flow  through  the  single  horizontal 
tube  by  securing  the  coolant  pumps. 

4.  Secure  the  water  flow  to  the  coolant  sump  tank. 

5.  To  return  the  system  to  atmospheric  pressure,  slowly  open 
the  vent  valve  on  the  auxiliary  condenser.  To  avoid 
contaminating  the  system,  keep  any  foreign  material  away 
from  the  vicinity  of  the  vent  valve  during  this 
procedure . 

6.  If  an  emergency  should  arise,  such  as  an 
over-pressurization  or  breakage,  ensure  the  heater  power 
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is  secured  first.   Then,  let  the  system  cool  down  prior 
to  checking  for  damage. 

C.  DATA  COLLECTION  PROCEDURE 

After  the  tube  has  been  installed  and  the  system  has  been 
operating  for  a  minimum  of  30  minutes,  data  collection  may 
begin.  Each  run  consisted  of  data  points  taken  at  coolant 
flow  rate  settings  from  80%  to  20%  and  back  up  to  80%,  in  10% 
increments.  This  allowed  for  each  data  point  to  be  checked 
for  repeatability  during  the  run.  Before  a  data  point  was 
accepted,  sample  readings  were  taken  to  ensure  the  coolant 
temperature  difference,  the  overall  heat  transferred,  and  the 
saturation  temperature  were  stable.  The  temperature 
difference  between  the  inlet  and  outlet  coolant  flow  must  be 
stabilized  (±0.01  °C)  prior  to  recording  data.  After  setting 
the  coolant  flow  rate,  the  aforementioned  parameters  required 
between  5  and  15  minutes  to  stabilize. 

D.  DATA  REDUCTION  PROCEDURE 

For  a  slice  of  the  condenser  tube,  the  overall  thermal 
resistance  from  vapor  to  coolant  can  be  represented  by  the  sum 
of  the  coolant  side  resistance  (R^  ,  the  fouling  resistance 
(Rf ) ,  the  wall  resistance  (Rw) ,  and  the  vapor  side  resistance 
(RQ)  ,  as  shown  in  equation  (4.1)  . 
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Only  clean  tubes  were  tested  so  the  fouling  resistance  was 
set  to  zero,  therefore: 

Being  convective  in  nature,  the  vapor  and  coolant  side 
resistances  can  be  expressed  as  follows: 


^ 


o 


R>  -  i^  <4-4> 

The  wall  resistance  term  comes  from  the  general  radial 
conduction  equation  for  a  concentric  cylinder: 

RW-^4M  (4.5, 

The  total  thermal  resistance  can  be  expressed  in  terms  of  the 
overall  heat  transfer  coefficient  (UQ)  and  the  effective 
outside  condensing  area  (AQ) ,  as  follows: 

R->  -ik  <4-6) 

The  effective  outside  condensing  area  of  the  tube  can  be 
calculated  using  the  following  equation: 

AQ  =  nD0L  (4.7) 

For  the  inside  effective  area,  two  different  conditions 
are  present.  The  area  corresponding  to  the  outside  effective 
area,  where  heat  transfer  takes  place  radially,  and  the  inlet 
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and  outlet  portions  of  the  tube,  where  heat  removal  occurs 
axially.  To  account  for  this  'fin'  effect,  the  extended  fin 
assumption  and  associated  fin  efficiencies  were  used  in  the 
following  expression: 

Ai   =  KDi(L+L1Tl1+L2Tl2)  (4.8) 

Substituting  equations   (4.3),   (4.4),   and   (4.6)   into 
equation  (4.2)  yields  the  following: 

11  1 

+  R»   +  -r^r-  (4.9) 


From  the  data  collected,  the  total  heat  transfer  rate 
across  the  tube  can  be  calculated  directly  from  an  energy 
balance  using  the  coolant  mass  flow  rate  and  the  coolant 
temperature  difference  (from  inlet  to  outlet)  through  the 
tube : 

Q  =  ihCpi^-T^  (4.10) 

The  total  heat  transfer  rate  can  also  be  expressed  in  terms  of 
the  overall  heat  transfer  coefficient,  the  effective  outside 
condensing  area,  and  the  log  mean  temperature  difference 
(LMTD) : 

Q  =  UaA0(LMTD)  (4.11) 

where: 

(r2-r,) 

LMTD  = 

4.12 


In 


T      -T 

■^  sat     ±  l 

T       -T 
1  sat     12 
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The  measurements  for  LMTD  were  obtained  in  the  following 
manner : 


•  T-l  was  measured  from  a  quartz  thermometer  installed  at  the 
inlet  section  of  the  coolant  flow  for  the  condenser  tube. 

•  T2  was  measured  with  a  quartz  thermometer  inserted 
downstream  of  the  mixing  chamber  at  the  outlet  section  of 
the  coolant  flow  for  the  condenser  tube.  A  correction 
factor  was  applied  to  T2  in  order  to  compensate  for 
viscous  heating  of  the  coolant  through  the  tube.  The 
corrections,  which  depend  on  the  inside  geometry  of  the 
tube,  the  coolant  flow  velocity,  and  the  type  of  insert, 
are  listed  in  Appendix  A. 

•  Tsat  was  measured  using  the  vapor  thermocouple  (channel 
40)  . 


Substituting  equation  (4.10)  into  equation  (4.11)  provides 
determination  of  the  overall  heat  transfer  coefficient 
directly  from  experimental  measurements: 


.  *CP(T2-Tr) 
°  AQ(LMTD) 

With  AQ,  Ai ,  Rw,  and  UQ  all  known  quantities,  the  only 
remaining  unknowns  in  equation  (4.9)  are  hQ  and  hit  the 
outside  and  inside  heat  transfer  coefficients.  The  modified 
Wilson  plot  technique  was  used  to  calculate  these 
coefficients . 

E.   MODIFIED  WILSON  PLOT  TECHNIQUE 

In  1915,  Wilson  [Ref.  63],  developed  a  method  for 
indirectly   determining   the   inside   and   outside   thermal 
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resistance  from  an  overall  resistance.  Wilson's  method 
requires  the  heat  flux  of  the  system  to  remain  constant. 
While  this  application  proves  valid  in  theory,  difficulties 
arise  when  trying  to  maintain  a  constant  heat  flux  while 
varying  the  cooling  water  velocity,  Marto  [Ref .  64] . 

Khartabil  et  al .  [Ref.  65]  provide  a  comprehensive  listing 
of  modifications  to  the  original  technique  developed  by 
Wilson.  For  a  detailed  development  of  the  modified  Wilson 
plot  technique,  the  reader  is  referred  to  Marto  [Ref.  64]  and 
Marto  and  Nunn  [Ref.  66] . 

Using  the  modified  Wilson  plot  method  to  determine  the 
inside  and  outside  heat  transfer  coefficients  requires  one  to 
choose  a  general  form  for  each  correlation.  The  program  DRPML 
enables  the  user  to  select  from  several  options.  The  inside 
heat  transfer  coefficient  correlation  can  be  written  in 
simplified  form  as  follows: 

^   =  C2Q  (4.14) 

where  Q  varies  with  the  particular  correlation  used.  Using 
the  Sieder-Tate  correlation  [Ref.  49]  (equation  (2.19)): 

Q  =  iEs  Rex  Pi^Hv- r/nJ°'14  (4.15) 

Di  c 

where  the  program  allows  the  Reynolds  number  exponent  (x)  to 
be  varied.   Using  the  Sleicher-Rouse  correlation  [Ref.  50] 
(equation  (2 .21)  )  : 
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Q  =  —  (5  +  0.  015  Refc  Pi. 


(4.16) 


where  c  and  d  are  defined  in  equation  (2.22)  and  (2.23 
respectively.  Using  the  Petukhov- Popov  correlation  [Ref.  51 
(equation  (2 .24) )  : 


e/8)  Re  Pi 


Kx   +  K2(e/Q)  l/2    (Pr2/3-l 


(4.17) 


where  e,  K1#  and  K2  are  defined  in  equation  (2.25),  (2.26), 
and  (2.27)  respectively. 

Similarly,    the   outside   heat   transfer   coefficient 
correlation  may  be  written  in  simplified  form: 


hQ  =  a  Z 


(4.18) 


where  Z  can  be  selected  from  equation  (2.1)  (Nusselt's  theory 
[Ref.  23] ) : 


Z  = 


kf2  g  pf2  hfg 
Hf  Do    (Tsat-Two) 


1/4 


(4.19) 


or    from  equation    (2.62)     (Fujii's    correlation    [Ref.    29]) 


Z  =    h    FUs   R&2U2 


(4.20) 


F   and   Re2(^   are   defined   in   equation   (2.4)   and   (2.5 
respectively. 
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Substituting  equation   (4.14)   and   (4.18)   into  equation 
(4.9)  yields: 


tfcA 


which  can  be  written  as 


1  1 

+  R,.  + 


C1QA1 


a.ZA, 


(4.21) 


L7 


■R^o 


z  =  — 
c, 


(4.22) 


By    letting 


Y 


-jr-SA 


Z 


(4.23) 


AoZ 

X  =  ^Q 


(4.24) 


m,   = 


_1_ 


(4.25) 


and: 


b  =   ± 


(4.26) 


Equation  (4.22)  can  be  written  in  simple  linear  form  as 


Y  =  m^  +  b 


(4.27) 


Due  to  the  temperature  dependence  of  the  parameters  Q  and 
Z,  the  program  utilizes  an  iterative  procedure  in  the  solution 
for  each  data  point.   A  least-squares  fit  of  equation  (4.27) 
then  determines  Ci  and  ot.      To  calculate  hQ,  the  known  value  of 


45 


a  can  be  used  directly  in  equation  (4.18)  ;  alternatively,  with 
the  known  value  of  Ci;  hi  can  be  calculated  from  equation 
(4.14),  and  applied  to  equation  (4.9)  in  the  following  form: 


a. 


h1A1 


*,A 


i-i 


(4.28) 


F.   ENHANCEMENT  RATIO 

The  enhancement  ratios  in  this  thesis  are  used  to  express 
a  comparison  of  enhanced  tube  results  to  smooth  tube  results, 
or  to  compare  wire-wrapped  corrugated  tube  results  to  plain 
corrugated  tube  results. 

From  Nusselt's  theory,  heat  flux  can  be  written  as: 


Q=    *     ^saC-Two)3/4 


(4.29) 


where 


a  =  a 


V  g  p/  hfg 


Hf  Do 


1/4 


(4.30) 


Heat  flux  can  also  be  represented  by: 


<ar  =  ho  (Tsac-Two) 


(4.31) 


Equating  equations  (4.29)  and  (4.31)  yields 
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h     =  a    (T      -T    )  'l/i  (4.32) 

11 0  a   \  j.  sac      ±  WQj 

Taking  the  ratio  of  the  outside  heat  transfer  coefficient 
from  an  enhanced  tube,  (hoe)  ,  to  a  smooth  tube  (hos)  ,  for  the 
same  (Tsat-Two)  yields: 

eT   =  ^  =  —  =  —  (4.33) 


Similarly,  the  enhancement  ratio  for  a  wire-wrapped  corrugated 
tube  to  a  plain  corrugated  tube  at  the  same  (Tsat-Two)  can  be 
represented  by: 


C  =  OWC       =  WC       _  WC  IA         TAX 

h  a.  a 

11oc  "c  c 


where  the  subscripts  wc  and  c  refer  to  wire -wrapped  corrugated 
and  plain  corrugated  tubes  respectively. 

For  the  purpose  of  this  study  all  enhancement  results  for 
wire-wrapped  corrugated  tubes  are  based  on  eTc  and  wire- 
wrapped  smooth  tube  enhancements  are  based  on  eT. 
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V.   RESULTS  AND  DISCUSSION 

A.   HEAT  TRANSFER  CORRELATIONS 

Lorenz  et  al .  [Ref .  52]  conducted  a  thorough  investigation 
of  several  inside  heat  transfer  correlations  to  determine 
which  correlation  most  accurately  predicted  their  data  for 
turbulent  flow  of  water  in  smooth  tubes.  Their  findings 
showed  that  the  Petukov- Popov  and  Sleicher-Rouse  correlations 
(equations  (2.24)  and  (2.21)  respectively)  were  the  most 
accurate  in  predicting  the  inside  heat  transfer  coefficient 
over  the  range  of  Prandtl  numbers  from  6.0  to  11.6. 
Consequently,  O'Keefe  [Ref.  18]  made  a  comparison  of  his 
results  obtained  using  both  correlations  to  process  data, 
which  was  collected  from  the  same  experimental  apparatus  used 
in  this  study.  His  results  showed  good  agreement  between  the 
outside  coefficients  obtained  for  smooth  copper  and  smooth 
titanium  tubes,  when  the  two  correlations  were  used  with  the 
modified  Wilson  plot  data  reduction  technique.  Based  on  his 
findings,  O'Keefe  [Ref.  18]  processed  most  of  his  data  using 
the  Petukhov- Popov  correlation.  All  data  collected  in  this 
study  were  processed  using  the  Petukov- Popov  correlation  so  a 
direct  comparison  could  be  made  with  the  results  obtained  by 
O'Keefe  [Ref.  181 . 
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The  data  collection  and  reduction  program,  DRPML,  listed 
in  Appendix  C,  allows  for  the  selection  from  two  possible 
outside  heat  transfer  correlations,  the  Nusselt  theory  or  the 
Fujii  correlation  (equations  (2.1)  and  (2.6)  respectively). 
Figure  5  shows  the  reduced  results  for  h0  obtained  from  using 
these  two  equations  for  a  vacuum  pressure  run  on  a  plain  MHT 
Korodense  tube.  It  can  be  seen  that  there  is  good  agreement 
between  the  two  equations.  This  indicates  that  the  leading 
coefficient,  ot ,  obtained  from  the  modified  Wilson  plot  method 
using  the  Nusselt  theory  compensates  for  any  vapor  velocity 
effect  inherent  in  the  Fujii  correlation.  O'Keefe  [Ref.  18] 
processed  all  of  his  data  using  the  Nusselt  correlation;  so 
for  comparison  purposes,  all  further  calculations  of  the 
outside  heat  transfer  coefficient  in  the  present  work  have 
also  been  made  using  the  Nusselt  theory. 

B.   VERIFICATION  AND  REPEATABILITY  OF  RESULTS 

Several  runs  were  made  using  the  same  tubes  tested  by 
previous  researchers  to  verify  system  operation  and  data 
collection  procedures.  The  five  tubes  that  were  tested 
included  two  smooth  titanium  tubes  and  two  LPD  Korodense 
titanium  tubes  tested  by  O'Keefe  [Ref.  18];  of  these,  one  of 
each  was  wrapped  with  a  0.5  mm  diameter  titanium  wire.  The 
fifth  tube  was  a  smooth  copper  tube  wrapped  with  a  0.5  mm 
diameter  titanium  wire  from  Mitrou's  work  [Ref.  14] .  Figures 
6   and   7  are   comparisons   of   the   outside  heat   transfer 
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Vacuum 


hQ  vs.  ATcf  for  a  Plain  MHT  Korodense  Tube  Under 
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Figure  7.   hQ  vs .  ATcf  for  a  Smooth  Titanium  Tube  at 
Atmospheric  Pressure 
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coefficient  for  a  vacuum  and  an  atmospheric  pressure  run 
respectively  for  data  collected  in  this  study  and  the  study 
conducted  by  O'Keefe  [Ref.  18]  for  the  smooth  titanium  tube. 
The  figures  show  good  agreement  for  the  two  sets  of  data.  The 
separation  in  the  data  sets  seen  in  Figure  6  is  due  to  a 
higher  coolant  water  inlet  temperature  when  O'Keefe  conducted 
his  runs.  Tables  II  through  VII  contain  a  listing  of  all  the 
data  collected  by  O'Keefe,  Mitrou  and  from  the  present  study 
for  comparison  purposes  in  this  study.  The  researcher 
initials  are  as  follows:  (0)=0'Keefe  [Ref.  18],  (L)=Long,  and 
(M) =Mitrou  [Ref.  14]  .  The  average  variation  for  a,  the 
leading  coefficient  for  the  outside  heat  transfer  equation, 
was  less  than  3.3%,  with  90%  of  the  results  differing  by  less 
than  2%.  The  Ci  values  obtained  in  this  study  differed  by  as 
much  as  8.5%  compared  with  values  obtained  by  previous 
researchers.  This  larger  difference  is  possibly  due  to  the 
uncertainty  as  to  which  HEATEX  insert  was  used  in  O'Keefe' s 
experiments.  Four  different  inserts  were  available  and  each 
had  a  slightly  different  length. 

C.   ANALYSIS  OF  CORRUGATED  TUBE  RESULTS 

To  establish  a  base  line  for  the  enhancement  ratios  of  the 
wire-wrapped  MHT  Korodense  tubes,  an  initial  series  of  runs 
were  conducted  to  determine  the  heat  transfer  characteristics 
of  the  plain  MHT  Korodense  tube.  The  data  obtained  are 
presented  in  Table  VIII.    Figure  8  shows  a  near  unity 
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Table  II.   SMOOTH  TITANIUM  TUBE  AT  ATMOSPHERIC  PRESSURE 


HEATEX  Insert 

Data  Run/Researcher 

Ci 

a 

FONMAHT1   (0) 

2.37 

0.78 

F0NMAHT2   (0) 

2.38 

0  .77 

F0NMAHT3   (0) 

2.39 

0.76 

F0NMAHT4    (0) 

2  .41 

0.75 

F0NMAHT5    (0) 

2.20 

0  .  77 

F0NMAHT6   (0) 

2.54 

0.76 

F0NMAHT7   (0) 

2.51 

0.79 

LSNTHA1     (L) 

2.90 

0.79 

No  Insert 

F0NMANT1   (0) 

1.21 

0.75 

F0NMANT2   (0) 

1.18 

0.76 

F0NMANT3   (0) 

1.17 

0.77 

F0NMANT4   (0) 

1.18 

0.76 

F0NMANT5   (0) 

1.24 

0.79 

LSNTNA1     (L) 

1.17 

0.82 

enhancement  of  the  outside  heat  transfer  coefficient  for  the 

MHT  Korodense  tube  (having  a  groove  pitch  of  8  mm)  over  smooth 

titanium  tube  results.   This  is  consistent  with  the  findings 

of  Rabas  [Ref.  41]  where,  in  a  comparison  of  an  extensive 

collection  of  data  on  spirally  indented  tubes  he  found: 

that  there  is  no  appreciable  enhancement  for 
condensation  on  roped  tubes  for  groove  pitches  greater 
than  4  mm. 
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Table  III 


SMOOTH  TITANIUM  TUBE  AT  VACUUM  PRESSURE 


HEATEX  Insert 

Data  Run/Researcher 

ci 

Oi 

FONMVHT3   (0) 

2.55 

0.75 

F0NMVHT4    (0) 

2.29 

0  .  79 

F0NMVHT5   (0) 

2.42 

0.76 

LSNTHV1    (L) 

2.87 

0.78 

LSNTHV2     (L) 

2  .57 

0  .76 

LSNTHV3     (L) 

2.43 

0.79 

LSNTHV4     (L) 

2.50 

0  .  77 

No  Insert 

F0NMVNT2   ( 0 ) 

1.09 

0.82 

F0NMVNT3   (0) 

1.08 

0.79 

F0NMVNT4    (0) 

1.07 

0.85 

F0NMVNT5   (0) 

1.11 

0.82 

LSNTNV1    (L) 

1.13 

0.83 

LSNTNV2    (L) 

1.13 

0.84 

LSNTNV4    (L) 

1.11 

0.83 

A  comparison  of  the  C±  values  for  the  no  insert  runs  on 
the  smooth  tube  (Tables  II  and  III)  and  LPD  tube  (Table  V) 
indicate  increases  in  hi  from  73  to  74%  for  vacuum  and 
atmospheric  pressure  runs  respectively.  As  discussed  earlier, 
Wolverine  [Ref .  2]  reported  a  25%  increase  in  the  inside  heat 
transfer  coefficient  for  the  MHT  Korodense  tube  when  compared 
to  the  LPD  Korodense  tube.  Comparison  of  the  C±  values  for  the 
no  insert  runs  on  the  MHT  tube  (Table  VIII)  and  the  LPD  tube 
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Table  IV.   0 . 5  mm  WIRE-WRAPPED  SMOOTH  TITANIUM  TUBE  WITH 
HEATEX  INSERT 


Atmospheric  Pressure 

Data  Run/Researcher 

c, 

a 

F0NMAH6T2 

(0) 

2.39 

0.99 

FONMAH5T3 

(0) 

2.37 

1.01 

LSW1THA1 

(L) 

2.67 

1.04 

Vacuum  Pressure 

FONMVH6T1 

(0) 

2.16 

0.95 

FONMVH6T2 

(0) 

2.21 

0.95 

LSW1THV1 

(L) 

2.59 

0.94 

(Table  V)  indicate  increases  in  hi  comparable  to  those 
reported  by  Wolverine  [Ref.  2],  21%  for  the  atmospheric 
pressure  runs  and  25%  for  the  vacuum  runs.  This  indicates 
that  as  the  groove  depth  increases  from  the  smooth  tube  (zero) 
to  the  LPD  tube  value  («*0.33  mm)  most  of  the  increase  in  the 
inside  heat  transfer  coefficient  due  to  corrugation  has 
already  been  realized.  The  additional  groove  depth  of  the  MHT 
tube  provides  only  a  21  to  25%  further  increase  of  the  inside 
heat  transfer  above  the  LPD  results  (MHT  tube  groove  depth 
«0 . 5mm) ,  while  the  pressure  drop  doubles  [Ref.  2]. 

Figures  9  and  10  show  the  effect  of  using  the  HEATEX 
insert  on  the  overall  heat  transfer  coefficient  (U0)  .  The 
general  trend  for  U0  to  increase  with  increasing  coolant  flow 
velocity  is  due  to  improved  coolant  mixing  and  a  lower  mean 
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Table  V 


PLAIN  LPD  KORODENSE  TUBE 


Atmospheric  Pressure  HEATEX  Insert 

Data  Run/Researcher 

Ci 

a 

F0NMAHLT2    (0) 

2.90 

0.90 

F0NMAHLT3    (0) 

2.89 

0.92 

LLNTHA1     (L) 

3.11 

0.95 

Atmospheric  Pressure  No  Insert 

F0NMANLT2    (0) 

2.06 

0.92 

F0NMANLT3    (0) 

1.99 

0.94 

LLNTNA2      (L) 

2.17 

0.95 

Vacuum  Pressure  HEATEX  Insert 

F0NMVHLT1   (0) 

2.72 

0.94 

F0NMVHLT2   (0) 

2.67 

0.95 

LLNTHV2      (L) 

3.02 

0.97 

Vacuum  Pressure  No  Insert 

F0NMVNLT2    (0) 

1.87 

0.95 

F0NMVNLT3   (0) 

1.86 

0.96 

LLNTNV3      (L) 

1.99 

0.96 

bulk  fluid  temperature.  Increases  range  from  3%  at  a  coolant 
flow  rate  of  3.83  m/s  up  to  9%  at  a  flow  rate  of  1.02  m/s. 
This  increase  is  attributed  to  the  greater  turbulence  of  the 
coolant  flow  due  to  the  presence  of  the  insert  and  can  be  seen 
in  the  higher  Ci  values  for  the  HEATEX  runs  for  every  tube 
tested.  A  comparison  of  O'Keefe's  data  [Ref.  18]  for  smooth 
titanium  and  LPD  Korodense  titanium  tubes,  tested  with  and 
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Table  VI.   0 . 5  mm  WIRE-WRAPPED  LPD  KORODENSE  TUBE  WITH 
HEATEX  INSERT 


Atmospheric  Pressure 

Data  Run/Researcher 

c, 

a 

FONMAHL3T3    (0) 

2.84 

0.97 

LLW0  5THA1    (L) 

3  .03 

0.98 

LLW05THA2    (L) 

3.06 

0.98 

Vacuum  Pressure 

F0NMVHL3T1   (0) 

2.61 

0.99 

LLW0  5THV1     (L) 

2  .  78 

0.98 

LLW0  5THV2     (L) 

2  .78 

0.98 

without  a  HEATEX  insert  indicates  a  20%  increase  in  U0  for  the 
smooth  tube  and  a  7.9%  increase  in  U0  for  the  LPD  tube  when 
tested  at  a  coolant  velocity  of  «2.5  m/s.  For  a  coolant  flow 
rate  of  =2.5  m/s  the  present  study  shows  a  7.2%  increase  in  U0 
for  the  MHT  tube  when  a  HEATEX  insert  is  used.  This  trend 
indicates  that  a  more  severe  corrugation  groove  depth  (zero 
for  the  smooth  tube,  0.3  3  mm  for  the  LPD  tube  and  0 . 5  mm  for 
the  MHT  tube)  results  in  a  smaller  increase  in  U0  with  the 
addition  of  an  insert.  The  grooves  present  in  the  MHT  and  LPD 
tubes  cause  greater  turbulence  in  the  coolant  flow  than  the 
inside  walls  of  the  smooth  tube.  Therefore,  the  addition  of 
an  insert  for  a  corrugated  tube  (which  al ready  has  turbulent 
flow)  has  a  less  pronounced  effect  on  U0  compared  to  the 
effect  realized  for  a  smooth  tube. 
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Table  VII.   0 . 5  mm  WIRE-WRAPPED  SMOOTH  COPPER  TUBE  WITH 
HEATEX  INSERT 


Atmospheric  Pressure 

Data  Run/Researcher 

C. 

0i 

F0NMAH71C1   (0) 

3.07 

1.50 

LSW71CHA1    (L) 

3  .15 

1.28 

S71A305*      (M) 

2.69 

1.50 

S71A314*      (M) 

2.69 

1.50 

Vacuum  Pressure 

F0NMVH71C2    (0) 

2.76 

1.28 

LSW71CHV1    (L) 

2  .94 

1.23 

S71V296*     (M) 

2.50 

1.27 

Runs  conducted  with  wire-wrapped  insert 

D.   ANALYSIS  OF  WIRE-WRAPPED  CORRUGATED  TUBE  RESULTS 

Seven  different  titanium  wire  diameters  ranging  from  0.1 
mm  to  1.6  mm  were  used  to  wrap  separate  MHT  Korodense  tubes. 
Tests  were  conducted  at  atmospheric  and  vacuum  conditions  and 
the  results  were  compared  to  the  plain  MHT  Korodense  tube  data 
to  determine  the  enhancement  effects  attributable  to  wire- 
wrapping.  Wire -wrapped  smooth  tube  data  obtained  by  previous 
researchers  were  also  used  for  comparison  to  establish  if 
similar  trends  exist. 

The  enhancement  due  to  wire -wrapping  can  be  seen  directly 
in  Figure  11  where  the  averaged  U0  (averaged  for  all  data  runs 
at  each  flow  rate)  are  plotted  for  all  the  wire  diameters. 
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Table  VIII 


PLAIN  MHT  KORODENSE  TUBE 


Atmospheric  Pressure  HEATEX  Insert 

Data  Run 

C. 

Of 

LMNTHA7 

3.34 

0.79 

LMNTHA8 

3.31 

0.79 

LMNTHA9 

3.32 

0.80 

Atmospheric  Pressure  No  Insert 

LMNTNA7 

2.49 

0.80 

LMNTNA8 

2  .47 

0.80 

LMNTNA9 

2.56 

0.80 

Vacuum  Pressure  HEATEX  Insert 

LMNTHV7 

3.07 

0.78 

LMNTHV8 

3.01 

0.79 

LMNTHV9 

3  .06 

0.78 

Vacuum  Pressure  No  Insert 

LMNTNV7 

2.36 

0.77 

LMNTNV8 

2.39 

0.76 

LMNTNV9 

2.43 

0.77 

Here  the  0.4  mm  wire-wrap  appears  to  give  the  best  heat 
transfer  with  an  enhancement  of  about  3%.  Similar  trends  for 
the  outside  heat  transfer  coefficient  can  also  be  seen  in 
Figures  12  and  13.  The  plain  MHT  Korodense  tube  results  are 
included  to  show  the  relative  enhancements  for  each  size  of 
wire.  Figure  12  shows  an  increasing  trend  with  decrease  in 
wire  diameter  starting  from  the  1.6  mm  wire-wrap  (Pw/Dw=5)  and 
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Figure  8.   hQ  vs .  ATcf  for  a  Smooth  Titanium  and  a  Plain  MHT 
Korodense  Tubes  Under  Vacuum 
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vs.  Vw  for  Wire-Wrapped  MHT  Korodense  Tubes 


Under  Vacuum  with  HEATEX  Insert 
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Figure  12.   hQ  vs.  Tcf  for  Wire-Wrapped  MHT  Korodense  Tubes 
Under  Vacuum 
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Figure  13.   hQ  vs.  Tcf  for  Wire-Wrapped  MHT  Korodense  Tubes 
Under  Vacuum 
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an  enhancement  ratio  of  0.88  to  the  0.4  mm  wire -wrap 
(Pw/Dw=20)  with  an  enhancement  ratio  of  1.08.  Figure  13 
completes  the  data  and  shows  the  trend  returning  to  the  plain 
MHT  tube  results,  starting  with  the  0.4  mm  wire- wrap  and 
ending  with  the  0.1  mm  wire-wrap  (Pw/Dw=80).  Two  separate 
figures  were  used  to  help  eliminate  the  clutter  from  the 
overlapping  data.  Table  IX  contains  the  averaged  C^  and  a 
values  for  all  the  conditions  tested  (vacuum  and  atmospheric 
pressure,  with  and  without  the  HEATEX  insert)  for  the  plain 
MHT  tube  and  all  the  wire-wrapped  MHT  tubes.  Two  to  four 
tests  were  conducted  at  each  condition  to  ensure  repeatability 
of  the  data. 

The  wire-wrapped  MHT  tube  data  were  also  checked  to  see  if 
there  was  a  relationship  between  the  fractional  coverage  and 
the  pitch  to  wire  diameter  ratio  to  the  enhancement  over  the 
plain  MHT  tube  (ETc)  .  Figure  14  shows  the  relationship 
between  the  enhancement  ratio  and  the  pitch  to  wire  diameter 
ratio  and  Figure  15  shows  enhancement  ratio  versus  fractional 
coverage.  Data  obtained  by  O'Keefe  [Ref .  18]  for  wire-wrapped 
LPD  Korodense  tubes  are  also  included  in  the  figures. 
However,  the  number  of  his  tests  are  limited  to  only  3  wire 
diameters  (0.5mm,  1.0mm  and  1.6mm) .  The  data  clearly  indicate 
that  wrapping  with  a  large  diameter  wire  causes  a  degradation 
in  performance.  As  the  wire  diameter  is  decreased,  the 
enhancement  passes  through  a  maximum  of  about  8%, 
corresponding  to  a  Pw/Dw<=20  and  a  fractional  coverage  of  «5%. 
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Table  IX.   WIRE -WRAPPED  AND  PLAIN  MHT  KORODENSE  TUBES 


Atmospheric  Pressure 

Tube 

No  Insert 

Heatex  Insert 

Ct 

01 

c< 

a 

Plain 

2  .51 

0.80 

3  .32 

0.79 

1.6mm  wire-wrap 

2.31 

0.72 

2.92 

0.73 

1.0mm  wire-wrap 

2.35 

0.78 

3  .07 

0.79 

0.5mm  wire-wrap 

2.50 

0.85 

3  .19 

0.85 

0.4mm  wire-wrap 

2.53 

0.84 

3  .38 

0.85 

0.3mm  wire-wrap 

2.52 

0.84 

3.27 

0.85 

0.2mm  wire-wrap 

2.53 

0.82 

3.33 

0.82 

0.1mm  wire-wrap 

2.50 

0.81 

3  .36 

0.82 

Vacuum  Pressure 

Plain 

2.39 

0.77 

3.05 

0.78 

1 . 6mm  wire  -  wrap 

2.17 

0.67 

2.84 

0.69 

1.0mm  wire-wrap 

2.20 

0.72 

2.99 

0.73 

0 . 5mm  wire  -  wrap 

2.34 

0.82 

3.05 

0.81 

0.4mm  w^re-wrap 

2.31 

0.83 

3.03 

0.84 

0.3mm  wire-wrap 

2.31 

0.82 

3.14 

0.83 

0.2mm  wire -wrap 

2.39 

0.81 

3.16 

0.80 

0 . 1mm  wi  re - wrap 

2.42 

0.78 

3.23 

0.80 

As  the  wire  diameter  is  further  decreased,  the  enhancement 
ratio  returns  to  unity,  i.e.  back  to  the  plain  tube 
performance.  The  results  from  O'Keefe  [Ref.  18]  for  the  LPD 
Korodense  tube  definitely  show  the  same  trends  as  the  MHT 
tube,  but  the  results  are  somewhat  inconclusive  due  to  lack  of 
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Figure  14.   Enhancement  Ratio  vs.  Pw/Dw  for  Wire-Wrapped 
Korodense  Tubes 
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Figure  15.   Enhancement  Ratio  vs.  Fractional  Coverage  for 
Wire -Wrapped  Korodense  Tubes 
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data.   His  smallest  wire  diameter  was  only  0.5  mm. 

The  enhancement  results  obtained  for  the  wire-wrapped  MHT 
Korodense  tube  were  also  compared  to  a  wire-wrapped  smooth 
tube.  Marto  et  al .  [Ref.  31]  reported  an  optimum  wire  pitch 
to  diameter  ratio  ranging  from  5  to  7  for  wire-wrapped  copper 
tubes.  Sethumadhavan  and  Rao  [Ref.  34]  realized  a  maximum 
enhancement  at  a  fractional  coverage  of  21%,  corresponding  to 
a  Pw/Dw«5 .  This  value  for  Pw/Dw  is  much  smaller  than  the 
optimum  pitch  to  diameter  ratio  of  20  obtained  in  this  study. 

To  further  understand  this  relationship,  data  were 
compiled  for  wire-wrapped  smooth  tubes  studied  at  NPS  and 
compared  with  the  current  data  in  Figures  16  and  17.  Both 
sets  of  data  were  plotted  with  constant  pitch  curves  for  a 
better  comparison  with  the  constant  pitch  results  of  the 
corrugated  tubes.  It  can  be  seen  that  the  wire-wrapped  smooth 
copper  and  wire -wrapped  smooth  titanium  tubes  with  Pw=4  mm 
have  roughly  the  same  slope.  The  same  trend  can  be  seen  for 
the  wire -wrapped  smooth  titanium  and  the  wire -wrapped  LPD  and 
MHT  tubes  with  a  pitch  of  about  7  mm.  Due  to  the  limited  data 
it  is  difficult  to  predict  the  optimal  fractional  coverage  of 
the  wire-wrapped  smooth  tubes  for  a  constant  wire  pitch. 
However,  it  does  appear  to  be  less  than  the  21%  obtained  by 
Sethumadhavan  and  Rao  [Ref.  34]  .  It  can  also  be  seen  that  the 
magnitude  of  the  enhancements  tend  to  vary  with  tube  material 
and  tube  geometry.  Mitrou's  data  [Ref.  14]  for  wire -wrapped 
smooth  copper  tubes  show  results  45  to  55%  higher  than  the 
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Figure  16.   Enhancement  Ratio  vs.  Pw/Dw  for  Wire-Wrapped 
Korodense  and  Wire -Wrapped  Smooth  Tubes,  NPS  Data 
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Figure  17.   Enhancement  Ratio  vs.  Fractional  Coverage  for 
Wire-Wrapped  Korodense  and  Wire-Wrapped  Smooth  Tubes,  NPS  Data 
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wire-wrapped  smooth  titanium  tube  results  for  an  equal  pitch 
(Pw=4  mm)  obtained  by  O'Keefe  [Ref .  18] .  This  indicates  that 
tubes  made  of  materials  with  high  thermal  conductivities  (such 
as  copper)  may  provide  additional  enhancement  characteristics. 
The  effect  of  tube  geometry  can  also  be  seen  in  the  comparison 
of  the  wire-wrapped  smooth  titanium  tube  and  the  wire-wrapped 
LPD  and  MHT  tubes  (all  with  pitches  =7  mm) .  The  wire-wrapped 
smooth  titanium  tube  achieved  an  enhancement  -10%  higher  than 
the  corrugated  tubes  of  similar  pitch.  Figures  18  and  19  are 
identical  to  Figures  16  and  17  with  the  addition  of  the  data 
obtained  by  Sethumadhavan  and  Rao  [Ref.  34]  .  The  tube 
material  used  in  their  tests  is  not  known,  but  appears  to  have 
a  thermal  conductivity  between  copper  and  titanium  (such  as 
copper/nickel  or  aluminum) .  The  most  noticeable  differences 
in  their  data  and  the  data  from  research  at  NPS  are  the 
fractional  coverage  of  the  apparent  maximum  enhancement  and 
the  opposite  slope  of  the  constant  pitch  curves.  This  is  most 
likely  due  the  wire  material  used  for  wrapping  their  tubes. 
They  used  copper  wire  and  all  the  NPS  tubes  were  manufactured 
with  titanium  wire.  The  use  of  a  highly  conductive  material 
may  have  produced  more  of  a  'fin'  effect  with  the  wire  and 
caused  the  peak  enhancement  to  shift.  Further  study  on  the 
effect  of  the  wire  material  is  required  to  determine  the 
effects . 

Possible  effects  on  the  inside  heat  transfer  coefficient 
due  to  wire -wrapping  both  smooth  and  corrugated  tubes  were 
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Figure  18.   Enhancement  Ratio  vs.  Pw/Dw  for  Wire -Wrapped 
Korodense  and  Wire -Wrapped  Smooth  Tubes 
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Figure  19.   Enhancement  Ratio  vs.  Fractional  Coverage  for 
Wire -Wrapped  Korodense  and  Wire -Wrapped  Smooth  Tubes 
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also  investigated.  Figure  20  shows  a  possible  decreasing 
trend  in  C^  with  increasing  fractional  coverage.  However,  the 
magnitude  of  the  decrease  is  minimal  and  the  overall  variation 
in  C^  is  within  the  uncertainty  of  the  plain  tube  results. 

In  summary,  the  enhancement  of  the  outside  heat  transfer 
coefficient  for  a  titanium  MHT  Korodense  tube  over  a  smooth 
titanium  tube  was  close  to  unity.  The  primary  enhancement  in 
the  overall  heat  transfer  coefficient  was  found  to  be  due  to 
an  increase  in  the  inside  heat  transfer  coefficient,  21  to  26% 
above  the  LPD  Korodense  tube  coefficient  and  110  to  117%  above 
the  smooth  tube  coefficient.  The  effect  of  wire-wrapping 
corrugated  tubes  provides  minimal  additional  enhancement  to 
the  outside  heat  transfer  coefficient,  a  maximum  of  =8%  with 
a  Pw/Dw<=20  and  a  corresponding  fractional  coverage  of  5%. 
However,  the  addition  of  wire-wrap  into  the  grooves  of  a 
corrugated  tube  may  have  greater  benefit  when  used  in  a  tube 
bundle  arrangement  due  to  improved  condensate  drainage. 
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Figure  20.    Ct   vs.  Fractional  Coverage  for  Wire-Wrapped 
Korodense  and  Wire -Wrapped  Smooth  Tubes  at  Atmospheric 
Pressure 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  Enhancement  of  the  outside  heat  transfer  coefficient  due 
to  the  corrugation  of  the  MHT  Korodense  titanium  tube  is 
near  unity  when  compared  to  smooth  tube  results. 

2.  The  corrugation  of  the  MHT  Korodense  tube  produces 
enhancements  in  the  inside  heat  transfer  coefficient 
from  110%  to  117%  above  smooth  tube  results,  and  21%  to 
26%  above  LPD  Korodense  tube  results. 

3.  Wire-wrapping  of  MHT  Korodense  tubes  yields  minimal 
additional  enhancement,  «=8%,  when  used  in  a  single  tube 
arrangement . 

4.  The  optimal  fractional  coverage  for  the  wire-wrapping  of 
MHT  Korodense  tubes  is  =5%  (Pw/Dw  ~2  0)  when  used  in  a 
single  tube  arrangement. 

5.  The  thermal  conductivity  of  the  tube  material  effects 
the  enhancement  characteristics  of  wire-wrapped  tubes. 

6.  Some  relationship  appears  to  exist  between  fractional 
coverage  and  the  inside  heat  transfer  coefficient. 

7.  The  modified  Wilson  plot  data  reduction  technique 
compensates  for  vapor  velocity  effects  when  using  the 
Nusselt  theory  in  determining  the  outside  heat  transfer 
coefficient . 

B .  RECOMMENDATIONS 

1.  Manufacture  and  test  additional  wire-wrapped  LPD 
Korodense  tubes  with  smaller  wire  diameters  ranging  from 
0.1  mm  to  0.4  mm  to  validate  the  trends  found  for  the 
MHT  Korodense  tube. 

2.  Manufacture  and  test  additional  constant  pitch  wire- 
wrapped  copper  smooth  tubes  with  wire  diameters  greater 
than  1.6  mm  and  less  than  0.5  mm  and  determine  the 
optimal  wire  diameter  for  specific  wire  pitches. 

3.  Test  a  set  of  wire-wrapped  corrugated  copper  or 
copper/nickel  tubes  to  determine  the  possible  effect  of 
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tube  material  on  condensation,  when  compared  to  the 
titanium  corrugated  tubes. 

4.  Develop  a  numerical  model  to  approximate  the  outside 
heat  transfer  characteristics  for  corrugated  and  wire- 
wrapped  corrugated  tubes. 

5 .  Fabricate  a  set  of  corrugated  tubes  with  outside 
electrical  heater  elements  and  conduct  tests  to  validate 
the  accuracy  of  the  inside  heat  transfer  correlations 
available  for  processing  data  on  corrugated  tubes  with 
and  without  a  HEATEX  insert. 

6.  Input  the  Wolverine  inside  heat  transfer  correlation  to 
the  existing  program.   Reprocess  selected  corrugation 
tube  data  and  compare  results  with  those  obtained  by  the 
Petukhov- Popov  correlation. 

7.  Test  wire-wrapped  corrugated  tubes  in  a  tube  bundle 
arrangement.   Determine  if  additional  overall  bundle 
enhancement  is  possible  with  inundation  effects 
included,  and  if  the  optimal  fractional  coverage  shifts 
from  that  found  for  single  tube  tests. 

8.  Fabricate  and  test  a  set  of  multistart  corrugated  tubes 
with  and  without  wire-wrapping  to  determine  if 
additional  enhancements  can  be  realized  using  the 
multistart  corrugation  geometry. 

9.  Conduct  additional  tests  on  wire-wrapped  smooth  and 
corrugated  tubes  to  determine  if  a  relationship  exists 
between  fractional  coverage  and  the  inside  heat  transfer 
coefficient . 

10.  Test  the  enhancement  effects  of  various  wire  materials 
when  used  for  wire -wrapping  smooth  and  corrugated  tubes. 

11.  Investigate  the  possibility  of  bonding  the  wire-wrap  to 
the  tube  surface  and  test  the  enhancement  effects  for 
smooth  and  corrugated  tubes. 
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APPENDIX  A.   SYSTEM  CORRECTIONS 

A.   FRICTIONAL  TEMPERATURE  CORRECTIONS 

Coolant  flowing  through  a  tube  experiences  a  temperature 
increase  due  to  frictional  heating.  The  magnitude  of  the 
temperature  rise  depends  on  the  internal  tube  geometry  and  the 
velocity  of  the  coolant.  With  the  cooling  water  sump  and  the 
test  section  at  ambient  temperature,  measurements  were  taken 
to  determine  the  temperature  increase  due  to  frictional 
effects  at  flow  rates  ranging  form  10%  to  100%  on  the 
rotameter.  Two  separate  runs  were  conducted  for  the  MHT 
Korodense  tube,  one  with  the  HEATEX  insert  installed  and  the 
second  with  no  insert.  Figure  A.  1  show  the  results  for  the 
two  runs.  The  data  was  curve  fitted  to  a  third  order 
polynomial  as  shown  in  Table  A.  1. 
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Figure  A.l   Frictional  Temperature  Rise  for  the  MHT 
Korodense  Titanium  Tube  with  and  without  a  HEATEX  Insert 
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Table  A.  1   FRICTIONAL  TEMPERATURE  RISE  EQUATIONS  FOR  MHT 
KORODENSE  TUBE 


Inside 
Condition 

Polynomial  Equation 

HEATEX 
Insert 

Trise  =    2.564X10"4  +  6.263xlO"4Vw 

+  2.603xl0"3 (Vw) 2  +  7.830x10" 
6  (V  )  3 

No 
Insert 

Trise  =    6.790X10"5  +  1.632xl0~3Vw 

+  8.409X10"4 (VJ2  +  1.111x10" 
4(V  )3 

where:   Trise  =  Temperature  Rise  (K) 
Vw   =   Fluid  Velocity  (m/s) 
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APPENDIX  B.   UNCERTAINTY  ANALYSIS 

All  derived  quantities  have  an  inherent  uncertainty- 
associated  with  them  when  their  results  are  obtained  from 
measured  values.  The  uncertainty  arises  due  to  inaccuracies 
in  the  measurement  of  the  physical  quantities  used  in  their 
calculations.  The  difference  between  the  actual  value  and  the 
measured  value  depends  on  the  accuracy  of  the  measuring 
device,  the  calibration  of  the  device,  and  the  experience  of 
the  operator.  Eventhough  the  error  associated  with  any  single 
measurement  may  seem  insignificant,  it  may  have  a  profound 
effect  on  the  accuracy  of  the  derived  quantity  due  to  the 
calculations  made. 

Kline  and  McClintock  [Ref.  67]  developed  a  method  to 
estimate  the  uncertainty  in  the  derived  quantity,  R,  which  is 
a  function  of  several  measured  variables,  xx ,  x2 ,  -..,  xn,  as 
follows : 


W    = 


ax/1 


ax,  2 


dx„    n 


1/2 


(B.l) 


where : 

WR  is  the  uncertainty  of  the  desired  dependent  variable 
xx,  x2 ,  ...,  x^  are  the  measured  independent  variables 


wi >      W2  > 
variables 


Wn  are  the  uncertainties  in  the  measured 
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Georgiadis  [Ref .  68]  presents  a  complete  discussion  on  the 
uncertainty  analysis  used  for  this  experiment.  Based  on 
Georgiadis'  developments,  Mitrou  [Ref.  14]  wrote  a  program  for 
calculating  the  uncertainties  in  the  mathematically  derived 
quantities  obtained  from  data  collected  using  the  apparatus 
utilized  in  this  research.  Modifications  were  made  to  the 
program  which  allow  for  the  data  collected  on  MHT  Korodense 
tubes  to  be  processed.  Sample  uncertainty  results  for  the 
following  data  collection  conditions  are  presented  in  this 
appendix. 

1.  Vacuum  pressure  with  no  insert  and  high  flow  rate  (3.83 
m/s)  . 

2.  Vacuum  pressure  with  no  insert  and  low  flow  rate  (1.03 
m/s)  . 

3.  Vacuum  pressure  with  HEATEX  insert  and  high  flow  rate 

(3.83  m/s) . 

4.  Vacuum  pressure  with  HEATEX  insert  and  low  flow  rate 

(1.02  m/s) . 

5.  Atmospheric  pressure  with  no  insert  and  high  flow  rate 

(3 .82  m/s) . 

6.  Atmospheric  pressure  with  no  insert  and  low  flow  rate 

(1.02  m/s) . 

7.  Atmospheric  pressure  with  HEATEX  insert  and  high  flow 
rate  (3 . 82  m/s) . 

8.  Atmospheric  pressure  with  HEATEX  insert  and  low  flow  rate 

(1.02  m/s) . 
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1.  DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  LMNTNU8 

Pressure  Condition:    Vacuum 

Uapor  Temperature  =  48.574         Deg  C 

Water  Flow  Pate  (%)  =  30.00 

Water  Velocity  =»  3.83        'm/s> 

Heat  Flu...  =  I.48IE  +  05   'W/m  2) 

Tube-metal  thermal  conduc.  =  21.0         (  W  /  m  .  K  ) 

Petukov -Popov  constant           =  2.3857 


UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 

Mass  Flow  Pate  ,  Md  0.80 

Reynolds  Number,  Re  1.09 

Heat  Flux  ,  q  1  . 37 

Log-Mean-Tem  Diff,  LMTD  1.63 

Wail  Resistance,  Ru  4.37 

Overall  H.T.C. ,  Uo  2.48 

Water-Side  H.T.C.  ,  Hi  7.47 

Vapor-Side  H.T.C,  Ho  6.77 

2.  DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  LMNTNV8 

Pressure  Condition:    Vacuum 

Vapor  Temperature  = 

Water  Flow  Rate  (  7.  ) 

Water  Velocity  = 

Heat  Flux 

Tube-metal  thermal  conduc.     = 

Petukov -Popov  constant 

UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 

Mass  Flow  Rate,  Md  2.99 

Reynolds  Number,  Re  3.08 

Heat  Flux  ,  q  3.07 

Log-Mean-Tem  Diff,  LMTD  .56 

Wall  Resistance,  Rui  4.97 

Overall  H.T.C.  ,  Uo  3.13 

Water-Side  H.T.C.  ,  Hi  7.81 

vapor-Side  H.T.C,  Ho  13.98 


48.691 

(Deg  C 

20.00 

1  .03 

(m/s  ) 

1  .  149E+05 

( W/m"2  ) 

21  .0 

(  W  /  m  .  K  ) 

2.3857 
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3.  DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  LMNTHV8 

Pressure  Condition:    Vacuum 

Vapor  Temperature  =  48.513        ( Deg  C) 

Water  Flow  Rate  (X)  =  30.00 

Water  Veioc  ity  =  3 . ~3         (m/s  ) 

Heat  Flu>  =  1 .526E+0S   C  U / m " 2  ) 

"jte-^etai  tr.ernal  c  o  n  d  u  c  .  =  21. 0         (  W  /  m  .  K  ) 

Petuk  ov-Popov  constant  =    3.0112 


UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 


Mass  p  1  c  w  Pate,  M  d 
Reynolds  Number ,  Pe 
Heat  Flu-  ,  q 
Log-Mean-Tem  Diff,  LMTD 
Wall  Resistance,  Rw 
Overall  H.T.C. ,  Uo 
Water-Side  H.T.C. ,  Hi 
V a p o r -  5 i d e  H.T.C.,  Ho 


4.  DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  LMNTHV8 

Pressure  Condition:    Vacuum 

Vapor  Temperature  =    48.551        (Deg  C) 

Water  Flow  Rate  (X)  =    20.00 

Water  Velocity  =    1.02         (m/s) 

Heat    Flu.-  =  1.202E  +  05       <W/m'2> 

Tube-metal  thermal  conduo.  =    21.0         (W/m.K) 


0 

.30 

1 

.08 

1 

.33 

1 

.53 

4 

.97 

2 

.42 

5 

.92 

6 

.35 

p, 


tul ov-Popov  constant  =    3.01 


UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 

Mass  Flow  Rate  ,  Md  3.00 

Reynolds  Numoer  ,  Re  3.09 

Heat  Flux  ,  q  3.08 

Log-Mean-Tem  Diff,  LMTD  .54 

Wall  Resistance,  Rw  4.97 

Overall  H.T.C.  ,  Uo  3.12 

Water-Side  H.T.C.  ,  Hi  6. 35 

Vapor-Side  H . T . C .  ,  Ho  1  1  .54 
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5.  DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 


File  Name:  LMNTNA9 

Pressure  Condition:    Atmospheric 

'Japcr  Temperature 

Water  Flow  Pate  <  %  ) 

Water  Velocity 

Meat  Flu.-. 

Tube-metal  thermal  c endue.     = 

Petukov-Popov  constant 


,1(31  I  Pa  ) 

99.954  Deg  C  ) 

80.  00 

3  .  9  Z  m  /  s  > 

3.959E+05  <w/m  :;< 

21.0  ( .  W / m . K  ' 
Z.4S59 


UNCERTAINTY  ANALYSIS: 


VARIAELE 


PERCENT  UNCERTAINTY 


Mass  Flow  Pate  ,  Md 
Reynolds  N  u  m  b  e  r  ,  P  e 
Heat  F 1 u  a  ,  q 
I_  :g-Mean- rem  Diff,  LMTD 
Wa  i  1  ces i 5 tance ,  Pw 
0 vera  11  H.T.C.,  Uo 
'wat  er- 5  i  be  H  .  T  .  C  .  ,  Hi 
Uapor-5ide  H . T . C . ,  Ho 


0.S0 
1.11 
1.10 
.5! 
4  .  97 
'  .2G 

".Eg 

4.97 


6.  DATA  FOR  THE  UNCERTAINTY  ANALY9IS 


File  Name:  LMNTNA8 

Pressure  Condition:    Atmospheric  (101  hPa  > 

Vapor  Temperature  =    99.981        ( Deg  C) 

Water  Flow  Pate  <  ".  )  =    20.00 

Water  Velocity  =    1.02         (m/s) 

Heat    Flux  =         3.M4E  +  05       (W/m*2) 

Tube-metal  thermal  conduc.     =    21.0         (W/m.K) 

Petukov-Popov  constant  =    2.4659 


UNCERTAINTY  ANALYSIS 
VARIABLE 


PERCENT  UNCERTAINTY 


Mass  F low  Rate  ,  Md 
Reynolds  Number,  Re 
Heat  F  lux  ,  q 
Log-Mean-Tem  Dif f  ,  LMTD 
Wall  Resistance,  Rw 
Overall  H.T.C.  ,  Uo 
Water-Side  H.T.C,  Hi 
Vapor-Side  H.T.C.,  Ho 


13 


00 
10 
04 
21 
97 
05 
04 
03 
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0 

.30 

] 

.  1  1 

1 

.09 

.53 

4 

.97 

1 

.24 

5 

.63 

4 

.25 

7.  DATA  FOR  The  UNCERTAINTY  ANALYSIS: 

File  Name:  LMNTHAB 

Pressure  Condition:    Atmospheric  M01  kPa  ) 

Vapor  Temperature             =  99.324        (  Deg  C 

Water  Flow  Rate  ('/.)                              =  80.00 

Water  Velocity                =  3.32        (m/s) 

Heat  Flu/.                      =  4.083E+0S   (W/m'2) 

Tube-metal  thermal  conduc.     =  21.0        (W/m.K) 

F  e  t  u  y   o  v  -Popov  constant           =  3.3112 


UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 


Mass  F lew  Ra t e  ,  Md 
Reynolds  Number,  Re 
Heat  Flu...,  q 
Log-Mean-Tem  Di f f ,  LMTD 
wall  Resistance,  R ui 
Overall  H. T.C .  ,  Uo 
Water-Side  H. T.C .  ,  Hi 
Vapor-Side  H.T.C.,  Ho 


8.  DATA  FOR  THE  UNCERTAINTY  ANALYSIS: 

File  Name:  LMNTHA8 

Pressure  Condition:    Atmospneric  (101  k  P  a  ) 

Vapor  Temperature  =    39.999        (Deg  C) 

Water  Flow  Rate  (7.)  =    20.00 

Water  Velocity  =    1  .02         (m/s  ) 

Heat  Flu.-.  =    3.367E  +  05   (W/m'2  > 

Tube-metal  thermal  conduc.     =    21.0         (W/m.K) 

Petut-ov -Popov    constant  =         3.3112 


UNCERTAINTY  ANALYSIS: 

VARIABLE  PERCENT  UNCERTAINTY 

Mass  Flow  Pate,  Md  3.00 

Reynolds  Number,  Re  3.11 

Heat  Flux  ,  q  3.04 

Log-Mean-Tem  Diff,  LMTD  .19 

Wall  Resistance  ,  Rw  4.97 

Overall  H.T.C. ,  Uo  3.05 

Water-Bide  H.T.C.  ,  Hi  6.14 

Vapor-Side  H.T.C.  ,  Ho  9.73 
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APPENDIX  C.   DRPML  PROGRAM  LISTING 

The  data  collection  and  reduction  program,  DRPML,  which 
was  used  throughout  this  study  follows: 


90 


1000 
I  00S 
1007 

:  003 
1009 
1010 


DRPML  .LONG) 
REVISED  FROM  ORPOI 


SEP  '992  (O'KEEFE) 


3E  USED  WITH  NGN- INSTRUMENTED  TUEES  ONLY 
TAKES  DATA  IN  THE  FORMAT  OF  SUENSEN/O ' KEEFE/LONG 
CAN  REPROCESS  ANY  NON-INSTRUMENTED  DATA 


1  3 

1  4 
15 


1018 

101  9 

1020 

1021 

1022 

1023 

1024 

1025 

1025 

1027 

102S 

1029 

1030 

1031 

1032 

1033 

1034 

1035 

103E 

1037 

1033 

1039 

1040 

104  I 

1042 

104  3 

cu  ,Re 

1044 

1045 

1047 

1  048 

1049 

1050 

1051 


THIS  PROGRAM  WAS  USED  TO  COLLECT  ALL  THE  NON- 
INSTRUMENTED  DATA  TAKEN  BY  LONG  (JAN-JUN  199: 

MEANING  OF  ALL  FLAGS  IN  PROGRAM 


FOR  TITANIUM  TUBES 


IFT:    FLUID  TYPE 

ISO:    OPTION  WITHIN  PROGRAM 

IM:     INPUT  MODE 

IUIL:   VALUE  OF  Ci    USED 

IFG:    FINNED  OR  SMOOTH 

INN:    INSERT  TYPE 

IWT:    LOOP  NO.  WITHIN  PROGRAM 

IWTH:   ALTERNATIVE  CONDENSER  TUBES 

IMC:    TUBE  MATERIAL 

ITDS:   TUBE  DIAMETER 

IPC:    PRESSURE  CONDITION 

INF:    DIMENSIONLESS  FILE  REQUIRED 

IFF:    FLOT  FILE  REQUIRED 

IOU:    OUTPUT  REQUIRED 

IHI:    INSIDE  HTC  CORRELATION 

IOC:    OUTSIDE  HTC  THEORY /CORRELATION 
COM  /Cc/  C(7) 
COM  /Cc55/  T55(5) 
COM  /Cc5G/  T5G(5  ) 
COM  /Cc57/  T57( 5  ) 
COM  -'Cc58/  T58<5) 
COM  /Fid/  If t  ,  Istu 
DIM  Emf ( 20  )  ,Tw< 6  ) 

COM    /Pr/    Qpa< 42  )  ,Tfm(42  )  ,Tfitr  ,Ipc  ,Qpr 

COM    /Wi  1/    Nrun.Itm.Iuith.Imc.Ife.Ijob.IuidJfg.Ipco.Ifto.Iujil  ,  Ihi  ,  loc  ,  Inam  ,K 
p  ,  R  m  ,  A  x 

COM  /Geom/  D1  ,02  ,Di  ,Do  ,L  ,L1  ,L2 
COM  /Fric/  I s tuo  ,  Inn  ,  Iwt ho  ,Vu  ,Vws , Inamo 
DATA  0.10086091  ,25727.94369,-767345.8295,78025595.81 
DATA  -924748B589  ,B.97638E+1 1  ,-2 . 66 1 92E+ 1 3  ,3.94078E+14 
READ  C( *  ) 

DATA  273. 15.2.5943E-2  ,-7.267lE-7  ,3.2941E-1 1  .-9.7719E-16  .9.7121E-20 
READ  T55( *  ) 
DATA  273. 15 .2.5378E-2  .-5.98S3E-7  .-3.1242E-11  ,1  .3275E-14  ,-1  .0188E-18 
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053  FEAD  T56( *  ) 

054  DATA  273. 15  .2.5923E-2  .-"7.3933E-7  .2.36  2SE-1  1  ,1  . 97 1 7E- 1 5  ,-2 . 2  486E- 1 9 

055  l5EAD  T57' *  ) 

05G  DATA  273.  15  ,2. 5931 E-2  , -7.S232E-7  .4.0657E-1  1  ,-l  .279  IE- I  5  .6.4402E-20 

057  READ  T5S( ♦  ) 

058  Dr=  . 015575    '  Outside  diameter  of  the  outlet  end 

059  Q5sp=.1524    '  Inside  diameter  of    stainless  steel  test  section 
050  Ax=PI*Ds5p  2/4 

0G 1  Alp2=0. 

052  L=. 13335      I  Condensing  length 

063  Lt=.0G0325    !  Inlet  end  "fin  length" 

064  l_2  =  . 034925    '  Outlet  end  "fin  length" 
0G5  PRINTER  IS  I 

06G  BEEP 

067  PRINT  USING  " 4X , ""Select  option: 

069  PRINT  USING  "EX."   0  J a^ e    data  or  re-process  previous  data 

09O  PRINT  USING  "GX,""  2  WILSON  Analysys 

06  INPUT  I  so 

08  Iso=Iso+l 

11  IF  Iso. 1  THEN  3094 

12  EEEF 

15  INPUT  "SELECT  FLUID  ( 0=UATER ,  I=R-1I3,  2=EG)",Ift 

16  Ifto-Ift 

17  BEEP 

18  ljob=0 

19  INPUT  "ENTER  INPUT  MODE  <  0  =  3054A  ,  1  =FILE  ) "  ,  In 

20  Im=Im+ 1 

23  BEEP 

24  IF  Im=1  THEN 

26      INPUT  "ENTER  MONTH,  DATE  AND  TIME  (  MM  :  DO  :  HH  :  MM  :  SS  ) "  ,Date$ 
29      OUTPUT  709; "TD" ;0ate$ 

32  OUTPUT  ^09; "TD" 

33  ENTER  709;Date$ 
Z5  END  IF 

Z5  IF  Ijob-I  THEN 

33      BEEP 

41      INPUT  "SKIP  PAGE  AND  HIT  ENTER"  ,0k 

44  END  IF 

45  PRINTER  IS  701 

46  IF  Im«1  THEN 

48      ENTER  "OSiDateS 

50  PRINT  "  Month,  date  and  tine  :";Date$ 

51  END  IF 
53  PRINT 

56  PRINT  USING  "10X,""NOTE:  Program  name  :  DRPML 

71  IF  I job=l  THEN  1 189 
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NAME  FOR  THE  RAW  DATA  FILE" 

,0_f ile$ 

IGX  ,"  "Fi  le  name     :  " "  , 1 4A" 

;D_fi le$ 

_f  iie$  ,32 

TO  D  file$ 

174  SEEP 

186  INPUT  "SELECT  ( Ci : 0=FIND  ,  1 =ST0RED  Cx)"  .Iwxl 

139  IF  1,^=1  THEN 

192  BEEP 

135  INPUT  "GIVE 

1  98  PRINT  USING 

201  CREATE  BOAT  D_ 

Z04  ASSIGN  @Fxle 

207  BEEP 

210  INPUT  "ENTER  GEOMETRY  CODE  (  1 -FINNED  ,0-PLAIN >"  ,1  fg 

21  1  lnn  =  0 

212  PRINTER  IS  1 

216  BEEP 

217  PPINT  "    ENTER  INSERT  TYPE:" 

213  PRINT  "  0=NONE  (DEFAULT)" 

219  FRINT  "  1 =  T  i J  I  5  T  E  D  TAPE" 

220  PRINT  "  2  =  UIIRE  WRAP" 
22'  FRINT  "       3=HEATEX" 
222  INPUT  Inn 

225  OUTPUT  @Fiie; Ifg  ,Inn 

227  lwt=0   '  FOR  UNINSTRUMENTED  TUBE 

228  Fh  =  0 
231  Fp=0 

234  Fu=0 

235  lstu=0 

236  l5tuo=l5tu 

238  IF  lfg=0  THEN  1241 

239  INPUT  "FIN  PITCH,  HEIGHT  AND  WIDTH  IN  MM,  Fp  ,Fh  ,Fuj"  ,Fp  ,Fh  ,Fuj 

241  OUTPUT  SFile; Iut  ,Fp ,Fw,Fh 

242  ELSE 

249  BEEP 

250  PRINTER  IS  1 

252  PPINT  "   STUDENT'S  DATA  TO  BE  REPROCESSED:" 

253  PRINT  "   0=SWENSEN/O'KEEFE/LONG  (DEFAULT)" 

254  PRINT  "   1=UAN  PETTEN/MITROU/COUMES/GUTTENDORF " 

255  INPUT  Istu 

256  I  s t uo=  Istu 
258  BEEP 
259 
260 
2B1 
262 
263 
264 
265 
266 


PRINT 

IF  Istu=l 

THEN 

PRINT 

"  STUDENT  NAME: 

PRINT 

"  0  =  UAN  PETTEN" 

PRINT 

"  1=MITROU" 

PRINT 

"  2=COUMES" 

PRINT 

"  3=GUTTEND0RF" 

ELSE 
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I2G7 
1  268 
1269 

1271 


1  275 
1  275 
1277 
1278 

EN 

1279 

1280 

1281 

RF 

1232 

1283 

1284 
1285 
12SG 
1287 
1288 
1289 
1290 
1291 
1292 
1293 
1294 
1295 
129G 
1297 
12  98 
1299 
1300 
1301 
1302 
1303 
1306 
1320 
1321 
1323 
1324 


PRINT  "  4=SUENSEN" 

PRINT  "  5=0'KEEFE" 

PRINT  "  6=L0NG" 
END  IF 
INPUT  I  nam 
I namo  = 1  nam 
BEEP 

INPUT  "GIVE  ~HE  NAME  OF  THE  EXISTING  DATA  F ILE "  ,D_f i ieS 
PRINTER  IS  701 
IF  lnam=0  THEN  FRINT  USING  " 16X 


Data  taken  by 
Data  taken  by 
Data  taken  by 
Data  taken  by 


IF  Inam=1  THEN  PRINT  USING  "1EX 

IF  Inam=2  THEN  PRINT  USING  "!GX 

IF  Inam=3  THEN  PR  ITJT  USING  "  1  6 X 

IF  Inam=4  THEN  ,L  INT  USING  "lGx,""Data  taken  by 

IF  Inam=5  rHEN  PRINT  USING  "1SX,""Data  taken  by 

IF  Inam=6  THEN  PRINT  USING  "l6X,""Data  taken  by 

PRINT  USING  "16X,,",This  analysis  done  on  file  : 

PRINTER  IS  1 

BEEP 

INPUT  "ENTER  NUMBER  OF  DATA  SETS  STORED"  ,Nrun 

ASSIGN  ©File  TO  D_f l le$ 

ENTER  BFilei Ifg  ,Inn 

IF  lstu=0  THEN 

ENTER  ©File;  Iwt  ,Fp  ,Fuj  ,Fh 

ELSE 

IF  Ifg-0  THEN  ENTER  @File;Iwt 

IF  Ifg=1  THEN  ENTER  @Fi le ; Fp ,Fu ,Fh 

END  IF 
END  IF 

IF  [job-1  THEN  1537 
IF  Ift0  THEN  1345 
BEEP 

PRINTER  IS  1 
PRINT  USING  "4X 
PRINT  USING  "6X 
PRINT  USING  "6X 
PRINT  USING  "SX 
PRINT  USING  "6X 
INPUT  Iuith 
Iwt  ho=  Iuit  h 


UAN  PETT' 

MITROU"" 

COUMES"" 


GU" 


:ndo 


SWENSEN' 

O'KEEFE' 
LONG 


0A";D  file$ 


Select  tube  type : 

0  Copper  Tube " " " 

1  Wolverine  Korodense  LPD  Titanium  Tube" 

2  Smooth  Titanium  Tube 

3  Wolverine  Korodense  MHT  Titanium  Tube" 
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I3ZG 

1323 
I  329 

1333 

1331 

1333 

1333 

I  334 

1335 

1336 

133S1 

1  338 

1339 

1340 

1  341 

1342 

1343 

1345 

1  343 

1349 

1351 

1353 

1357 

1360 

1363 

1366 

1367 

1359 

1373 

1375 

1373 

1383 

1394 

1399 

1402 

1405 

1457 

1458 

1459 

1460 

1461 

1462 

1463 

1464 

1466 


6EEP 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


USING  "4X 
USING  'GX 
USING  '5X 
■J  SING  "6a 
USING  '5X 

PRINT  USING  "6X 

PRINT  USING  "6X 

PRINT  USING  "BX 

INPUT  Hyp 

PRINTER  IS  701 

BEEP 

PRINTER  IS  1 

PRINT  USING  "4X 

PRINT  USING  "6X 

PRINT  USING  "6X 

PRINT  USING  "6X 

INPUT  Imc 

PRINTER  IS  1 

BEEP 

Itds-1 

IF  luth=0  THEN 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
INPUT  Itds 

END  IF 

PRINTER  IS  701 

IF  lwth=0  THEN 
Di=.0127 
Do=.01 905 

END  IF 

IF  Iwth=1  THEN 
Di=. 01347 
Do=.  01585 

END  IF 

IF  Iuith  =  2  THEN 
Do-. 01585 
Di=.  01386 

ENO  IF 

IF  Iwth-3  THEN 
Do=. 01607 
Di=. 01353 

END  IF 

D1=.01 905 


Select  tube  Enhancement  used:" 

0  SMOOTH  TUBE 

!  FINNED  TUBE 

2  WIRE-WRAPPED  SMOOTH  TUBE 

3  LPD  KORODENSE  TUBE 

4  WIRE-WRAPPED  LPD  KORODENSE  TUBE' 

5  MHT  KORODENSE  TUBE 

6  WIRE-WRAPPED  MHT  KORODENSE  TUBE' 


Select  Material  Code:""" 

0  Copper     1  Stainless  steel" 

2  Aluminum   3  90:10  Cu-Ni ' 

4  Titanium 


4X  ," "SELECT  TUBE  OIA  TYPE 

6X,""0   SMALL 

1  BX  ,"  "  1   MEDIUM  (DEFAULT  )"  ' 
6X,""2   LARGE 


i  ID  OF  MEDIUM  8,  LARGE  TUBES 
"  OD  OF  MEDIUM  TUBE 
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467 

468 

469 

471 

474 

477 

478 

479 

484 

485 

487 

489 

492 

495 

497 

498 

501 

504 

507 

=  08 

510 

557 

538 

543 

546 

549 

552 

556 

557 

558 

559 

560 

56  1 

563 

563 

564 

565 

567 

558 

569 

570 

571 

573 

573 

574 

575 


D2-. 01587 

IF  ltd 5=0  THEN 
Oo=.0127 
ji=. 009525 

END  IF 

IF  I tds=2 
Iuth-1 
Iwth-1 
Iuth-2 
Iwth-2 


11- 
IF 
If 

IF 
IF 
IF 
IF 
IF 
IF 


THEN  Do=.02S 

OR  Iwth»3  THEN  D!  -=.01535 

OR  Iwth-3  THEN  02=. 01 585 

THEN  01 =.3 1537 

THEN  02=. 01567 
lmc  =  0  THEN  Kcu=335 
Inc-I  THEN  Kcu=lG 
Imc  =  2  THEN  Kcu-167 
Imc  =  3  THEN  Kcu  =  45 
Imc  =  4  THEN  Kcu  =  2 1  .0 
Rm=Qo*LOG( Do/Di  i/(2*Kcu 
BEEP 

INPUT  "ENTER  PRESSURE  CONDITION  :0 
Ipco=  Ipc 
I  n  f  =  0 
BEEF 
Ife-I 

PRINTER  15  701 
PRINT  USING  '■  1  6 X 
PRINT  USING  "16X 
PRINT  USING  " 16X 
PRINT  USING  "  1BX 
Ihi-0 

PRINTER  IS  1 
PRINT 


Wall  resistance  based  on  outside  area 


Ipc 


This  analysis  includes  end-fi 
Thermal  conductivity       = 
Inside  diameter,  Di        = 
Outside  diameter,  Do       = 


n  effect" 

"  " 

"  ,3D.D,"" 

( W/m.K  ) 

"  ;  Kcu 

"  ,  DD.DD," 

"  (nn  ) 

'  ;Di* 1000 

"  .DD.DD," 

(  mm  ) " 

'  ;Do»1000 

PRINT 
PRINT 
PRINT 
INPUT 


SELECT  INSIDE  CORRELATION:" 
0=SIEDER-TATE  (DEFAULT)" 
1=SLEICHER-P0USE" 
2=PETUKH0U-P0P0U" 


Ihi 
IF  lhi=0  THEN 

BEEP 

INPUT  " 
END  IF 
loc=0 
EEEP 
PRINT 
PRINT  " 
PRINT  " 
PRINT  " 
INPUT  Ioc 
EEEP 


SELECT  REYNOLDS  EXPONENT"  ,Re.<p 


SELECT  OUTSIDE  THEORY/CORRELATION  FOR  WILSON  ANALYSIS:" 
0=NUSSELT  THEORY  (DEFAULT)" 
1 =FUJ 1 1  (1979)  CORRELATION" 
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1576   I  tin- I 
157"   PRINT 

3   PRINT  '     SELECT  COOLANT  TEMPERATURE  RISE  MEASUREMENT:" 
1579   IP  l5tu=0  THEN  PRINT  "         0=SINGLE  TEFLON  T/C" 
15S0   PRINT  "         l=QUARTZ  THERMOMETER  'DEFAULT)" 
1581   PRINT  "         2=I0-JUNCTION  THERMOPILE" 
I53Z   INPUT  Itm 
!583   PRINTER  IS  701 

1584  IP  Itm-0  THEN  PRINT  USING  " 1 6X .'"This  analysis  uses  the  SINGLE  TEFLON  T/C 
readings  

1585  IF  Itm—  1  THEN  PRINT  USING  "16X,""Thi5  analysis  uses  the  QUARTZ  THERMOMETER 
readi ngs 

155G   IF  Itm-2  THEN  PRINT  USING  "!6X,""This  analysis  uses  the  10-JUNCTION  THERMO 

PILE  read l ngs 

1537   Iic=1  i  FOP  MODIFIED  WILSON 
1588   IF  I  h  i  =0  THEN  Ci=.027 
1595   IF  Ihi-1  THEN  Ci  =  l . 
159G   IF  Ihi-2  THEN  Ci  =  l  . 
1599   IF  Iwil=1  THEN 
1801      BEEP 

1602  INPUT  "ENTER  Ci  IF  DIFFERENT  FROM  STORED  VALUE"  ,Ci 

1603  END  IF 

1604  PRINTER  IS  701 

1G05   IF  lhi=0  THEN  PRINT  USING  " 1 6X , " "Modi f led  Sieder-Tate  coefficient    =  "  "  , Z 
.4D" ;Ci 

1606  IF  lhi=0  THEN  PRINT  USING  " 1 6X  ,"" Chosen  Reynolds  No.  exponent        =  " "  ,0 
.DD" ;Rexp 

1607  IF  Ihi-1  THEN  FRINT  USING  " 1 6X , " " Mod i f l ed  S leicher-Rouse  coefficient    =  " 
"  ,Z.4D";Ci 

1609   IF  Ihi-2  THEN  FRINT  USING  " 1 6X  , " " Modi f led  Petuk hov-Popov  coefficient    =  " 
"  ,Z.4D"  ;Ci 

1609  IF  lnn  =  0  THEN  PRINT  USING  " 1 6  X  ,  ""Using  no  insert  inside  tube 

1610  IF  Inn=2  THEN  FRINT  USING  "16X, ""Using  wire  wrap  insert  inside  tube 

1611  IF  Inn  =  5  THEN  FRINT  USING  "16X," "Using  HEATEX  insert  inside  tube 

1612  IF  lstu=0  THEN 

1615  IF  Inn=1  THEN  PRINT  USING  "16X, ""Using  twisted  tape  insert  inside  tube" 

1614   ELSE 

1515      IF  Inn=l  THEN  PRINT  USING  "16X, ""Using  wire  wrap  insert  inside  tube 

1616  END  IF 

1619  IF  Iic-0  AND  Ife-1  THEN  Ac=26.4 

1620  IF  Iic-1  THEN  Ac  =  0. 

1621  BEEP 

1622  IF  Ijob-1  THEN  1648 

1623  PRINTER  IS  1 

1624  1  INPUT  "NAME  FOR  A  TEMPORARY  PLOT  PILE  (TO  BE  PURGED  )"  ,P_f i 1 e$ 
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1  625 

1626 
1634 

164a 
1648 
1651 
1654 
1657 
1660 
I  661 
1  666 
1667 
1672 
1673 
1674 

TUBE 
1675 


P_f l ie$="DUMMY" 

5EEP 

CREATE  8DAT  P_f lie*  ,  10 

ASSIGN  ®Fiiep  TO  P_file$ 

IF  Ijob=l  THEN 

I  o  v  =  1 

GOTO  1657 
END  IF 
BEEP 

INPUT  "SELECT  OUTPUT   0=SHQRT  ,  1=L0NG;",Iov 
Iov=Iov+l 
PRINTER  IS  701 

IF  ltyp=0  THEN  PRINT  USING  "16X 
IF  Ityp-1  THEN  PRINT  USING  "15X 
IF  Ityp=2  THEN  PRINT  USING  "16X 


Tube  Enhancement 
Tube  Enhancement 

Tube  Enhancement 


IF    I  t  y  P  : 


THEN  PRINT  USING 


6X  ,  "  Tube  Enhancement 


SMOOTH  TUBE 

FINNED  TUBE 

WIRE-WRAPPED  SMOOTH 

LPD  KOPODENSE  TijEE" 


1678  IF  Ityp=4  THEN  PRINT  USING  "16X,""Tube  Enhancement    :  WIRE-WRAPPED  LPD  KO 
RODENSE  TUBE 

1679  IF  Ityp-5  THEN  PRINT  USING  " 1 6X , " " Tube  Enhancement    :  MHT  KORODENSE  TUBE" 


1660 
PODEN 

i6s: 

1  633 

1684 

1685 

1686 

1687 

1688 

168  9 

15901 

Z.DD" 

1591 

169C 

1693 

1696 

1699 


1712 
1723 
1723 
I  724 

Ts 
17  25 


IF  Ityp=6 

SE  TUBE 

BEEP 
IF  lmc=0 
IF  lmc=1 
IF  Imc=2 
IP  Imc=3 


'HEN  PRINT  USING 


6X,""Tube  Enhancement 


THEN  P»INT  USING 

THEN  PRINT  USING 

THEN  PRINT  USING 

THEN  PRINT  USING 

IF  Inc=4  THEN  PRINT  USING 

IF  lpc=0  THEN  PRINT  USING 

IF  Ipc=1  THEN  PRINT  USING 
PRINT  USING  " 16X  ," "Fin  pi1 
;Fp  ,Fui  ,Fh 

IF  CIuul=0  OR  Iwil-2)  AND  Im=2  THEN 
I job-1 
Iwd=  I 

CALL  Wi lson<Ci  ) 
END  IF 

J=0 

IF    Iov»1     THEN 
PRINT 
IF    Ihi-1     THEN 

PRINT    USING    " I0X 
Re.<p 

PRINT    USING    " 10X 


WIRE-WRAPPED    MHT    KO 


iTEEL"" 


16X  , 

"Tube  material 

COPPER 

16X  , 

"Tube  material 

STAINLESS-S 

16X  , 

"Tube  matenial 

ALUMINUM 

16X  , 

"Tube  material 

90/10  CU/NI 

16X  , 

"Tube  matenial 

TITANIUM 

16X  , 

'  "Pressure  condit 

l  on 

UACUUM 

16X  , 

'"Pressune  condit 

ion 

ATMOSPHERIC 

"1  ,  Ul 

l  dt h  ,  and  he lght 

<-  mm  > 

" "  ,QD.DD  ,2X 

7    nn    "• 


DD  ,2X  , 


Data      Uoj  Uo 

»    (  m  /  s  )       (  W  /  m    i 


Ho 


■K  )       (W/m 


Qp 
■K  )  (W/m "2 


Tcf 
)  (C) 
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£-P  ;  '  "  ' 

!  '26 

ElSE 

i  7  2  8 

c  f 

PRINT  USING  '  10X," "Data    Vw       Uo           Ho 

172  3 

i  5 

-PINT  USING  '  10X,""  tt     !w/s>    (W/m-2-h      CW/m"2-K) 

p  ...  ii  H 

!  730 

1   U  z1 

END  IF 

!  740 

END  IF 

1747 

L  .' 

=  0 

1750 

z  ■• 

2  =  0 

1753 

- 

y'=0 

1755 

£■  y 

=  0 

1755 

S' 

=  0 

I7G2 

2  y 

=  0 

1  755 

3  ■ 

5=0 

1753 

3  * 

y=0 

1771 

Go 

_on=l 

1774 

Rep 

eat  :  ' 

1777 

J  = 

J  +  1 

1780 

IF 

I  «i=l  THEN 

17S3 

6EEP 

1785 

INPUT  "LIKE  TO  CHECK  NG  CONCENTRATION  (  1 =Y  ,0  =  N  )? "  ,Ng 

173  3 

IF  J=l  THEN 

1732 

OUTPUT  703; "AR  AF40  AL41  URS" 

1735 

OUTPUT  703; "AS  SA" 

1738 

END  IF 

1801 

6EEP 

1804 

INPUT  "ENTER  FLOWMETER  READING" ,Fm 

1807 

OUTPUT  709; "AR  AF60  AL62  URS" 

1510 

OUTPUT  709;  "AS  SA" 

1813 

ENTER  709;Etp 

IS  i  6 

OUTPUT  709; "AS  SA" 

1313 

BEEP 

1  c  22 

INPUT  "CONNECT  UOLTAGE  LINE" ,0k 

1325 

ENTER  709;6vol 

132  3 

BEEP 

1831 

INPUT  "DISCONNECT  UOLTAGE  LINE" ,0k 

1  334 

IF  6 vol  <.  1  THEN 

1837 

BEEP 

1  840 

BEEP 

1843 

INPUT  "INVALID  UOLTAGE,  TRY  AGAIN i ",0k 

'  346 

GOTO  1319 

1849 

END  IF 

1858 

OUTPUT  709;  "AS  3A" 

1361 

ENTER  709;Bamp 

1862 

Etp=Etp*1 .E+6 

Op 


U/m 
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IGE5  OUTPUT  709;"AR  AF40  AL47  UR5" 

1374  Nn=7 

1876  FOR  [=0  FO  Nn 

'  B79  OUTPUT  '709;  "AS  SA" 

1885  5e=0 

1883  FOR  K=l  TO  10 

1691  ENTER  709 ;E 

1394  5e=5e+E 

18  97  NEXT  K 

1 900  Enf (  I  >=ABS< Se/10  ) 

191 G  Emf < I  )=Emf< I  >*1  .E  +  6 

1918  NEXT  I 

192  1  OUTPUT  709; "AS  SA" 

1 924  OUTPUT  713; "T1R2E" 

192":  WAIT  : 

1 930  ENTER  7  I  3  ;  T  1  1 

1933  OUTPUT  -713;  "  T2P2E  ' 

1936  WAIT  2 

1939  ENTER  713;T2 

1342  OUTPUT  713;  "TIR2E" 

1945  WAIT  2 

1  948  ENTER  713;T12 

1 951  T1«(T1 1+T12  )*.5 

1 954  OUTPUT  713; "T3R2E" 

1960  BEEP 

1970  INPUT  "ENTER  PRESSURE  GAGE  READING  (Pga)",Pga 

197!  Pvapl=Pga*6894.7  I  PSI  TO  Pa 

1972  OUTPUT  709; "AR  AF64  AL64  UR5" 

1973  OUTPUT  709; "AS  SA"    I  PRESSURE  TRANSDUCER 

1974  Ss=0 

1975  FOR  K-1  TO  20 
1575  ENTER  709;Etran 
1977  5s=Ss+Etran 
1973  NEXT  K 

1979  Ptran  =  ABS(Ss/20  ) 

1980  BEEP 

193P  PRESSURE  IN  Pa  FROM  TRANSDUCER 

1 982  Pvap2  =  (  -2 . 93604 *Pt rant  14 . 7827  •♦6  394.7 

1985  ELSE 

1986  IF    lstu=0    THEN 

1989  ENTER    @F i le ; Bvo i  ,Bamp  ,Etp  ,Fm  ,  T 1  ,  T2  ,Pvap 1  ,Pvap2  ,Enf ( *  ) 

1990  ELSE 

1992  ENTER    @File;Bvol  ,Bamp .Utran ,E t p ,Emf ( 0 ) ,Emf (  1  )  ,Enf  (  2 ) ,Emf  (  3 )  ,Emf  (  4  > ,F 

m  ,T 1  ,T2  ,Phg  .Pwater 

1994  END    IF 

1996  IF    J=t'   OR    J  =  20    OR    J=Nrun    THEN 
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1997 
I  998 

1999 
2000 
2002 
2003 
2008 
2009 


2012 
2015 
2023 
2038 
2039 
2042 
2043 
2044 
2045 
204G 
2048 
2050 
2053 
2063 
2081 
2082 
2083 
2084 
2085 
2087 
2088 
2089 
2092 
2098 
2104 
2121 


-1 

124 

2 

125 

2 

126 

~i 

127 

n 

128 

n 

131 

-> 

134 

2 

137 

Ng=  1 

EL5E 

Ng  =  0 

END  IF 
END  IF 
IF  lstu=0  THEN 

Tsteam  1  =FNT  ,svS7(  Emf  (  0  >  ) 

T5teaml=T5teaml-273.  15 

Tsteam2=FNTvsv5G(Emf <  I  I 

Ts  team2  =  T5 1  eam2-273  .  15 

Ts  t  eam  =  l"5  t  earn  1 

Troom=FNTvsv53< Emf  (  2 )  ) 

Troom=Troom-273. '5 

Tcon=FNTvsv58( Emf ( 7  ' ) 

Tcon=Tcon-273. 15 
ELSE 

rsteam=FNTv5v(Ercf(0  )  ) 

Troom=FNTvs v( Emf ( 3  )  ) 

Tcon=FNTvsv  (  Erlf  (  4  )  ) 
END  IF 

F5at=FNPvst< Tsteam  ) 
Rohg=l 352  9-1 22 *( T room- 2 6 .85  )/50 
Rowat er  =  FNRhoui<  Troon  ) 
IF  lstu=0  THEN 

Ptest l=Pvapl 

Ptest2=Pvap2 
ELSE 

Ptest2="  Phg»Rohg-Puiater»Rouater  :*9.31  /  1  000 
END    IF 

PI-  s=Psat*1  .E-3 
Pi  p=Ptest2» 1  .E-3 
Pi  t=Pks 

Fsat=FNTvsp< Psat  ) 
Vst=FNVvst  <  Tsteam  ) 
Ppng=(Ptest2-Psat  )/Ptest2 
P p s  t  =  1 -Fpng 
Mwv=1 3. DIG 

IF  Ift-1  THEN  Mwv-137    i  TO  BE  CORRECTED 
IF  Ift=2  THEN  Mwv=62 
L'fng  =  (Pte5t2-Psat  ) /Ptest  2 
Mfng=l /(  1 +(  I /Vfng-1  >»Mwv/28.97 ) 
Mfng=Mfng* 100 
*Jfng  =  Vfng*100 
BEEP 
IF  Iov=2  THEN 

PRINT 
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2138' 

RECORD  TIME  OF  TAKING  DATA 

2139 

IF  Im=1  THEM 

2140 

OUTPUT  ^09; "TD" 

2141 

ENTER  709;Told$ 

2142 

END  IF 

21  44 

PRINT  USING  "i0X,""Data  set  nurocer            =  ""  ,DD,4X 

, 1 4A" ; J  ,To 

las 

2145 

OUTPUT  709;  'AR  AF40  AL40  UR5" 

2  1  46 

OUTPUT  709;  "AS  5A" 

2  149 

END  IF 

2152 

IF 

Iov=2  AND  Ng-1  THEN 

2  155 

PRINT  USING  "10X,""    Psat     Ptran    Tmeas     Tsat 

NG  '/."" 

" 

2159 

PRINT  USING  "10X,""  <A-Pa)         ( kPa )      (C)      (C) 

Molal  "" 

" 

2151 

PRINT  USING  "  1  IX  ,1  (3D.DD  ,2X  )  ,1  (  3D .  DD  ,3X  )  ,2(3D.DD,2X  )  ,2X  ,1 

( M3D.D  ,2X  i 

"  :P^ 

Pl-p  ,T 

5  t  earn  ,  T  5a t  ,M  f  ng 

2164 

PRINT 

2167 

END  IF 

2170 

IF 

Mfng  .5  THEN 

2173 

BEEP 

21^6 

IF  In-1  AND  Ng-1  THEN 

2179 

BEEP 

2132 

PRINT 

2  135 

PRINT  USING  " 1 0X  ,  " "Energi ze  the  .acuun  system 

2188 

EEEP 

2191 

INPUT  "OK  TO  ACCEPT  THIS  RUN  (  1 =Y  ,0=N )? "  ,0k 

2194 

IF  Ok=0  THEN 

2197 

BEEP 

2200 

DISP  'NOTE:  THIS  DATA  SET  'JILL  BE  DISCARDED'  '  " 

2203 

WAIT  5 

2206 

GOTO  1780 

2203 

END  IF 

">"?\n 

END  IF 

2215 

END  IF 

221  8 

IF 

Im=1  THEN 

2221 

IF  Fn\10  OR  Fm  100  THEN 

2224 

lfm=0 

->  n  n  — « 

BEEP 

2230 

INPUT  'INCORRECT  FM  (  1 =ACCEPT  ,0  =  DELETE  ) "  ,  I f m 

2233 

IF  lfm  =  0  THEN  1804 

END  IF 

2239 

END  IF 

2242' 

ANALYSIS  BEGINS 

2243 

IF 

lstu=0  THEN 

Ti l=FNTvsv58(Emf < 3  )  ) 

2262 

Ti2=FNTvsv55(Enf< 5  )  ) 

2272 

Tol =FNTvsv58( Enf ( 4  )  ) 

~i  -)  o  ~t 

To2  =  FNTvsv55(Enf(6  )  ) 
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::g: 


Z34Z 

Z353 

Z3S4 

Z0IE- 

2355 

Z353 

Z353 

Z3GI 

PILE 

Z3GZ 

Z364 

Z3G5 

Z367 

i.  JDO 

Z369 
Z370 
Z375 
Z377 
Z378 
2379 
2380 
2381 
Z33Z 
Z383 
Z3S4 
Z385 
2336 
Z387 
2338 
2338 
2331 
2292 
2393 
2394 
2395 
2396 
2397 
2398 
2399 
2400 


Ti  l=Ti ]  -273.  15 

r;Z  =  TiZ-Z73.  15 
Tol=Tci-Z'73.iS 
TcZ  =  ToZ-Z73.  15 
Tdel ! =T3l -T: 1 
Tdel2=To2-Ti2 
Tdei3=T2-Tl 
Etpl-Enf ( 3  )+Etp/20. 
Dtde  =  2  .5931E-2-I  .  5046  J-E-5  *  E  t  p  1  +  I  .2 
:  t  p  I  ■•  4 

Tri5=Dtde*Etp/10. 
Fo3  =  Ti 1 +Tr  i s 
IF  Iov  =  2  THEN 

PRINT  USING  "  IX  ,"  "   TIN1    TOUT 


'0IE-10*Etpl  "2 -5.  1  1  64E-  1  5*E tp  1  2^-Z 


TIN3 


T0UT3   DELTI    DELT3 


PRINT    USING    "IK,""  'TEFLON)  'QUARTZ)    

PRINT    USING    "2X  ,7(3D.0D,2X  ) "  ;  Ti  1  ,To1  ,T1  ,T2  ,Tdel  1  ,  Tdel  3  ,  Tns 

END    IF 

Erl =  ABS( Ti I  -T1  ) 

Er3=AES(Ti2-T1  ) 

PRINTER    IS    I 

BEEP 

Er2=ABS<  <  T2-T1  >-<  Tns  )  )/(  T2-T1  ) 

IF  Er2  .05  AND  Im=l  THEN 
BEEP 

PRINT  "QCT  AND  T-PILE  DIFFER  BY  MORE  THAN  57." 
Qk2=l 
IF  OK2=0  AND  Er2  .05  AND  Im=1  THEN  1780 

END  IF 

PRINTER  IS  701 
ELSE 

Ti=FNTvsv<  En f < 2  )  ) 

Grad=FNGrad<  < T 1 +T2  )* .5  ) 

To  =  Ti+ABS(Etp  )/(  10*Grad  >* I  .E  +  6 

Til=Ti 

T  o  3  =  T  o 
END  IF 
IF  lstu=0  AND  Itm-0  THEN 

TI i=Ti I 

T2o=To 1 
END  IF 
IF  Itn=1  THEN 

TI i=Tl 

T2o=T2 
END    IF 
IF    Itn-2    THEN 
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240 1  T1i=Til 

2402  T2o=Tq3 

2403  END  IF 

2  4  04  Tavg=(T1i  +  T2o)*.5 

2405  I  ft  =0 

2406  Cpw  =  FNCpw-  Tavg  ■ 
24^7  Rhow=FNRhow( Tavg  ) 
2408  IF  lstu=0  THEN 

24  10     Md=(6.7409*Fm+13.027)/ 1000. 

241  1      Md=Md*(  1  .0355-!  .95644E~3*T1  i+5.252E-6*T1  i"2  )/1  .0037 

2412  ELSE 

2413  Md=1 .04805E-2+5.30532E-3*Fn 

24  1 4     Md=Md»<  1  .0355-1  . 96644E-3*T1 i+5.252E-S»T1 i*2  >/ .995434 

2415  END  IF 

2417  Mf=Md/Rhow 

241 8  Uu=Mf / < PI*Oi '  2/4  ; 

2419  'Jw5  =  tJu»<  Di/ 1  .27E-2  )    2 
2  42  0  Tcor=FNTfnc<  T2o  ) 

2444  T2o=Tcor 

2445  Q=Md»Cpu*< T2o-T1 i  ) 
244G  Qp=Q/( PI*Do*L  ) 

2447  lft=0 

2448  Kw=FNKw(Tavg  ) 

2449  Muuj=FNMuuj(  Tavg  ) 

2450  Rei=Rhouj»l7w*Di  /Muu       '     ASSUMED    SAME    FOR    KORODENSE 
245!  Prw  =  FNPruj(  Tavg  ) 

2452  Fel=0. 

2453  Fe2=0. 

2454  Cf-1 . 

2455  Prwf=Pru 
245G  Reif-Rei 

2451  IF  Ihi-0  THEN 

2453     Ome=Pei"Pe>;p*PruT  .3333*Cf 

2465  END  IF 

2466  IF  Ihi-1  THEN 

2467  Sra=.88-<  .24/(4. +Prwf  )  ) 

246  3     Srb=. 333333-. 5 *EXP( -.6*Prwf  ) 

2470  0ne='5.  +  . 0  I 5»Rei f"Sra*Pruf "Srb  ) 

2471  END  IF 

2472  IF  Ihi=2  THEN 

247  3      Epsi  =  (  1  .82*LGT(  Rei  )-1  .64  )  ' ',  -2  ) 

2474  Ppk  1  =  1  .+3.4*Epsi 

2475  Ppk2=M  .7+1  .8*PruT<-1/3  ) 

2476  Ppl=(Epsi/8 >*Rei»Prw 

2477  Pp2  =  (Ppk !+Ppk2*<Epsi/8  )' .5»(Pru' .6666-1  )) 

2478  0me=Pp1/Pp2 
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12479  END    IF 

2461  Hi=Kui/Di»Ci»Ome 

:4g:  if  :►-  =  ,!  ihen      i9> 

2433  P  1  =  P  I  -  •'  Do  +  Qi 

248  4  A1=<  Ql-Di  MPI*' 01 +Di  *.5 

2465  M1=(Hi*Pl   Kcu*AI  )  '  .5 

2436  PZ=PI*« 0i+u2  i 

2487  A2=(D2-0i  >*PI*(Di+D2  ♦ . E 

2488  M2  =  <Hi»P2AKcu*A2  )  >"  .5 
2439  Fel=FNTanh(M1 *L!  )/<Ml *L1  ) 
2  4  90  Fe2=FNTanh(M2*L2  )/(M2*L2 ) 

2431  Dt=Q/(PI*Di*<L+L1*Fel+L2»Fe2 )*Hi  ) 

2492  IF  lhi=0  THEN 

2494  Cfc  =  <  Muu/FNMuw(  Tavg+Dt  !  )  '  .  1  4 

2495  IF  ABS( (Cfc-Cf  )/Cfc  )    .001  THEN 
2497        Cf  =  ( Cf  +  Cfc  )*.  5 

2500  GOTO  246  1 

2501  END  IF 

2503  END  IF 

2504  IF  Ihi-I  THEN 

2505  Pru)fc  =  FNPru<  Tavg+Dt  ) 

2506  Reifc=Vw«Di*FNRhow< Tavg+Dt )/FNMuw< Tavg+Dt  ) 

2507  IF  A6S(  (Pruifc-Pruf  )/Prwfc  )>  .001  OR  ABS(  (  Rei  f  c-Re  i  f  )/Re  i  f  c  ); 

2508  Pruif  =  ( Prwfc+Pruif  )/2  . 

2509  Reif=' Rei fc+Reif  )/2. 

2510  GOTO  2451 
251  1      END  IF 

2513  END  IF 

2516  Ift-Ifto 

2517  Lntd=(T2o-T1  l  )/L0G(  i  Tsteam-TI  i  )/<  Tsteafi-T2o  )  ) 
2513  Uo=Q/CLntd»PI*Do»L) 

251 9  Ho= 1/(1 /Do-Do *L/( D l * i  L  +  L 1 »Fe 1 +L2 *Fe2  ) *H l  >-Rn ) 

2520  Tcf=Qp/Ho 

252!  Cp5c  =  FNCpuj(/Tcon  +  T5team;*.5) 

2522  Hfg  =  FNHfg^ Tsteam  ) 

2524  Two=Tsteam-Qp  -'Ho 

2527  Tf ilm=Tstean/3+Tuo*2/3 

2530  Kf=FNKw(Tf iln ) 

2533  Rhof=FNRhow< Tf ilm ) 

2536  duf  =  FNMuw( Tf iln  ) 

2539  Hpq=.651 *hf  *<  Rhof "2*9.31 *Hf g/ ( Muf *Do*Qp  )  >" .3333 

254  1  Hnu5  5=.728*(Kf  '3*9.81  »Hf  g*Rhof  '2/  (  Muf  *Do*Tc  f  )  ) "  .  25 

2542  Alp1=.728*Ho/Hnuss 

2548  Tfm<  J-1  )=Tf iln 

2551  Qpa( J-1 )-Qp 

2554  Y=HDq*Qp ' .3333 


001  THEN 
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i  TO  BE  MODIFIED 


2557  <=Qp 

2560  Sx=Sx+X 

255  3  Sy  =  Syi-Y 

2555  S*s  =  S/.  s  +  X  2 

255  3  Sxy=S  '.y  +  X*Y 

2E72  01=502 

25  ~  5  Qlo55=01/ (  100-25  )*( Tateam-Troom 

2575  Hfc  =  FNHf < Tcon  ; 

2584  Mdv=0 

2587  Bp  =  <  Bvol* 100  I  '2/5.75 

2590  Hsc=Cpsc*( Tsteam-Tcon ) 

2593  Mdvc=< (Bp-Qloss  )-Mdv#Hsc  )/Hfg 

259b  IF  ABS( (Mdv-Mdvc  //Mdvc  )  .01  THEN 

2599      Mdv=< Mdv+Mdvc  )» .5 

2B02     GOTO  2595 

2605  END  IF 

2508  Mdv=i Mdv+Mdvc  i*.5 

25  1 1  Ug=FNUvst< Tsteam  > 

2514  Uv=Mdv»Ug/Ax 

2  520  F=i  9.81 *Do*Muf »Hf g  )/( VvA2*Kf *( Tstean-Tuio  )  ) 

2623  Nu=Ho*0o/Kf 

2525  Ret=Uv*Rhof«Do/Muf 

2529  Nr=Nu/Ret" .5 

2630  Hfuj  =  .96*(9.81  *Hfg/Tcf  )A  ,2*Kf  A  .8*VvA.  1  *Rhof  A.5/(Do*Muf  )A.2j 

2635  IF  Iov=2  THEN 

2645      PRINT 

2547      FRINT  USING  "EX.""  Uw       Rei  Hi  Uo  Hfuj(DT) 

Hnu<  Q  ) 

PRINT  USING  "5X  ,Z.DD,!X  ,3(MZ.3DE  ,1X  >  ,3X  ,2(MZ.3DE  ,3X  )"  ;Uui,Rei  ,Hi  ,Uo  ,Hfuj 


2550 
,Hpq 
2651 
2653 


PRINT 

PRINT  USING  "SX,""  Uv        Ho  q  Tcf       NuRe       F 

Hnu( DT  ) 

2G55      PRINT  USING  " SX  ,  Z . DD  ,  1 X  ,2( MZ . 3DE  ,  1 X  )  ,2X  ,3D . DO  ,2X  ,3( MZ . 3DE  ,  1 X  ) " ; Uv  ,Ho ,0p 
,Tcf  ,Nr  ,F  .Hnuss 
2656     PRINT 
2558   END  IF 

2659  IF  Iov=1  THEN 

2660  IF  Ihi-1  THEN 

2561      PRINT  USING  " 1 1 X  ,DD ,2X , Z . DD , 1 X ,3( MD . 3DE  ,  1 X  )  ,2 ( 3D . DD  ,  1 X  )  ,D . DDD" ; J  ,Uu  ,Uo , 
Ho  ,Qp  ,  Tc  f  ,  Ts t earn  ,Sra 
2662     ELSE 

2568    PRINT  USING  " 11 X  ,DD  ,4X  ,  Z . DD  ,2X  ,2( MD . 3DE  ,2X  )  ,  Z  .  3DE  ,3X  ,3D . DD  ,2X  ,3D . DD" ; J  ,Uu 
,Uo  ,Ho  ,Qp  ,Tcf  .Tsteam 
2571     END  IF 
2674   END  IF 
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1675 
:G84 

:594 
:707 

:_2  = 

!709 

:7i  i 


•  7  i  q 


[F  [m=2  THEN 

[F   IuuL=0  ANC  [job=1  )  OR  lull  0  THEN  OUTPUT  i?F  i  i  ep  ;  Op  ,Ho 
END  IF 
BEEF 
[F  Im=l  THEM 

[F   lwil=0  AND  Ijob=!  )  OP  I ui 1 1=  1  THEN  OUTPUT  @Filep;Qp,Ho 

INPUT  'CHANGE  TCCOL  RISE1?  I  =Y  ,  2=N",Itr 


l  ir- 


'HEN  SOTO  2384 


IG 

25 


2^52 

2755 

2753 

2  75  ! 

2"54 

2^57 

2^_0 

-i  -i—  -7 

2775 

2779 

27P2 

2735 

2'738 

279! 

2^97 

2S00 

2303 

2305 

2309 

231  2 

2  815 

23  13 

_  ~l  ._  *_ 

2524 

TQ-n-7 

i_  3*.  i 

2530 

2333 

2343' 

2353 

286  3 

2373 

2983 

BEEP 

INPUT  'OK  TO  5T0RE  THIS  DATA  SET  <  1 =Y  ,0=N  )?"  ,Qks 

IF  01  5=1  THEN 

OUTPUT  raFileiBvol  ,6amp  ,Etp  ,Fm  ,TI  ,T2  ,Pvap1  ,Pvap2  ,Enf ( *  ) 

AlpZ=Alpl+Aip2 
ELSE 

J  =  J  -  1 
ENO  IF 
BEEF 

:rlr'_"  'WILL  r^EFE  BE  ANOTHER  RUN  (  1  =Y  ,0=N  )? "  ,Go_on 
Nrun=J 
IF  Go__on   0  THEN  Repeat 

Else 

IF  J  Nrun  THEN  Repeat 
END  IF 

IF  Ijob=l  THEN  2809 
IF  Iwil=a  THEN 

ASSIGN  @File  T0  * 

I  j  o  D  =  I 

I  uj  d  =  1 

I  m  =  2 

ASSIGN  SFile  r0  0_fileS 

OOTO  I  1 36 
END  IF 
IF  Ifg=3  THEM 

PRINT 

Sl  =  v  Nrur.*S.<y-Sy»S.<  •  '  (  Nrun*Sxs-Sx  '2  ) 

Ac  =  <  5y-Si  *S.«.  >  Nrun 

FRINT  USING  "1 0X ,"" Leas t -Squares  Line  for  Ho  vs  q  curve:" 

PRINT  USING  "10X,""   Slope      =  " "  ,MD . 4DE" ; SI 

PRINT  USING  "10X,""   Intercept   =  ""  ,MQ . 4DE" ; Ac 
END  IF 
6EEP 

INPUT  'ENTER  SAME  TEMPORARY  PLOT  FILE  NAME",FplotS 
ASSIGN  §Fi lep  TO  P_f ile$ 
FOP  1=1  TO  Nrun 

ENTER  @Filep;Qp  ,Ho 

Xc  =  LOG'Op/Ho  ) 
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2884 

2885 

2386 

2  38  7 

2838 

Z389 

2690 

2381 

2892 

2893 

2394 

Z895 

2995 

2  897 

2898 

2399 

2902 

2905 

2908 

291  1 

2314 

2917 

2920 

2923 

2326 

2  92  9 

2930 

2931 

2  9w'  j 

2934 

2  339 

2  947 

2355 

2959 

29B2 

2955 

2958 

2971 

2974 

2977 

2390 

2989 

2995 

2998 

3001 

A"  ;  J 


.east -squares  line  for  q  -    a*de 1 t a-T "b " 
" "  .MZ.40E" ;Aa 
"'  .MZ.4DE"  ;Bb 


rc  =  L0G>  Op  ) 
Z-.  =  Z-  <-Xc 
Z  a  2  =  Z  *  2  +  X  c  2 
Z  .<.  y  =  Z  x  y  +  X  c  ♦  V  c 
Z/=Zy+Yc 
NEXT  I 
Bb=.75 

Aa=EXPv  (  Zy-Bb*Z.v  )/Nrun  ) 
PRINT 

PRINT  USING  " 1 0X  , 
PRINT  USING  " 12X  , "  "  a 
PRINT  USING  "  12X  ,"  "b 
IF  lft=0  THEN 
IF  lpc=0  THEN 
Qps=2.5E+5 

IF  lic  =  0  THEN  Hop=932G 

IF  Iic-1  THEN  Hop=10IG5*(  .01905/Do  )'. 3333: 
END  IF 

IF  Ipc=1  THEN 
Qps=7.SE+S 

IF  Iic-0  THEN  Hop=7176 

IF  Iic-1  THEN  Hop  =  75S9*(  .01 905/Do  )"  .33333 
END  IF 

Hos=Aa" i  1 /6b  )*Qps  ( <  6b-  1  )/Bb  ) 
IF  lpc  =  0  THEN  Aas  =  2.32EJ-4 
IF  Ipc«1  THEN  Aas=2.59E+4 
Alpsm=.876    '  SWENSEN  DATA 
IF  luiil=0  THEN  GOTO  2959 
Enrat=Alp2/Alpsm 


Enhancement  ratio  (Del-T)  =■ 


( V  =  0.4  m / s  ) 
33333   i  WAN  PETTEN 


Enr=Ho5/Hop 

PRINT 

FRINT  USING  " I2X 
ELSE 

PRINT 

IF  Ift-I  THEN 

Aas=2687.2   I  ZEBROWSKI 
Mas  =  2557.0»'>  .01  905/Do  )  ' 

END  IF 

IF  Ift-2  THEN  Aas  =  92S9.7*(  .01905/Do >  '  .3 

Edt=Aa/ Aas 

Eq=Edt"<  4/3  > 
END  IF 

6EEP 

PRINT 

PRINT  USING  "  10X  ," "NOTE 
,D_f ile$ 


:D.3D" ;Enrat 


(U  =  0.25  rn/s  ) 


,11,""    data  points  were  stored  in  file 
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3007  BEEP 

3013  PRINT 

3016  PRINT  USING  '  I  0X  , ''  NOTE :  "",ZZ,""  X-Y  pair-  were  5tcred  in  data  file  "" 

Pi "  ;  J  ,  P  1  o  t  $ 

3031  BEEC 

3073  ASSIGN  @F  i  ie  TO  ♦ 

3079  ASSIGN  ©Filep  TO  ♦ 

3030  PURGE  "DUMMY" 

00  IF  Iso=3  THEN  CALL  Wilson(Ci) 

16  END 

1 8  DEF  FNPvst( Tc  ) 

21  COM  /Fid/  I  ft  ,Istu 

3  4  D I M  K  (  8  ) 

37  IF  Ift-0  THEN 

30  DATA  -7.691 3  3456  4  ,-35.080236  96  ,-168. 1706546  ,64.2  3235504  ,-11 8. 96  462  25 
33     DATA  4. 1671  1732  ,20.9750676  ,1  .E9  ,6 
36     READ  K( *  ) 
33     T=<Tc  +  273.  15;/G47.3 
42      Sun=0 
45      POP  N=0  TO  4 
48        5um=Sum+K(N  )*( 1 -T  )"< N+1) 

ji 51      NEXT  N 

3154     Br-Sum/(  T*<  1+K(5)»<  l-T  >+K(B)*(  1  -T  )  "2  )  )-<  1-T  )/(K(7  )*<  l-T)A2+K(8)  ) 

57     Pr=EXP(Br) 
3160     P=22120000»Pr 

3163  END  IF 

66  IF  I-ft-1  THEN 

69     Tf=Tc*1 .8+32+459.6 

72     P=10-  ( 33.0655-4330. 99/ Tf -9. 26 35 *LGT( T f  >+2 . 0539E-3*Tf  ) 

75     P=P* 1 0 1 325/ 1 4 . 69G 

"9  END  IF 

31  IF  Ift=2  THEN 

84     A  =  9. 3946 85- 3066. 1 /( Tc  +  273. 15  ) 
87     P=133.32* !0'A 

90  END  IF 

93  RETURN  P 

96  FNEND 

99  DEF  FNHfg(T  ) 

3202  COM  /Fid/  Ift  ,Istu 

3205  IF  lft=0  THEN 

3208      Hfg=2477200-2450*' T-10) 

321 1  ENO  IF 

3214  IF  Ift-I  THEN 

3217     Tf=T»l.S+32 

3220  Hfg=7.0557857E+1 -Tf*< 4.838052E-2H  . 26 1 9048E-4*T f  ) 

3223  Hfg=Hfg*2326 
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3225  END  IF 

322  9  IF  I  ft =2  THEN 

3232  Tk=T+273.15 

3Z35     Hfg=l  . 35264E+6-T1  *< 6  .  38263E+2  +  Tk *. 747462  ) 

3236  END  IF 

324  1  RETURN  Hfg 

32  44  FNEND 

3247  DEF  FNMuw< T  ) 

3250  COM  /Fid/  I f t  ,Istu 

3253  IF  lft=0  THEN 
3256     A  =  247.3/(  T-H33.  15  ) 
3259     Mu=2.4E-5*l0  A 

3262  END  IF 

3265  IF  Iff  I  THEN 

3268     Mu=8.96  298l 9E-4-TM  1  . 1 094609E-S-T»S . 56682 9E-3  ) 

3271  END  IF 

3274  IF  Ift=2  THEN 
3277      Tk=l / ( T+273. 15 ) 
3230      Mu  =  EXP( -!  1  . 0179  +  Tk  »(  1  .  744E  +  3-T1  *<  2 . S0335E+5-Tk » I  .  I  256 1 E  +  8  )  )  ) 

3233  END  IF 
3286  RETURN  Mu 
3289  FNEND 

3292  DEF  FNVvst< Tt  ) 

3295  COM  /Fid/  Ift  ,Istu 

3298  IF  lft=0  THEN 

3301      P=FNPv5t ( Tt  ) 

3304      T=Tt+273.1S 

3307      X= I 500/T 

331 0     Fl-1/<  1+T*1 . E-4 » 

3313      F2  =  (  1-EXP<  -X  )  )"2.5*EXP(  X  )/X"  .5 

33! 6      B=.0015*F1 -. 000942 *F2- . 0004882 *X 

3319     K=2*P/ i 461  .52*T  ) 

3322     U  =  <  1  +  <  1  +2*B*K  )"  .5  >/K 

3325  END  IF 

3328  IF  I  ft  —  1  THEN 

3331      Tf=Tt»1 .3+32 

3334     'J=l3.9553S7-Tf  *(  .  161  27262-Tf»5.  17261  90E-4) 

3337     U  =  U/ IS .0  I  8 

3340  END  IF 

3343  IF  Ift-Z  THEN 

3346      Tk=Tt+273.1S 

3349      P=FNPvst( Tt  ) 

3352     U=133.95*Tk/P 

3355  END  IF 

3358  RETURN  U 

3361  FNEND' 
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3354  DEF  FN  iu  ' 

3367    3  M   :  i  d /  [ f  t  ,  I  s  t  u 

\:  \)  if   [ft«a  '  K  '• 

3  3"Z      ,  j-4.21  i20353-T*i  2.26B26E-3-T*(4.42361E-5  +  2.7UZ3E-7*T)  ; 

33"=  END  IF 

33~=  :r  Ift=l  FKEN 

3382     Cpui=9.2507273E-l  +  T*<  9  .  3400433E-4+ 1  . 7207792E-b*T  ) 

3385  END  IF 

3388  IF  I  ft =2  THEN 

3391      Tk=T+273.l5 

3394     Cpw=4.  1868*H  .6834E-2  +  U  *<  3 . 35083E-3-Tk*<  7  .  224E-G-Tk*7  .  G  1  748E-9  >  >  > 

3397  END  IF 

3400  RETURN  Cpu*1000 

3403  FNEND 

340G  DEF  FNRhouh T  ) 

3409  COM  /Fid/  If t , Istu 

3413  IF  lft=0  THEN 

34  15     Rg=399.52945+T*<  . 0 1 2S9-T*{ 5 . 4825 1 3E-3-T* 1  .2  34  147E-5  )  ; 

3418  END  IF 

342  1  IF  Ift  =  i  THEN 

3424     Ro=l  .6207479E+3-T*'  2 . 2 1 36346+T»2 . 3573291 E-3  ) 

3427  END  IF 

3430  IF  I  ft— 2  THEN 

3433     Tk=T+273. 15-338. 15 

3436      Uf=9.24848E-4+Tk*< 8 . Z73GE-7  +  Tk *  -  9 . 2  4  44E- 1 0  +  Tk *3 . 057E- 1 2 ) ) 

3439      Ro=l/Uf 

3442  END  IF 

3445  RETURN  Ro 

3448  FNEND 

3451  DEF  FNPrw< T  ) 

3454  Pr^  =  FNCpw( T  >»FNMuun  T  )/FNKw( T  ) 

3457  RETURN  Prw 

3460  FNEND 

3463  DEF  FNKw( T  ) 

3466  COM  /Fid/  Ift  ,Istu 

34G9  IF  lft=0  THEN 

3472      X=(T+273. 15>/273. 15 

3475      Kw  =  -.32247  +  X*<  2.8395-X*(  1  . 8007-X*(  . 52577- . 07344 *X  )  )  ) 

3478  END  IF 

3481  IF  Ift=l  THEN 

3484     Ku=8.2095238E-2-T*( 2 . 22 1 4286E-4+T*2 . 3809524E-8 ) 

3487  END  IF 

3490  IF    Ift=2    THEN 

3493  Tk=T+273.15 

3496  Kui=4.  1868E-4*(51  9.442+ .32092 0*Tk  ) 

3499  END    IF 
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3502  RETURN  Ku 

5505  FNEND 

3508  DEF  FNTanrv  « 

351  1  P=EXP( X  ) 

351  4  Q=EXP< -X  ) 

3517  Tanh=(P-Q)/(P+Q  ) 

35Z0  RETURN  Tanh 

3523  FNEND 

352G  DEF  FNTvsvOJ  > 

3529  COM  /Cc/  C(7  ) 

3532  T=C(0) 

3535  FOR  1=1  TO  7 

3538  T  =  T  +  C(  I  >*ITI 

3541  NEXT  I 

3547  RETURN  T 

3550  FNEND 

3553  DEF  FNHfi T  ) 

3556  COM  /Fid/  I  ft  ,Istu 

3559  IF  I f t =0  THEN 

3562  Hf  =  T*(4.203S49-T*(5.88l32E-4-T»4.55160317E-6  )  ) 

3565  END  IF 

3568  IF  Ift-1  THEN 

3571  Tf=T«1 .8+32 

3574  Hf=8 .2073S71 +Tf *( . 1 9467857+Tf*1 .3214286E-4) 

357"  Hf=Hf*2.326 

3580  END  IF 

3583  IF  I  ft =2  THEN 

3586  Hf=250   i  TO  BE  VERIFIED 

3589  END  IF 

3592  RETURN  Hf *  1000 

3595  FNEND 

3598  DEF  FNGrad( T  ) 

3601  Grad-37.9853+. 104388*T 

3604  RETURN  Grad 

3607  FNEND 

3610  DEF  FNTvspCP > 

3613  Tu=190 

3616  Tl=10 

361  9  Ta  =  ( Tu  +  Tl  )*.B 

3622  Pc  =  FNPvst( Ta  ) 

3625  IF  ABS< ( P-Pc )/P  ).   .0001  THEN 

3628  IF  Pc<P  THEN  Tl=Ta 

3631  IF  Pc;P  THEN  Tu=Ta 

3634  GOTO  3619 

3637  END  IF 

3540  RETURN  Ta 
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5B43  FNEND 

■S646  DEF    FNSigma    : 

66*4  9  <=547.3-T-273.iS 

S65Z  5-.  11 60935507   I  +  . S3*X  )+ I  .  1  2  1  404688E-3-5 . 752305  I 3E-5 *X  +  I  . 236Z746SE-3* <  "2-  I 

. 14971 9Z9E-1 1 *X' 3 

6655  RETURN  S* .001 *x    2 

5653  FNEND 

6325  SUB  Wi  lsom  Ci  ) 

63Z3  COM  /  U)  1 1  /  Mrun  ,  I  tm  ,  lath  ,  I  roc  ,  I  f  e  ,  I  job  ,  Iwd  ,  I  f  g  ,  Ipco  ,  I  f  to  ,  Iwi  1  , 1  hi  ,  loc  ,  Inam  ,K 

:u  ,Rexp  ,Rm  ,  Ax 

633Z  COM  /Fid/  If t  ,Istu 

6333  COM  . Geom/  D 1  ,02  ,Di  ,Do  ,L  ,L 1  ,L2 

6334  COM  /Fnc/  I  5  t  uo  ,  Inn  ,  Iut  ho  ,Uw  ,(Jw5  , 1  nano 

6336  DIM  Emf  (20)  ,Bvo(42  ) ,Bam( 4Z  )  ,Eata(42  )  ,Ear( 42  ,7  )  ,Fma< 42  )  ,T1 a( 42  )  ,T2a( 42  ) 

5345  IF  Ioc-0  THEN 

6847  PRINT  USING  " 1 6X  , " "Nussel t  theory  is  used  for  Ho 

6848  ELSE 

6349      PRINT  USING  "  1 6X  , " "Fu j 1 1  correlation  used  for  Ho 

5850  END  IF 

6853  BEEF 

6856  INPUT  "RE-ENTER  DATA  FILE  BEING  PROCESSED"  ,D_f i le$ 

6859  BEEP 

686Z  INPUT  "GIVE  A  NAME  FOR  XY  PLOT-DATA  FILE",Plot$ 

63B5  CREATE  BOAT  Plots  ,  10 

6863  ASSIGN  @Io_path  TO  Plots 

6371  Jj=0 

6874  ASSIGN  @File  TO  D_fiieS 

6577  ENTER  ©File; I fg.Inn 

6373  IF  lstu=0  THEN 

6383     ENTER  §F i le ; Odd  ,0dd  ,Ddd  ,Ddd 

6364  ELSE 

6535      IF  Ifg-0  THEN  ENTER  ©File;  Iut 
•6666      IF  [fg=1  THEN  ENTER  £F  i  le  ;  Fp  ,Fui  ,Fh 

6887  END  IF 

6388  IF  J  j  =0  THEN 
6395      IF  lhi=0  THEN  Ci=.0Z7 

6896  IF  Ihi=1  THEN  Ci-1 .00 

6897  IF  Ihi-2  THEN  Ci-1 .00 

6899  IF  Ifg-0  THEN  Alp-1 .2 

6900  IF  Ifg-1  THEN  Alp=2.6 

6901  IF  Ift  =  2  AND  Ifg-1  THEN  Alp=5.0 
6904  END  IF 

6907  J=0 

6910  Sx=0 

6913  Sy=0 

6916  5x5=0 
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691  9 

E92Z  i 

59Z5 

6925 

55Zi 

593: 

6934 

m.TI  , 

6936 

5938 

6939 

6940 

6  94  3 

6946 

6949 

6952 

6955 

6956 

6953 

5959 

6960 

636  1 

696Z 

6964 

6967 

6970 

6373 

5976 

6  97  9 

698: 

6985 

5983 

699! 

699: 

6994 

6995 

5996 

6397 

6998 

7000 

7003 

7004 

7006 

7007 

7009 

7010 


PEfiD  DATA  FROM  n  USER-SPECIFIED  FILE 
IF  jj=0  THEN 

IF  Istu=0  THEN 

ENTER  i5Fi!e;Bvol  ,6  amp  ,Etp  ,Fm  ,T 
ELSE 

ENTER  @Fi  ie;EvoI  ,5 amp  ,'Jtran  ,E t p  ,Emf  (  0 
:  ,Phg  , P water 
END  IF 
5  v  o  (  J  >  =  B  v  o  1 
Ban1.  J  )=6amp 
Eata( J  )=Etp 
Ear<  I  ,0  >=Emf  <  0  ) 
1  )=Emf (  1  > 
Z   ' = E  m  f  (  2  ) 
3 >=Emf < 3  ) 
J  ,4  >  =  Emf  (  4  ) 
stu=l  THEN  GOTO  636  I 
J  ,5  ;=Emf  <  5  ) 
J  ,5  )=Emf (6  ) 
J  ,7  >=Emf  (7  ) 
J  )=Fm 


IF  INPUT  MODt 


im  i 


T:  ,Ddd  ,Ddd  ,Emf 


.Emf 


,  Em  f  (  ?.  )  ,Em  f  (  3  )  ,Em  f  (  4  i  , c 


Ear(  J 

Ear^  J 

Ear>  J 

Ear 

IF 

Ear 

Ear 

Ear< 

Fma( 

T1a( J  )=T1 

T2a<  J  >=T2 
ELSE 

6voI=B^o( J  ) 

Bamm=6am( J  ) 

Etp=Eata( J  ) 

Emf  < 0  )=Ear< J  ,0  ) 

Emf  (  1  )=Ear<  J  ,1  ) 

Emf(2  )=Ear<  J  ,2) 

Emf  (  3  )  =  Ean  J  ,3  ) 

Emf ( 4  >=Ear<  J ,4  ) 

IF  Istu-I  THEN  GOTO  6397 

Emf  (5  »=Ear<  J  ,5) 

Emf  i  6  >  =  Eari  J  ,6  ) 

Emf < 7 )=Ear( J  ,7  ) 

Fm  =  Fma( J  ) 

TI=T1a( J  ) 

T:  =  T:a( J  ) 
END  IF 
IF  Istu-0  THEN 

Tsat=FNTvsv57<  (Emf ( 0  HEmf  ( 

Tsat=Tsat-273  .  1 5 

Ti=FNTv5v5S(Emf ( 3  )  ) 

Tol=FNTv5v58(Emf (4  )  ) 


)/:.  ) 
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7ai : 

7013 

7015 

7015 

201E- 

7017 

~0!  3 

701  g 

"020 

702  1 

7022 

7023 

7024 

7025' 

"02  G 

702  7 

7028 

702  9 

7020 

7031 

7032 

7033 

7034 

7035 

7035 

7035 

7039 

7040 

^041 

7042 

7043 

7044 

"04  9 

7048 

705  1 

7054 

-055 

7057 

7050 

706  1 

"705  2 

70S  3 

7065 

7066 

7067 

7068 

T  :  =  t  L -273 .  '5 

r :,  i  -  t  o  1-273.  15 

£  *  p  I  =  E  m  f  (  2    •-  E  t  p    2  0  . 
•-e  =  2.553iE-:-:  .  50464E-S*EtpH-1  .  2  1  70  1  E-  !  3*E  tp  1     2-5.1  !64E-15*Etp 
!  9  *  E  t  p  I    4 

Tr  I  5  =  0 1  de*Etp/ 1 0  . 

To=T ;+Tr i s 
ELSE 

Tsat=FNTv5v< En f ( 0 i  ) 

Ti=FNTvsv( Enf < 2  )  ) 

Grad=FNGrad<  ! 7  I +■ T;  i . . c  , 

fo»Ti+ABS(Etp  ■/'  10*Grad  )*!  .  E  +  6 
END    IF 

CALCULATE    THE    LOG-MEAN-TEMPEPATURE    DIFFERENCE 
EF    lstu  =  0    AND    ltm  =  0    THEN 

F  f  =  T  i 

Ti=Tol 
END    IF 
IF    1 1 m= 1     "HEN 

T  f = T 1 

T1  =  T2 
END  IF 
IF  Itm=2  THEN 

n=To 

END  IF 

Tavg=( Tf  +  Tl  i* .5 

Trise=Tl-Tf 

Lmtd=Trise/LQG<  (  Tsat-Tf  )/'  Tsat-Tl  )  ) 

If  t=0 

Cpw=FNCpui(  Tavg  i 

Rhouj  =  FNRhciAj(.  Tavg  ) 

-  k^FM  w  Tavg 

Muuja  =  FNMuu(  Tavg  ) 

prw=FNPrui(  Tavg  ) 

r  f t=r  f to 

IF  lstu=0  THEN 

Mdt  =  (6.7409»Fri+13.02'7>/l  000  . 

Md=Mdt*' 1 . 0365-Tf ♦( I . 96644E-3-Tf *5 . 252E-6 ) )/ 1 .0037 
ELSE 

Mdt = 1 .04805E-2+6.80932E-3*Fn 

Md  =  Ndt»(  1 .0365-Tf *(  1  .96644E-3-Tf »5.252E-5  )  )/. 995434 
END  IF 
Uf=Md/Rhow 
v"uj=Vf/<  PI*Di'  2/4  ) 
'Juts=Uw*(Di/l  .27E-2  )'  2 
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"7069 

7070 

"07! 

"07: 

"074 

"076 

"078 

7080 

70s: 

7084 

"085 

708G 

7088 

7080 

708: 

7094 

7096 

7098 

7100 

7108 

7  111 

rill 

71  1  4 

71  17 

7i:0 

7123 

712G 

7i:7 

7i:8 

7i:9 

7151 

7132 

7135 

7133 

7141 

7144 

7147' 

"1  48 

7150 

7153' 

7154 

7156 

7159 

7162 

7165 

7168 

7171 

■Q    BE    MODIFIED 


T  c  -j  r  =  F  N  T  f  r  1  c  >.  T  r  ;  3  e  ) 

rnse=Tcor 

Tcon=FNTv5v53< Enf « 7  '  ) 

Tcon=Tcon-273.  15 

Cpsc=FNCpui  ■■  T  :cn-*-Tsat  .;  *  .  5  ) 

Hf g=FNHfg! Tsat  1 

01=500 

Qloss=Ql/( 100-25  )*(Tsat-Troom ) 

Mdv=0 

Bp  =  < Bvol* 100  )"2/5.7S 

Hsc=Cpsc*( Tsat-Tcon  ) 

Mdvc  =  <  '.  Bp-Qloss  )-Mdv*Hsc  1 '  Hf  g 

IF  ABS(  (Mdv-Mdvc  )/Mdvc  )'-■  .01  THEN 

Mdv=( Mdv+Mdvc  )« .5 

GOTO  7086 
END  IF 

Mdv  =  (  Mdv  +  Mdvc  <*  .  5 
Vg=FNVvst<  Tsat  ) 
Vv=Mdv*Vg/i°ix 
0  =  Md*Cpuj*Tr  1  52 
Qp=Q/(PI»Do*L  ) 
Uo=Qp/ Lnt d 
Re  =  Rhou*tJw*Di/Muuja 
Fel  =0 
Fe:  =  0 
Cf-1  . 
Prwf  =Prui 
Reif=Re 
Ift-0 

Tuo=Tsat -5 

Tfiln=T5at/3+Tuio*2/3 
h  f=FNKui<  Tfiln  > 
Rhof =FNRhou( T  f i im  ) 
Muf=FNMuw< Tf 1 im  ) 

Hfgp=FNHfg<  Tsat  >+ . 6S*FNCpw< Tf llm  >*<  Tsat -Two  ) 
Neui=Kf  *<  Rhof  "3*9. 81  *Hfgp/(  Muf *Do*Qp  )  )  '  .3333 
New=( Kf "3*9.81 *Hfgp*Rhof "2/ (Muf *Do*<  Tsat -Two  )  ) ) 
IF    Ioc-1     THEN 

New  =  Kf *( ( 9.8  1 *Hfgp/Qp  )     .35  )♦<  < Muf *Do  )" 

New=< 9.81 *Hfgp/( Tsat -Two  )  )A ,2*KfA .8*Vv 
END    IF 
Ho  =  Pilp  *Neio 
Tuioc=Tsat-Qp/Ho 
IF    ABS( (Twoc-Two)/Twoc  )> .001    THEN 

Two  =  Tujoc 

GOTO    7132 


25 


-.375  )  )*<  Rhof A .625  >*(Uv 
.  1 *Rhof " .5/(Do*Muf )" .3 


!5) 
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M  '4 
71  75 
'184 

7186 
7187 

71 3  a 

7189 
71  90 
7191 
7192 
7193 
7194 
7  195 
719G 
7197 
7198 
7199 
"7200 


7204 
7205 
7206 

7Z07 
7208 
7209 
7210 
7213 
72  16 
72  17 
721  9 
7222 
7225 
7223 
7231 

"^23^ 
7235 
7236 
7237 
7239 
7240 
7241 
7242 
7243 
7245 


END  IF 
Rexpi =  Re  vp 
Ic  [hi=0  T  ••  ■ 

m e g a = R e  Re.<pi»Pruj  .  3 3 3 3 * C  f 
END  IF 
IF  Ihi=1  THEN 

Sra=.88-(  .24  (4. fPrwf  )  ) 

Srb=.  333333+.  5*  EXP  <  -S*Pruif  ) 

Or,iega  =  <5.  +  .01S*Peif"Sra*Prujf  ' Srb  ) 
END  IF 
IF  Ihi=2  THEN 

Epsi  =  (  I  .S2»LGT( Re  )-l  .64  >"  < -2  ) 

Ppk  ! =1 . +3.4*Epsi 

Ppk2=1  1.7+1  .8*PruT  (  -1  /3  ) 

Pp  I  =  <  Ep  s  l  /  3  /*Re*F^'j 

Fp2  =  (Ppk  I  +Ppk  2»(Ep5i/3  )■  .S*'PruT  .6666-1  )  ) 

Omega=Fp  I  .Pp2 
END  IF 

H  i  =  K u / D l *  C i*0  mega 
IF  lfe  =  0  THEN  -72  16 
P1«PI*<  Di+Dl  ) 
P2=PI*( Di+D2  ) 
A1=«D1-Di  >*PI»(Di+D1  >*.5 
iqi2  =  <  D2-Di  »*PI  *(Di+02  >*.S 
M1=(Hi*Pl/(Kcu*n1  )  )-  .5 
M2  =  <  Hi  »P2.'  <  Kcu*A2  >  )'  .5 
Fe  i  =FNTanh( Ml *L1  »/(Ml *L1  ) 
Fe2=FNTanh(M2*L2  )/(M2*L2  ) 
Dt=Q/(PI*Dx*(L+Ll *Fe1 +L2*Fe2  )*Hi  ) 
IF  lhi=0  THEN 

Muwi=FNMuui<  Tavg+Dt  ) 

Cfc  =  <  Muwa/Muuii  )  '  .  1  4 

IF  n63>  ( Cfc-Cf  >/Cfc  i  .001  THEN 
Cf  =  < Cf+Cfc  >*.5 
GOTO  7  184 

END  IF 
END  IF 
IF  Ihi-1  THEN 

Prufc=FNPrw( Tavg+Dt  ) 

Reifc«Vui»Di*FNRhow<  Tavg  +  Dt  )/FNMuw(  Tavg  +  Dt  ) 

IF  ABS(  (Prwfc-Prwf  )/Prufc  )>  .001  OR  ABS(  (  Rei  f  c-Rei  f  )/Rei  f  c  >) 
Pruf  =  < Prw  f c+Pruf  )/2 . 
Rei  f  =  ( Rei fc+Reif  )/2. 
GOTO  7184 

END  IF 
END  IF 


001  THEN 
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72  46  I  f t  =  1 f to 

7247  X=Qo*New*L/< Omega»Ku*< L+L I *Fel +L2*Fe2 ) ) 

7243  r=New*<  I /Uo-Rn  ) 

72491  COMPUTE  COEFFICIENTS  FOR  THE  LEAST-SQUARES-FIT  STRAIGHT  LINE 

~250  IF  Jp=!  THEN  OUTPUT  @Io_path;X,Y 

•"'252  5*=Sa+X 

7255  Sy-Sy+t1 

'258  S  x  5  -  S  x  s  +  X  *  X 

""25  1  S*y=Sxy+X*Y 

72G4  IF  Im=1  nND  Jj=0  THEN  OUTPUT  @Fi le ; Bvol  ,Bamp  ,Etp  ,Fm  ,T1  ,T2  ,Pvap I  ,Pvap2  ,Eni 
*  ) 

7267  J=J+I 

7273  IF  J  Nrun  THEN  6925 

7273  3  1  =  <  Nrun*Sxy-Sy  *Sx  )/(  Nrun*Sxs-Sx  A2  ) 

7276  IF  Iwil-2  THEN 
?237      IF  lhi=0  THEN  51=1/ .027 
9      IF  I  n i  =  1  THEN  51  =  1  '  I  . 00 
7289      IF  I  hi =2  THEN  Sl-I / 1  .00 

72  91  END  IF 

7294  Ac=( Sy-Sl*Sx )/Nrun 

7297  Cic-1/Sl 

7300  Alpc=1/Ac 

7303  Jj=Jj+1 

7306  IF  Jp-|  THEN  Jp=2 

7309  Cerr=ABS( (Cic-Ci  )/Cic  ; 

7312  Aerr=ABS( ( Alpc-Alp  )/Alpc  ) 
7315  IF  Cerr  .001  OR  Aerr  .001  THEN 

731 3  Ci  =  (Cic+Ci  )*.5 
7321      Alp=(Alpc+Alp)».5 
7324      BEEP 

7327      IF  Ijob-1  THEN  6907 

"330  ELSE 

7333      IF  Jp=0  THEN  Jp=1 

"336  END  IF 

7339  IF  Jp=1  THEN  6874 

7342  Ci  =  < Ci+Cic  >*  .  5 

7345  PRINT 

7346  IF  lhi=0  THEN 

7348  PRINT  USING  "10X,""Ci  (based  on  Sieder-Tate)        =   "",Z.4D";Cl 

7349  END  IF 

7350  IF  Ihi-1  THEN 

7351  PRINT  USING  "10X,""Ci  (based  on  S le icher-Rouse )     =   "",Z.4D";Ci 
73521     PRINT  USING  "10X,""Re  exponent  for  S le  icher-Rouse   =   " "  ,D . ODD " ; Sra 

7353  END  IF 

7354  IF  Ihi-2  THEN 

7355  PRINT  USING  '10X,""Ci  (based  on  Petu^ hov-Popov  )     =   "",Z.4D";Ci 
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~Z5G  END  IF 

■357  IF  L3c=0  THEN 

7353      h-I'it    [N   '  IQX,   Alpha  'based  en  Nusselt  (  Tdel  >  !  =   '"  ,  Z  .  40 "  ;  A  Ip 

^359  E  ID  IF 

";Z50  IF  Ioc=!  'HEM 

73S1      -~":f."  USING   i0X,""Alpha   based  on  Fuj  1 1  (Tdel))    =   ""JZ.4D";ALp 

_Z5  2  END  IF 

"Z53  IF  Inam=5  OR  Inam=6  THEN 

"54      IF  lhi=0  THEN 

73EG         IF  lpro=0  HND  lnn  =  0  THEN  Alpsm=.8218   'NO  INSERT  .VACUUM  ,S-T 

7357  IF  Ipco=l  AND  lnn  =  0  THEN  Alpsm=.7793   i NO  INSERT  .ATMOSPHERIC  ,S-T 

7358  IF  lpco=0  AND  Inn  =  3  THEN  Alpsm=.7854   i HEATEX  .VACUUM  ,S-T 
73G9         IF  Ipco=l  AND  Inn  =  3  THEN  Alpsm=.77G9   i HEATEX  .ATMOSPHER IC  ,  3-T 
7371      END  IF 

"373      IF  Ihi=!  THEN 

7Z7Z         IF  [pco=0  AND  lnn  =  0  rHEN  Alpsm=.8613   'NO  INSERT  .VACUUM  ,S-R 

7374  IF  Ipco=l  AND  lnn  =  0  THEN  A lpsm= . 82 1 8   'NO  INSERT  .ATMOSPHER  I C  ,  3  -R 

7375  Ic  lpco=0  AND  Inn  =  3  THEN  Alpsm=.7791   ' HEATEX  .VACUUM  ,  S-R 
737G         IF  Ipco-1  AND  Inn  =  3  THEN  AIpsm=.7929   i HEATEX  .ATMOSPHERIC  , S-R 
■7-3" 3      END  IF 

7379  IF  I  h i  =  2  THEN 

7380  IF  lpco=0  AND  lnn  =  0  THEN  Alpsm=.8205   'NO  INSERT  .VACUUM ,P-P 

7381  IF  Ipco=1  AND  lnn  =  0  THEN  Alpsri=.75S4   'NO  INSERT  .ATMOSPHER IC  ,P-P 

7382  IF  lpco=0  AND  Inn  =  3  THEN  Alp5m=.7670   ' HEATEX  .VACUUM  ,P-P 
7333         IF  Ipco=1  AND  Inn  =  3  THEN  Alpsm= . 7708   ' HEATEX  .ATMOSPHERIC  ,P-P 
7385      END  IF 

738b  END  IF 

7387  IF  Inan  =  4  THEN 

90  IF  Ipco=l  THEN  Alpsm=.87G   ' SUENSEN  DATA  BASED  ON  DEL-T 

91  END  IF 

32  IF  lnam=0  OR  Inam=3  THEN 
93      IF  lpcc=0  THEN  AIpsm=.S3    ' VP  MISTV103 

'394      IF  Ipco=l  THEN  Alp5m=.S8    '  VP  SMTHSTAE5 

37      IF  Ift=l  THEN  Aipsm=.677    'VAN  PETTEN  (V  =  0.25  m/s  ) 
7399      IF  I  ft  =  2  THEN  Alpsm-1.2G2 

99  END  IF 

7401  IF  Inan-1  THEN   'MITROU  ALPHA  FOR  P-P  FROM  REPROCESSING 
"402      IF  lpco=0  THEN  Alpsm=.8437 

7403  IF  Ipco=1  THEN  AIpsm=.8418 

7404  END  IF 

7405  Et=Alp/AlP5m 
"406  Eq  =  Ef  1  .333333 

7407  PRINT  USING  " 1 0X  ,"" Enhancement  (q)  =  "",DD.3D";Eq 

7408  PRINT  USING  " 1 0X  ,"" Enhancement  (Del-T)  =  "",DD.3D";Et 

7409  ASSIGN  @File  TO  * 

7410  SUBEND 
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7721  DEF  FNTvsv55( U  ) 

^731  COM  /Cc55/  T55<5  ) 

7741  T  =  T55<0) 

7751  FOR  1=1  T0  5 

7761      ~=T*T5S;  I  )»W  I 

7771      NEXT  I 

7781  RETURN  T 

7791  FNEND 

7801  DEF  FNTvsv5G(V  ) 

781  1  COM  /Cc56/  T56<5  > 

7321  T=TS6(0) 

7831  FOR  1=1  TO  5 

7841      T=T+T56<  I  )*V  I 

7851  NEXT  I 

7861  RETURN  T 

7871  FNEND 

788  1  DEF  FNTvsvS7(U > 

7891  CON  /Cc57,'  T5?»  5  ) 

7901  T=T57(0) 

791 1  FOR  1=1  TO  5 

79Z1      T=T+T57(  I  )*VAI 

793  1  NEXT  I 

7941  RETURN  T 

7951  FNEND 

^96  1  DEF  FNTvsv58< U  ) 

7971  COM  /Cc58/  T58(S  ) 

7981  T=T58(0) 

7991  FOR  1=1  TO  5 

8001      T=T+T53^  I  )*U"I 

801 1  NEXT  I 

8021  RETURN  T 

3031  FNEND 

8100  DEF  FNTfncl Tcor  ) 

8110  COM   Fnc'  I  5  t  uo  ,  I  nn  ,  I  wt  ho  ,Vw  ,'Jws  ,  I  namo 

8115  IF  lstuo=0  AND  luitho  =  0  THEN  '  SWENSEN  FRICTION  FAC.  FOR  COPPER  TUBE 

8120      IF  lnn  =  0  AND  Uw>.S  THEN  Tcor*Tcor-( -2 . 73E-4+1  . 75E-4*Uw+9 . 35E-4*VuT2-l  .  9 

6  E  -  5  *  U  w '3 ) 

8130  IF    Inn-1     THEN    Tcor  =  Tcor-( -6  .  44E-5+ 1  .  7  1  E-3*Uw  +  4  .  45E-4#»Juj  ~  2  +  4  .  07E-5»Uw '3  ) 

8140  IF    Inn=Z    THEN    Tcor=Tcor-( -3 . 99E-4+2 . 75E-3*Vui  + 1 . 4SE-3«UuT2  +  8 . 1 6E-5*VwA3 ) 

8150  IF    Inn=3   THEN   Tcor=Tcor-< 8 .57E-5+1 . 23E-3*Uw+1 .08E-3*VuT2+8. 1 6E-5*VuA3 ) 

8160  END    IF 

8170  IF  lstuo=0  AND  Iwtho  =  2  THEN  iQKEEFE  FRIC.  FAC .  FOR  SMOOTH  TITANIUM  TUBE. 

8180      IF  lnn  =  0  AND  Vw>.5  THEN  Tcor  =  Tcor-< -4 . 62E-5-7 . 53E-4»Vw+1 . 80E-3»VuT2-S  .  8 

4E-5*Ww  3  ) 

8190  IF    Inn  =  3    AND    Inamo=5    THEN    Tcor  =  Tcor-(  2  .  09E-4  +  9  .  74E-4»Vu+2  .  1  2E-3*VuT2-3  . 

31E~5»Uw  '3  ) 


120 


8191  IF  Inn  =  3  AND  Inano=5  THEN  Tcor=Tcor-< 1 . 9555E-4+3 . 972 1 E-3*Vw+3 . 1 27E-4*Uw 
2+3.51 9E-4*Vw  Z  > 

3200  END  IF 

8210  IF  lstuo=0  AND  Iutho-1  THEN  'FRICTION  FACTORS  FOR  LFD  K0R0DEN5E 

5220  IF  lnn  =  0  AND  \)u    .5  THEN  Tcor  =  Tcor-<  -3  .  386E-4+  1  .  88E-3**Jw+6  .0  I  3E-4.«Vuj  >4 
.  I33E-5*V/u  3  ' 

5210  IF  Inn=3  THEM  Tcor=Tcor-( 2 .089E-4+9 . 202E-4*Vw+1 . 893E-3»VuT2-2 . 781 E-S»Vw 
'3  ) 

8240  END  IF 

8242  IF  l5tuo=0  AND  Iutho=3  THEN  'FRICTION  FACTORS  FOR  MHT  KORODENSE 

6244  IF  :nn  =  0  AND  Vu;  .S    THEN  Tcor  =  Tcor-( 6 . 79E-5+ 1  . 632E-3*Vw  +  8 . 409E-4»Uw  '  2+ 1  . 
1  1 lE-4*Vw  3  ) 

8245  IF  Inn  =  3  THEN  TCor  =  Tcor-  (  2  .  5S4E-4+6  .  263E-4*Uu+2  .  S03E-3  *Uuj  "  2  +  7  .  330E-S  »Uw 

3) 

3247  END  IF 

5250  IF  [atuo-l  THEN 

8250'  IF  lnn  =  0  THEN  Tcor=Tcor-(  . 01 38+ . 001 »VuT2  ) 

3270  IF  Inn=1  THEN  Tcar  =  Tcor- .  004*Uujs  "2 

8280  IF  Inn  =  2  THEN  Tcor-Tcor-.004»VwsA2 

S290  END  IF 

8300  RETURN  Tcor 

8310  FNEND 
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